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Abstract 

Traditionally,  logistics  analysts  have  treated  the  delivery  problem  separate  from 
the  inventory  and  repair  problems.  In  other  words,  each  of  these  three  problems — 
delivery,  inventory,  and  repair — are  treated  individually.  We  combine  the  vehicle  routing 
problems  and  inventory  allocation  problems  in  a  single  formulation.  Furthermore,  lateral 
re-supply  is  considered,  in  that  supplies  are  not  only  available  from  the  main  depot,  but 
also  from  other  bases.  After  we  get  the  lateral  re-supply  information,  we  set  a  delivery 
allocation  plan  out  of  the  main  depot.  Such  allocations  are  then  delivered  via  vehicles 
stationed  at  the  main  depot. 

There  is  a  repair  capability  for  each  base  other  than  the  main  depot  as  well.  There 
is  a  newsboy  inventory-cost  function  associated  with  each  base  to  the  every  other  base. 
We  set  up  a  model  for  repairable  items  only.  Each  base  has  two  options:  accept  the 
repairable  items,  or  deliver  the  repairable  items,  but  not  both.  If  the  base  chooses  the 
“accept”  option,  we  make  sure  that  that  base  has  used  all  its  available  resources,  initial 
inventoried  items  and  repair  capabilities,  before  it  receives  re-supplies.  In  the 
transportation/delivery  submodel,  we  place  limitation  on  the  crew  duty  hours  available. 
The  objective  function  is  to  satisfy  the  demands  by  minimizing  traveling  cost  and 
inventory  cost. 

We  solve  the  same  problem  by  using  a  generalized  Benders  decomposition 
technique.  The  decomposition  allows  us  to  attack  larger  problems  and  use  general  failure 
and  repair  functions  familiar  to  logistics  analysts.  Computational  experience  suggests 


xi 


that  efficiency  is  achieved  by  combining  the  delivery  and  inventory  functions.  Lateral  re¬ 
supply  and  computational  efficiency  is  particularly  useful  in  emergency  situations. 
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A  COMBINED  INVENTORY  AND  DELIVERY 
MODEL  FOR  REPAIRABLE  ITEMS 


Chapter  1-Problem  Background 


Overview 

Is  logistics  a  one-time  a  life  thing?  Is  it  ending?  Is  it  crucial?  Is  it  possible  to 
think  of  a  private  or  public  enterprise  that  does  not  need  to  consider  the  logistical 
support? 

I  would  like  to  answer  these  questions  by  pointing  out  the  importance  of  logistics 
with  a  Turkish  quote:  “No  matter  how  much  an  instructor  knows  in  his  area,  he  is  a  good 
instructor  only  as  much  as  his  students  comprehend  the  material.”  It  seems  there  is  no 
relation  between  this  quote  and  logistics,  but  I  will  clarify  the  relation  later  in  this  section. 
If  the  capacity  of  a  student  is  limited,  then  you  should  plan  and  apply  the  education 
program,  and  finally  evaluate  the  student’s  performance  by  placing  an  exam  accordingly. 

Think  about  a  Dessert  Storm  /  Dessert  Shield  scenario.  Instructors  correspond  to 
Staff  Commanders  from  all  countries  attending  this  operation,  and  students  correspond  to 
operational,  logistical,  and  medical  units,  etc.,  that  are  going  to  execute  the  plan.  It  takes 
some  time  to  be  ready  for  the  exam  that  corresponds  to  real  operations.  Students  are 
provided  with  pens,  pencils,  books,  notebooks  as  well  as  additional  class  notes. 

Instructors  ask  some  basic  questions  to  help  them  understand  the  material,  and  assign 
some  homework  and  projects  to  increase  the  readiness  level  for  the  exam.  Assuming 
students  will  pass  the  exams  can  mislead  us,  because  we  are  skipping  the  time  process 


and  the  availability  of  materials  to  get  ready  for  the  exam.  If  the  students  are  responsible 
for  exams,  we  should  provide  them  with  the  materials  they  need  and  time  to  study  for  the 
exam. 

When  we  go  back  to  the  Dessert  Storm  /  Dessert  Shield  Operation  example,  we 
assumed  that  all  units  are  trained  under  operational  conditions.  It  means  that  any  unit  can 
respond  to  changing  situations  very  quickly  and  efficiently  in  terms  of  answering 
questions  as  a  student  during  class.  They  can  get  the  mission  done  with  the  provided 
equipment  in  terms  of  turning  in  projects,  homework  and  passing  exams  as  a  student. 

As  a  result,  all  students  are  good  enough  as  long  as  the  instructors  know  their 
student’s  capacity  and  let  them  to  study  for  exams  with  the  appropriate  tools.  We  cannot 
assume  that  all  students  will  pass  the  exam  without  time  to  study  and  the  appropriate 
tools  to  apply.  That  is  the  reason  the  United  Nations  (UN)  brought  logistical  materials 
first,  months  before  bringing  any  operational  units.  Here  we  can  relate  this  situation  to 
our  example  as  follows.  Instructors  that  correspond  to  UN  staff  personnel  provide  the 
materials  for  students  that  correspond  to  any  kinds  of  operational  units. 

We  find  out  that  logistics  is  crucial  by  just  looking  at  the  preparation  time.  The 
instructor  and  student  interaction  is  forever  while  they  are  in  education  process.  This  is 
true  for  logistics  as  Bowersox  et  al.  states  out  that  logistics  is  unique:  it  never  stops! 
Logistics  is  happening  around  the  globe,  twenty-four  hours  of  every  day,  seven  days  a 
week  during  fifty-two  weeks  a  year.  (Bowersox  and  others,  1996:3) 
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Supply  Chain  Management 


If  logistics  is  a  crucial,  important,  on-going  process,  then  what  is  the  definition  of 

it?  We  can  see  a  lot  of  titles  used  for  logistics  partly,  or  completely  in  the  literature  such 

as  business  logistics,  physical  distribution,  materials  logistics  management,  materials 

management,  physical  supply,  logistics  of  distribution,  marketing  logistics,  inbound 

logistics  and  total  logistics.  In  1991,  Bowersox  et  al.  states  that 

the  Council  of  Logistics  Management  modified  its  1976  definition  of  physical 
distribution  management  by  first  changing  the  term  to  logistics  and  then  changing 
the  definition  as  follows:  “Logistics  is  the  process  of  planning,  implementing  and 
controlling  the  efficient,  effective  flow  and  storage  of  goods,  services  and  related 
information  from  the  point  of  origin  to  the  point  of  consumption  for  the  purpose 
of  conforming  to  customers  requirements.”  (Bowersox  and  others:4) 


After  we  defined  logistics,  we  can  figure  out  the  relationship  between  logistics 

and  supply  chain  management  by  looking  at  the  definition  Chan,  et  al.  presents: 

Integrated  logistics,  or  supply  chain  management  (SCM)  is  the  integration  of  all 
key  business  processes  from  end-user  through  the  original  suppliers  that  deliver 
the  right  products  to  consumers  at  the  lowest  possible  time  and  cost.  The 
elements  of  supply  chain  management  such  as  cost,  quality,  timeliness  and 
flexibility  should  be  accomplished  in  every  phase  of  the  supply  chain 
management.  (Chan  and  others,  Sustainment  Logistic  Proposal:!). 


SCM  allows  us  to  deal  with  different  ways  of  implementing  logistics  such  as 
logistics  planning,  supply,  production,  distribution,  transportation,  inventory  etc.  at  the 
same  time  in  contrast  to  earlier  conceptions  of  logistics.  Bowersox,  et  al.  lays  down  this 
change  in  the  conception  of  logistics  as  a  paradox.  Logistics  has  been  performed  since 
the  beginning  of  civilization:  it  is  hardly  new.  However,  implementing  best  practices  of 
the  integration  of  information,  transportation,  inventory,  warehousing,  materials 
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handling,  and  packaging,  etc.  has  become  one  of  the  most  exciting  and  challenging 
operational  areas  of  business  and  public  sector  management  (Bowersox  and  others  :4). 

We  can  find  another  definition  of  SCM  reflecting  the  production  aspect  of  the 
business  in  the  literature.  Domier  et  al.  calls  SCM  the  management  of  activities  that 
transform  raw  materials  into  intermediate  goods  and  final  products,  and  that  deliver  those 
final  products  to  customers  (Domier  and  others:  1998). 

Statement  of  the  Problem 

We  have  node  0  (main  depot)  and,  nodes  1, 2, 3  (3  bases)  and  2  vehicles.  In  the 
beginning  of  the  period,  each  base  depot  has  an  initial  inventory.  There  is  a  minimum 
and  maximum  inventory  level  for  each  base  node.  Failure  parts  arrive  at  the  bases  with  a 
Uniformly  distributed  mean  X-,  (i=l,  2,  3).  Every  node  have  a  repairing  capabilities  with 
an  Uniformly  distributed  service  rate  pi  (i=0, 1,  2,  3). 

We  will  schedule  the  vehicles  whenever  the  inventory  level  goes  below  its 
minimum  level  to  minimize  costs  such  as  inventory  cost,  and  shipping  cost.  The 
inventory  cost  is  uniformly  distributed  in  every  node.  In  other  words,  the  shortage  and 
surplus  cost  is  the  same  regardless  of  the  number  of  the  level  of  the  demand. 

There  is  a  cost  associated  with  inventory  such  as  holding,  shortage  and  surplus 
cost.  There  is  a  fixed  cost  associated  with  repairing  the  item  in  every  node.  There  is  a 
cost  associated  with  using  vehicles  for  shipping  serviceable  parts  in  the  inventory  from 
any  node,  to  nodej  where  i,  j  =  0,1, 2,3. 

We  assumed  that  we  have  one  type  of  cargo.  We  can  re-supply  from  any  node  to 
another  within  the  system.  But  we  can  only  provide  direct  shipments  between  nodes. 
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We  assumed  that  we  know  the  steady  state  failure  rates  associated  with  each  base, 
and  steady  state  service  rate  associated  with  each  node  (main  depot  and  three  bases).  In 
addition  to  that  we  know  the  shipping  cost  from  node  i  to  node  j  in  terms  of  time  where 
i,  j  =  0,1, 2, 3.  We  know  the  initial  inventory  level,  maximum  and  minimum  inventory 
level  for  each  base  depot  and  main  depot  as  well. 

We  decide  which  vehicle  to  use  and  how  many  parts  to  deliver  from  node  i  to 
node  j  including  both  repaired  items  at  repaired  shops  and  inventoried  items  at  depots  for 
each  node. 

Background 

When  we  examine  the  definition  of  logistics,  we  see  that  it  consists  of  several 
areas.  It  is  in  the  industrial  area  while  producing  goods. 

It  is  in  the  inventory  area  after  we  produce  but  before  we  transport  the  goods  from  the 
production  plant  (depot)  to  the  customers  (demand  points)  stationed  at  various  points. 
Initial  inventory  position,  maximum  inventory  position  which  is  equivalent  to  the 
depot/base  storage  capacity,  and  distribution  of  the  demand  at  each  demand  point  as  well 
as  various  kinds  of  costs  such  as  holding  cost,  surplus  cost,  and  shortage  cost  are  all 
important  for  both  inventory  and  transportation  issues.  Sherbrook  asks  some  important 
questions  such  as  “How  can  we  insure  that  95%  of  our  scheduled  aircraft  flights  will  not 
be  delayed  for  lack-of  spare  parts?”  Additionally  he  asks  some  more  supplementary 
questions  like  what  would  be  the  optimal  level  in  order  to  answer  this  question 
satisfactorily.  More  generally  “What  can  we  do  to  change  our  logistics  support  structure 
to  achieve  a  desired  availability  more  efficiently?”  Or  “  Is  it  economical  to  have  more 
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repair  capability  at  the  operating  sites ?”(Sherbrook,  1992:2).  All  these  types  of  questions 
are  important  when  we  examine  the  problem  in  the  inventory  area. 

It  is  in  the  transportation  area  -or  more  specifically  vehicle  routing  area  if  we 
decide  to  ship  the  goods  with  a  fleet  of  vehicles  to  the  demand  points  provided  that  there 
is  a  need  to  be  visited  by  any  vehicle  due  to  the  lack  of  inventory  level.  Various  time 
windows  associated  with  each  demand  point  can  restrict  delivery  times. 

When  we  add  the  routing  component  into  the  problem,  it  becomes  a  routing- 
inventory  problem.  Transportation  requires  some  consideration  about  the  inventory  level 
of  a  demand  point.  Stenger  explains  the  tradeoff  between  inventory  and  transportation  in 
his  seminar  held  at  the  University  of  Cincinnati  as  follows.  As  we  do  things  to  reduce 
transportation  costs  (bigger  shipment  sizes  and/or  slower,  less  reliable  modes),  we 
increase  the  need  for  inventories  (cycle  stocks  and  safety  stocks)  at  the  receiving  point 
(and  possibly  at  the  shipping  point)  (Stenger,  Seminar  at  UC).  It  is  of  great  of  interest  if 
we  examine  these  two  areas  together.  We  generally  try  to  answer  the  questions  as 
F.Baita  et  al.  explain  in  their  paper. 

When  shipments  have  to  be  made;  i.e.  when  trucks  have  to  be  loaded  and  when 
customers  have  to  be  visited;  how  much  each  truck  has  t  o  be  loaded,  in  terms  of 
quantity  of  each  item  under  consideration,  and  how  such  a  load  has  to  be 
distributed  between  the  requesting  customers;  which  route  has  to  be  followed  by 
each  truck  in  order  to  visit  its  customers  (Baita  and  others,  1998:586). 

Almost  all  of  the  inventory  and  vehicle  routing  problems  try  to  find  answers  for 
the  questions  above.  There  may  be  dynamic  attributes  associated  with  each  demand 
point.  In  other  words,  the  demand  for  parts  can  be  unknown  and  subject  to  change  due  to 
prior  decisions  about  the  three  main  questions  mentioned  above.  F.  Baita  et  al.  classify 
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dynamic  routing-inventory  (DRAI )  problems  in  terms  of  seven  main  groups.  These 
seven  main  groups  are  topology  of  the  problem  such  as  one-to-many,  the  number  of  items 
in  the  problem  such  as  one  item  in  our  problem,  the  attributes  of  demand  in  the  problem 
such  that  it  can  be  known  and  lies  between  two  boundaries  like  our  problem,  the  type  of 
constraints  in  the  problem  associated  with  trucks,  stock,  supply  capacities,  costs 
associated  with  inventory,  repairing,  and  distribution,  and  finally  solution  approach  to  the 
problem  such  as  mixed  programming  models  including  linear,  non-linear  and  integer 
programming  (Baita  and  others,  1998:587).  We  will  cover  each  element  in  greater  details 
in  chapter  two. 

We  will  assume  that  the  demands  for  parts  will  be  unknown  beforehand.  We 
cannot  know  what  kinds  of  malfunction  in  the  aircraft  will  cause  the  pilot  to  return  to 
base  without  completing  his  mission.  Even  though  the  demand  for  parts  is  not 
deterministic  (known),  we  can  take  advantage  of  uniformly  distributed  steady  state  values 
for  failure  and  repair  as  well  as  inventory. 

In  addition  to  that,  repair  comes  into  play  when  demand  points  have  a  limited 
repair  capabilities  as  well  as  preventive  and  corrective  maintenance  over  capital 
equipment  which  directly  affect  the  various  kinds  of  inventory  and  transportation  costs. 

As  McGrath  states,  capital  equipment  maintenance  is  more  important  if  you  have  large 
factories,  energy  generation  power  plants,  refineries,  and  organizations  that  operate  large 
fleets  of  land,  air,  or  seagoing  vehicles  (Me  Grath,  1999:3 1).  We  can  think  of  Air  Forces 
in  terms  of  large  factories  that  keep  huge  number  of  items  and  large  fleets  of  air  and  land 
vehicles  in  inventory.  In  our  case,  flight  units,  customers,  bring  failed  item,  and  request 
for  a  serviceable  item.  As  long  as  we  get  the  failed  item,  we  want  to  repair  it  by  using 
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local  resources  if  it  is  possible.  If  it  is  not  possible,  we  then  take  advantage  of  the  main 
depot  repair  capabilities.  Repair  capabilities  have  two  folds.  One  of  them  is  repairing  the 
coming  failed  items  out  of  aircraft  at  the  base  repair  shop,  and  the  other  one  is  repairing 
capital  equipment  such  as  land  vehicles  used  for  carrying  personnel  and  material  on,  off- 
base,  machines,  tools  and  equipment  used  for  production,  testing  and  calibrating  etc.  We 
are  concerned  with  the  former  one  in  our  model. 

When  we  include  inventory,  repair  and  vehicle  routing  with  time  windows  we 
come  up  with  a  breadth  representation  of  SCM.  Notice  that  repair  is  a  matter  of  military 
or  more  generally  government  sector  rather  than  that  of  private  sector.  It  is  unlikely  to  see 
such  an  example  that  involves  all  three  main  areas.  Although  we  do  not  want  to  impose 
that  we  have  to  accomplish  it  within  one  unique  model,  we  want  to  go  through  a  case  that 
requires  combining  all  of  them  in  SCM  with  time  windows.  If  this  is  the  case,  we  can 
have  a  breadth  understanding  of  the  issue.  When  we  consider  the  application  of  this  issue 
on  Turkish  Air  Force,  we  have  an  insight  what  would  be  advantages  and  disadvantages  if 
we  had  a  combined  model  as  we  developed  in  this  thesis. 

Scope  and  Limitations 

I  am  going  to  scope  the  problem  with  respect  to  the  nature  of  the  problem. 

Keeping  in  mind  that  if  we  choose  the  problem  very  broad  we  cannot  investigate  the 
problem  very  deep.  In  other  words,  there  would  be  a  tradeoff  between  going  breadth  and 
depth.  This  criterion  would  be  determined  by  the  problem  itself. 
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I  will  try  to  come  up  with  a  model  including  repair,  inventory  and  transportation 
areas.  The  assumptions  made  in  the  model  should  match  very  close  to  Turkish  Air  Force 
applications.  We  can  list  these  assumptions  as  follows. 

A.  The  model  will  deliver  whenever  bases  demand  for  items. 

B.  The  model  will  use  all  available  resources  at  base  before  placing  a  request  for 
item. 

C.  The  model  will  use  all  available  vehicles  if  possible. 

D.  If  there  exists  a  delivery  of  an  item  into  node  i,  then  no  delivery  goes  out  of 
node  i. 

E.  The  model  will  deliver  and  pick-up  only  “serviceable”  items.  We  assumed 
that  we  run  the  model  one  time  for  a  steady  state. 

The  purpose  of  this  model  can  be  summarized  in  three  steps.  In  the  first  step,  we 
will  analyze  and  comprehend  the  inventory/allocation  model  (Federgruen  and  Zipkin, 
1984: 1019).  Here  we  take  an  advantage  of  a  numerical  example  in  Chan’s  textbook 
(Chan,  1999-Draft:  9-6).  In  the  second  step,  we  will  apply  the  generalized  Benders 
decomposition  technique  to  the  problem  as  Federgruen  and  Zipkin  did.  Benders  cut  is 
playing  a  major  role  in  the  solution  process.  Since  the  Benders  cut  is  composed  of  dual 
variables,  we  will  perform  a  sensitivity  analysis  for  the  dual  variables  associated  with 
main  depot  capacity  and  vehicle  capacity  constraints.  We  should  come  up  with  the  same 
result  as  Federgruen  and  Zipkin  did.  It  is  summarized  as  follows:  Combining  vehicle 
routing  and  inventory  should  save  a  lot  of  resource  in  terms  of  time  or  money.  In  the 
third  step,  we  extend  the  starting  point  model,  which  corresponds  to  our  thesis  model,  in 
such  a  way  that  we  let  the  model  issue  the  lateral  supply  between  every  nodes  in  the 
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model.  We  track  the  items  individually;  namely,  they  may  come  from  either  inventory  or 
repaired  shop  of  the  node.  We  will  examine  the  results  whether  combining  vehicle 
routing  and  inventory  is  useful  in  this  context  or  not. 

The  stopping  point  would  be  a  model  related  with  all/some  of  the  areas.  After  we 
would  take  advantage  of  Air  Force  Material  Command  Studies  and  Analysis  Office  and 
Turkish  Air  Force  senior  officers  assigned  at  Wright  Patterson  AFB  to  validate  the  model 
we  would  perform  the  sensitivity  analysis  related  to  comparisons  and  relations  in  each 
area  described  above. 

One  of  the  biggest  limitations  for  this  study  is  to  find  a  real  data  for  failure  arrival 
rate  and  repair  service  time. .  The  model  assumes  that  we  have  enough  capabilities  and 
resources  so  that  we  have  a  very  past  data  to  obtain  a  mean  value  for  the  repairing 
process.  We  would  make  up  the  data  to  run  the  model.  In  addition  to  that  I  will  not  go 
through  repairing  process  itself,  as  well. 
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Chapter  2  Literature  Review 


Introduction 

This  chapter  presents  a  general  overview  of  the  literature,  which  is  closely  related 
to  the  application  of  vehicle  routing,  repairing,  and  inventory  problems  in  SCM.  We 
would  explain  the  components  of  SCM  first  and  then  give  a  general  terminology  and 
basic  principles  related  to  each  area. 

The  Relationship  between  Vehicle  Routing  and  Traveling  Salesman  Problem 

As  we  stated  in  chapter  1,  the  definition  of  SCM  requires  delivery  requirements 
from  origin(s)  (depot)  to  destinations  (bases)  via  a  fleet  of  vehicles.  We  recall  a  very 
well  known  type  of  problem  such  as  traveling  salesman  problem  (TSP)  when  we  are 
considering  visiting  bases.  Since  many  scientists  put  their  efforts  into  solving  TSP,  we 
can  easily  find  a  very  clear  and  concise  definition  of  TSP  in  the  literature.  According  to 
Chan,  if  the  origin  is  the  same  as  the  destination,  the  path  is  called  a  tour  with  the 
origin/destination  as  the  home  base  or  the  depot.  Such  a  tour  is  often  referred  to  as  a 
hamiltonian  circuit  or  cycle.  Generally,  we  want  to  find  the  shortest  tour  that  visits  all  the 
nodes  exactly  once  in  a  depot-based  continuous  tour.  This  is  formally  known  as  the  TSP. 
(Chan,  1999Draft:8-16).  Note  that  you  can  start  your  travel  from  an  origin  node,  which 
can  be  any  node  in  the  tour,  and  come  back  to  the  origin  node  provided  we  visit  each 
node  exactly  once  in  the  tour.  It  means  that  we  want  to  solve  a  TSP  in  order  to  come  up 
with  the  shortest  tour  satisfying  a  minimization  objective  function  that  can  be  expressed 
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in  terms  of  money,  time,  or  any  kinds  of  value  we  can  measure.  We  have  more  than  one 
traveling  salesman  in  the  generalized  TSP. 

We  can  take  a  look  at  TSP’s  in  terms  of  a  distance  matrix.  It  would  not  take  the 

same  amount  of  time  if  a  vehicle  passes  through  city  during  rush  hour  rather  than  in  the 

morning  at  two  o’clock.  This  idea  is  valid  for  aircraft  as  well.  Going  from  point  A  to  a 

point  B  with  the  back  wind  would  not  be  the  same  as  traveling  the  same  distance  with  the 

head  wind.  Chan  points  out  this  point  as  follows. 

If  the  distance  matrix  is  symmetrical,  TSP  does  not  depend  on  the 
direction  of  travel  and  every  node  has  two  arcs  incident  to  it.  For 
asymmetrical  distance  matrix,  however,  there  is  one  arc  entering  a  node 
and  another  out  of  a  node  specifically.  The  former  is  referred  as  a 
symmetrical  TSP,  and  the  latter  an  asymmetrical  TSP.  Computationally 
speaking,  symmetrical  TSP’s  are  much  easier  to  solve. 

(Chan,  1999Draft:8-19). 


It  seems  to  us  that  nothing  is  wrong  with  the  formulation  and  solution  of  TSP. 

But  complete  enumeration  takes  a  extremely  long  time  as  Reeves  points  out: 

As  starting  point  is  arbitrary,  there  are  clearly  (N-l)!  possible  solutions  (or 
(N-l)!/2  if  the  distance  between  every  pair  of  cities  is  the  same  regardless 
of  the  direction  of  travel).  Suppose  we  have  a  computer  that  can  list  all 
possible  solutions  of  a  20  city  problem  in  1  hour.  Then  using  the  above 
formula,  it  would  clearly  take  20  hours  to  solve  a  21 -city  problem,  and 
17.5  days  to  solve  a  22-city  problem;  a  25-city  problem  would  take  nearly 
6  centuries.  (Reeves,  1995:7). 


Because  of  the  exponential  growth  in  the  complete  enumeration,  we  should  be 
very  careful  in  selecting  the  number  of  nodes  in  the  problem. 

We  can  define  various  kinds  of  problems  originating  from  TSP.  Chan  gives  a 
definition  of  the  vehicle  routing  problem  (VRP)  related  to  TSP  as  follows.  If  the  delivery 
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requirements  are  placed  upon  the  various  bases  (demand  points)  of  a  TSP,  we  end  up 
with  a  VRP. 

LaPorte,  Louveaux,  and  Mercure  point  out  the  similarity  between  multiple  TSP 
and  VRP  as  follows: 

The  classical  VRP  consist  of  optimally  designing  vehicle  routes  from  one 
or  several  depots  to  a  set  of  customers  in  such  a  way  that: 

(i)  All  vehicles  start  and  end  their  journey  at  the  same  depot 

(ii)  All  customers  are  served  once  by  exactly  one  vehicle,  but  a  vehicle 
route  may  include  several  customers 

(iii)  Some  side  constraints  on  the  routes  are  satisfied 

(iv)  The  sum  of  vehicle  utilization  costs  and  of  routing  costs  is 
minimized  (La  Porte,  1989:71). 


LaPorte  et.al  points  out  the  assumptions  for  the  classical  VRP  as  follows: 


A.  The  demands  of  cargo  at  each  destination  are  fixed. 

B.  The  starting  /ending  locations  are  fixed.(La  Porte,  1989:72). 

We  can  relax  the  third  similarity  in  such  a  way  that  we  do  not  send  a  vehicle 
unless  there  is  a  need  at  the  bases.  In  addition  to  that  we  can  relax  the  first  assumption 
made  in  the  classical  VRP  in  such  a  way  that  demands  can  be  defined  by  a  random 
variable.  These  two  relaxations  make  the  classical  VRP  convenient  for  solving  real 
world  issues. 

The  Relationship  between  Vehicle  Routing  and  Inventory  Routing  Problem 

By  the  solution  of  TSP,  we  obtain  an  order  to  visit  each  demand  points  exactly 
once  so  that  we  minimize  the  cost  associated  with  each  traveling  time  between  depot  and 
demand  points.  If  we  turn  it  into  VRP,  then  our  objective  function  as  Haughton  et  al. 
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points  out  is  to  find  a  set  of  delivery  routes  that  simultaneously  satisfies  demand  at  each 

retail  outlet  (demand  points)  and  minimizes  total  transportation  cost 

(Haughton,  1999:25).  Note  that  we  add  the  delivery  requirements  to  various  demand 

points  in  VRP  without  considering  inventory  issues  at  the  demand  points.  Dror  and  Ball 

give  a  definition  for  inventory  routing  problem  (IRP)  as  follows. 

The  IRP  involves  a  set  of  customers,  where  each  customer  has  a  different 
demand  on  each  day.  For  example,  each  customer  uses  a  commodity  such 
as  heating  oil  or  methane  at  an  estimated  consumption  rate.  Each 
customer  possesses  a  known  capacity  (for  example,  the  size  of  his  tank  to 
hold  home  heating  oil).  The  objective  is  to  minimize  the  annual  delivery 
costs  while  attempting  to  insure  that  no  customer  runs  out  of  the 
commodity  at  any  time.  The  impetus  behind  this  distribution  system  is  the 
importance  of  maintaining  a  sufficient  supply  of  inventory  at  a  customer’s 
location  (Dror,  1987:891). 


When  we  examine  the  definition  of  IRP,  we  realized  that  we  should  consider  the 
annual  delivery  as  an  objective  function  to  minimize.  We  can  use  a  long-term  delivery  or 
steady  state  delivery  instead  of  annual  delivery.  We  will  mention  a  lit  bit  more  on  annual 
delivery  later. 

The  demand  of  each  customer  for  the  IRP  is  small  compared  with  the  capacity  of 
the  vehicle.  IRPs  have  inventory-related  costs  and  incentive  mechanisms  for  early 
delivery  to  a  customer. 

We  see  some  other  definitions  of  IRP  in  the  literature  under  the  name  of 
allocation/routing  or  assignment  routing  problems  (Dror,  1987:892). 

We  extend  classical  VRP  by  adding  periodic  horizon  time  in  order  to  get  a 
similarity  to  the  IRPs.  Christofides  et  al.  explains  the  Periodic  Vehicle  Routing  Problem 
(PVRP)  as  follows. 
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In  PVRP  the  problem  is  to  design  a  set  of  routes  for  each  day  of  a  given 
(p-day)  period.  Each  customer  may  require  a  number  of  visits  by  a  vehicle 
during  this  period.  If  a  customer  requires  k(say)  visits  during  the  period, 
then  these  visits  may  only  occur  in  one  of  a  given  number  of  allowable  k- 
day  combinations  (Christodifides,  1984:237). 


PVRP  is  similar  to  IRP  in  the  sense  of  long-term  period  and  satisfying  customer 
demand.  However,  it  is  different  in  the  sense  of  specifying  the  number  of  visits  during 
the  period.  The  decision-maker  has  a  power  to  determine  the  number  of  visits  during  the 
period  in  PVRP  whereas  IRPs  let  the  model  decide  the  number  of  visits  and  the  size  of 
delivery. 

According  to  Russel  et  al.  we  can  consider  the  routing  aspect  of  IRP  as  follows. 

Routing  design  problem  in  which  the  objective  is  to  assign  customer 
demand  points  to  days  of  the  week  in  such  a  way  that  the  resulting  node 
routing  problems  yield  a  near-optimal  solutions.  It  is  assumed  that 
demand  points  offer  some  flexibility  in  their  assignment  to  days  of  the 
week;  each  point  may  require  service  anywhere  from  one  to  seven  days 
per  week.  This  routing  design  problem  assumes  that  the  resulting  node 
problem  on  each  day  of  the  week  is  a  single  depot  vehicle  dispatch 
problem  whose  objective  is  to  minimize  the  distance  or  time  required  to 
service  customer  demand  points  and  to  minimize  the  number  of  vehicles 
required  (Russel,  1979:page  1). 


Routing  design  problem  (RDP)  is  similar  to  IRPs  in  the  sense  of  assigning 
vehicles  to  demand  points;  however,  since  the  service  day  of  any  node  in  the  model  is 
predetermined  it  is  different  from  IRPs  in  that  respect.  The  primary  objective  function  of 
the  model  is  to  minimize  total  distance  (time)  traveled  per  week  subject  to  some 
constraints  very  similar  to  classic  VRPs.  Demand  points  may  require  being  visited  one, 
two  or  up  to  six  days  of  the  week  assuming  six  business-day  in  one  week.  The  model 
decides  the  assignment  of  vehicles  after  demand  is  determined;  namely,  decision-maker 
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puts  a  control  mechanism  on  the  demand.  This  type  of  formulation  is  valid  for  either 
delivery  or  pick-up  type  problems  such  as  refuse  collection,  scheduled  retail  and 
wholesale  delivery  problems  as  Russel  et  al.  explain. 

According  to  Federgruen  et  al.,  we  can  consider  a  combined  vehicle  routing  and 
inventory/allocation  problem  very  similar  to  IRPs.  The  authors  point  out  this  problem  as 
follows. 

There  are  many  situations  where  the  vehicle  schedules  and  the  delivery 
sizes  are  (or  should  be)  determined  simultaneously.  Such  is  often  the  case, 
when  at  each  location  the  demand  for  the  resource  is  random.  Here, 
deliveries  serve  to  replenish  the  inventories  to  levels  that  appropriately 
balance  inventory  carrying  and  shortage  costs,  but  thereby  incur 
transportation  costs  as  well  (Federgruen,  1984:1019). 

A  combined  vehicle  routing  and  inventory/allocation  problem  is  very  similar  to 
IRPs  in  terms  of  delivering  items  to  demand  points  via  fleet  of  vehicles  over  a  period.  As 
Federgruen  states  out  this  type  of  problem  is  valid  for  internal  distribution  problems  such 
as  delivering  fuel  oil  to  automotive  service  stations.  All  decisions  are  made  centrally  in 
an  internal  distribution  problem.  Besides  that,  the  model  can  be  applied  to  some  external 
distribution  applications  such  as  gas  producers  themselves  install  tanks  at  customer 
locations  and  determine  the  replenishment  frequency  and  delivery  sizes  as  Federgruen 
states  out  as  well  (Federgruen,  1984:1019).  However,  decision-maker  does  not  have  any 
control  mechanism  on  demand  as  Russel  et  al.  have  in  this  type  of  model.  In  addition  to 
that,  the  model  does  not  need  requiring  k- visits  to  node  i  per  period  as  Christofides  et  al. 
do.  The  period  can  be  day,  week,  month,  or  even  a  year  in  Federgruen’s  model.  But 
Dror  et  al.  state  out  some  reasons  not  to  consider  one  period  as  a  year  as  follows. 
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1.  The  validity  of  a  formulation  with  an  annual  time  base  is  questionable 
because  some  parameter  values  are  uncertain  over  that  long  a  period, 
others  change. 

2.  The  number  of  constraints  and  variables  for  such  a  formulation  is 
prohibitively  high  (Dror,  1987:891). 


As  a  starting  point  of  understanding  the  concepts  of  delivery  models  with 

inventory  concern,  we  decided  to  solve  a  toy  problem  that  includes  a  main  depot,  two 

capacitated  vehicles,  and  three  bases  embellished  with  newsboy  inventory  (Chan,  1999 

Draft:  9-6).  The  concepts  and  assumptions  of  the  toy  problem  is  the  same  as  what 

Federgruen  et  al.  cover  in  their  studies  as  follows. 

The  initial  inventory  (perhaps  supply  remaining  from  the  previous  day)  for 
each  location  is  reported  to  the  depot.  This  information  is  used  to 
determine  for  the  following  day  the  allocation  of  the  available  product 
among  the  locations.  The  assignment  of  locations  to  the  vehicles  and  the 
routes  are  set  at  the  same  time.  After  deliveries  are  made  (  say  at  the  end 
of  day)  the  demands  occur,  and  inventory-carrying  and  shortage  costs  are 
incurred  at  each  location  proportional  to  the  end-of-the-day  inventory 
levels  (Federgruen,  1984:1020). 


We  try  to  minimize  the  transportation  and  inventory  cost  in  this  toy  problem. 

The  results  show  us  that  the  total  cost  is  driven  by  inventory  cost.  A  small  portion 
of  the  total  cost  is  due  to  transportation  costs.  Therefore,  we  place  a  close  examination 
into  the  inventory  aspect  of  the  issue. 

There  are  a  lot  of  delivery  models  integrated  with  inventory  concern  in  the 
literature  as  well  as  Baita  et  al.  and  Herer  et  al.  cover  in  their  studies,  respectively.  We 
will  touch  the  routing-and-inventory  models  in  great  detail  in  the  following  section. 
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Dynamic  Routing-and-Inventory  problems 


Baita  et.al.  present  three  basic  aspects  such  as  routing,  inventory,  and  dynamicity 
that  we  can  embed  in  VRP  and  IRP  to  form  dynamic  routing-and-inventory  problems 
(DRAI).  It  includes  routing  because  VRP  requires  organizing  the  physical  movement  of 
goods  between  different  geographic  sites,  such  as  depots,  warehouses,  production  and 
retail  points,  etc.  It  includes  inventory  because  elements  such  as  quantities  and  values  of 
the  goods  being  moved  are  relevant  to  the  definition  and  assessment  of  organizational 
and  operational  strategies.  It  has  a  dynamicity  property  because  routing  and  inventory 
are  imbedded  in  a  dynamic  environment.  In  other  words,  repeated  decisions  have  to  be 
taken  at  different  times  within  some  time  horizon,  and  earlier  decisions  influence  later 
decisions  (Baita,  1998:585). 
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Table  1:  The  elements  for  a  classification  of  DRAI. 


Classification  element 

Attribute 

Alternatives 

Topology 

Endpoints 

One-to-one,  Many- 
to-many 

Items 

Number 

One-to-one,  Many 

Demand 

Knowledge 

Known,  uncertain 

Time  behavior 

Constant,  variable 

Distribution 

Uniform,  not 
uniform 

Decision  type 

Domain 

Type,  frequency 

Constraints/objective 

Truck  capacity 

Stock  capacity 
Supply  capacity 

Equal,  different 

Vehicle  number 

Given 

Costs 

Inventory 

Holding 

Transportation 

Fixed 

Solution  approach 

Decompositon 

Time,  Cluster  route 

Clustering 

Time,  frequency 

Model  /  Algorithm 

Exact,  approximate 

MP  models 

LP,  IP,  NLP 

After  we  lay  out  the  general  attributes  of  DRAI,  we  would  like  to  explain  the 
classification  elements  of  DRAI  problems  that  we  presented  in  chapter  one. 

The  first  classification  element  is  topology,  which  can  be  one-to-one,  one-to- 
many,  or  many-to-many.  For  example;  one-to-many  denotes  that  we  have  one  origin 
(depot  or  distribution  center)  and  many  destinations  (demand  points  or  bases). 

The  second  one  is  number  of  items  we  consider  can  be  a  single  item  or  several.  F. 
Baita  et  al.  points  out  the  importance  of  dealing  with  several  items  independently.  This 
idea  may  be  an  important  even  though  we  have  one  item  in  the  problem  setting  as  Bauer 
et  al.  point  out.  Here  we  assumed  that  one  item  corresponding  to  one  type  of  item  and 
multi-item  corresponding  to  multi-type  items.  While  we  are  studying  the  problem  that 
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Bauer  et  al.  presented  like  a  part  of  our  thesis  problem,  we  see  that  we  should  track  the 
items  independently.  For  instance,  we  may  want  to  know  what  percentage  of  items 
delivered  from  depot  and  demand  point  1  to  demand  point  2  if  we  assumed  that  the 
vehicle  1  follows  the  route  in  the  sequence  of  depot-demand  point  1 -demand  point  2. 
Vehicle  1  pick-ups  some  items  out  of  main  depot  to  deliver  node  1,  pick-ups  some  items 
out  of  node  1  to  deliver  to  node  2.  Besides  that  there  is  another  possibility  that  vehicle  1 
pick-ups  some  items  out  of  main  depot  to  deliver  node  2  or  pick-up  out  of  node  1  to 
deliver  main  depot.  All  these  possibilities  are  valid  for  delivery/pick-up  problems  as 
Bauer  et  al.  solved  in  their  working  paper.  When  we  consider  the  following  situation  in 
terms  of  aircraft  parts,  it  is  not  acceptable  to  deliver  to  node  1  and  pick-up  out  of  node  1 
simultaneously  on  the  same  tour.  For  example,  if  we  deliver  out  of  main  depot  to  node  1, 
then  we  cannot  deliver  out  of  node  1  to  node  2  or  main  depot.  It  is  clear  that  the  desire  to 
deal  with  items  independently  will  cause  some  difficulties  as  the  problem  size  grows. 

The  remaining  properties  of  DRAI  problems  such  as  demand,  decision  type, 
constraints/objective  function,  and  costs  excluding  the  solution  approach  are  the  basic 
among  inventory  models  as  well  as  Tersine  state  out  (Tersine,  1994: 12).  We  will  cover 
these  along  with  DRAI  problems. 

The  third  one  is  the  characteristics  of  demand  that  consist  of  three  elements  like 
knowledge,  time  behavior  and  distribution.  Since  demand  is  a  major  element  in  DRAI 
problems  and  inventory  models,  we  will  cover  this  point  in  depth.  When  we  say 
knowledge  we  mean  that  the  demand  can  be  either  known  or  uncertain.  It  seems  clear  if 
we  have  a  known  demand.  We  can  view  the  uncertain  demand  either  as  a  random 
variable  associated  with  an  upper  and  lower  bound  or  unknown  completely.  Demand  can 
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be  a  constant  value  or  random  value  across  the  time  horizon  when  we  talk  about  time 
behavior  of  demand.  Baita  et  al  points  out  the  distribution  property  of  demand  such  that 
it  can  be  equally  (uniformly)  distributed  or  not  when  several  customers  or  several  items 
are  involved.  We  note  that  we  may  not  have  a  uniformly  distributed  demand  even  though 
we  consider  one  type  of  items.  For  instance,  the  distribution  function  related  to  demand 
points  could  be  Poisson,  Exponential  etc. 

We  can  look  at  the  demand  in  such  a  way  that  Tesine  does.  Demands  can  be 
examined  according  to  their  size,  rate,  and  pattern  (Tersine,  1994: 12).  Demand  size  is 
basically  what  we  need/order  in  terms  of  quantity.  It  is  important  because  it  specifies  the 
nature  of  the  problem  like  deterministic  model  or  stochastic  model.  Then  the  problem 
boils  down  to  a  point  whether  we  know  the  demand  size  or  not.  Taha  specifies  the 
distinction  as  follows: 

In  the  deterministic  case,  it  is  assumed  that  the  quantities  needed  over 
subsequent  periods  of  time  are  known  with  certainty.  This  may  be 
expressed  over  equal  periods  of  time  in  terms  of  known  constant  demands 
or  in  terms  of  known  variable  demands.  These  two  cases  are  referred  to  as 
static  and  dynamic  demands,  respectively.  Probabilistic  demand  occurs 
when  the  requirements  over  a  certain  period  of  time  are  not  known  with 
certainty  but  their  pattern  can  be  described  by  a  known  probability 
distribution.  In  this  case,  the  probability  distribution  is  said  to  be  either 
stationary  or  nonstationary  over  time  (Taha,  1976:390). 

We  know  that  probability  distributions  can  be  discrete  or  continuous.  As 
Wackerly  et  al.  state  that  the  nature  of  the  distribution  function  associated  with  a  random 
variable  determines  the  variable  is  continuous  or  discrete  (Wackerly,  1996:136).  We  will 
go  from  backward  to  understand  the  distribution;  namely,  first  define  what  are  continuous 
or  discrete  random  variables  and  then  the  distributions  associated  with  them.  A  random 
variable  that  can  take  on  at  most  a  countable  number  of  possible  values  is  said  to  be 
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discrete  (Ross,  1997:25).  Wackerly  et  al.  defines  the  probability  distribution  associated 
with  discrete  random  variable  as  follows:  The  probability  that  Y  takes  on  the  value  of  y, 
P(Y=y),  is  defined  as  the  sum  of  the  probabilities  of  all  sample  points  in  S  that  are 
assigned  the  value  y  (Wakerly  et  al.,  1996:77).  It  means  we  can  assign  nonzero 
probabilities  to  only  a  countable  number  of  distinct  y  values.  On  the  other  hand,  the 
random  variable  whose  set  of  possible  values  is  uncountable  is  said  to  be  continuous 
(Ross,  1997:31).  Wackerly  et  al.  defines  the  probability  distribution  associated  with 
continuous  random  variable  as  follows:  Let  Y  denote  a  random  variable  with  distribution 
function  F(y).  Y  is  said  to  be  continuous  if  the  distribution  function  F(y)  is  continuous, 

for  -oo <  Y  <oo  (Wackerly,  1997:139).  As  an  example,  demand  for  an  aircraft  part 

cannot  be  a  fractional  value  that  corresponds  to  discrete  random  variable  such  as  Poisson 
and  probability  distribution  function  of  Poisson  whereas  demand  for  repairing  the  failure 
item  cannot  be  restricted  to  an  integer  number  that  corresponds  to  continuous  random 
variable  such  as  Exponential,  Weilbull  and  continuous  probability  distribution  function  of 
Exponential  etc.  As  Tersine  states  out  that  we  generally  use  standard  probability 
distributions  such  as  Normal,  Exponential,  and  Poisson  for  demand  (Tersine,  1994:12). 

As  a  result,  if  we  know  the  demand  size  with  certainty,  the  inventory  model  is 
deterministic;  otherwise,  it  is  probabilistic. 

After  we  discuss  the  demand  size,  we  can  explain  the  demand  rate  based  upon 
demand  size.  Demand  rate  is  the  demand  size  per  unit  of  time.  Demand  pattern  is 
corresponding  to  the  pattern  that  how  we  are  pulling  the  items  out  of  inventory  such  as  at 
the  beginning  of  the  period,  at  the  end  of  the  period,  uniformly  throughout  the  period. 
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The  fourth  one  is  the  decision  type  involves  two  elements  like  frequency  domain 
and  time  domain.  Baita  et  al.  point  out  that  frequency  domain  decision  type  is 
appropriate  with  synchronous,  periodic  operations.  In  other  words,  decision  variables  are 
replenished  between  frequencies.  On  the  other  hand,  time  domain  decision  type  is 
appropriate  with  asynchronous  operating  models,  possibly  using  feedback  (Baita  and 
others,  1998:586-587).  In  other  words,  we  can  explain  it  as  follows.  The  schedule  of 
shipment,  the  quantities  and  routes  are  decided  at  fixed  time  intervals.  Once  we  fixed 
them  at  the  beginning  of  the  period,  we  keep  them  the  same  until  the  end  of  that  period. 
Herer  et  al.  points  out  this  difference  related  to  his  literature  review  in  the  paper  as 
follows.  The  VRP  is  solved  for  a  single  period  with  known  demands  that  correspond  to 
time  domain  decision  type,  whereas  inventory  routing  problem  (IRP)  is  solved  over 
several  periods  with  stochastic  demands  that  correspond  to  frequency  domain  decision 
type  (Herer,  1997:71). 

We  can  look  at  the  difference  as  Tersine  does  as  well.  He  uses  the  term 
replenishments  for  decision  type.  Replenishments  can  be  examined  according  to  size, 
pattern,  and  lead-time  (Tersine,  1994:13).  Replenishment  size  is  basically  the  quantity  or 
size  of  the  order  to  be  received  into  inventory,  which  can  be  either  constant  or  variable 
depending  on  the  type  of  inventory  model.  If  the  organization  uses  deterministic  type 
model  then  the  replenishment  size  is  known,  on  the  other  hand  if  the  organization  uses 
stochastic  type  model,  then  the  replenishment  size  becomes  variable.  The  replenishment 
pattern  is  all  about  the  way  we  added  the  replenishment  sizes  into  inventory.  It  can  be 
instantaneous,  uniform,  or  batch.  We  should  define  first  the  term  lead  time  in  order  to 
explain  replenishment  lead-time.  The  definition  of  lead-time  is  the  time  between  placing 
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an  order  and  its  receipt  in  the  inventory  system  as  Anderson  et  al.  define  (Anderson  and 
others,  1994:503).  Replenishment  lead-time  is  the  lead-time  of  replenishments.  It  can  be 
constant  or  variable  the  same  as  demand.  We  should  use  probability  distribution 
functions  to  define  the  variable  replenishment  lead-time. 

The  fifth  one  is  the  decision  variables  involved  in  constraints  and  objective 
function  such  as  vehicle  (truck-airplane)  capacity,  stock  capacity  associated  with  demand 
points,  and  supply  capacity  related  to  origin  points.  We  may  have  either  homogeneous 
set  of  vehicles  or  not.  In  addition  to  that  the  number  of  vehicles  we  have  can  be  a 
decision  variable  or  it  is  known  beforehand.  We  usually  assume  that  we  know  the 
capacity  constraints  related  to  supply  points,  demand  points  and  vehicles.  Herer  et  al. 
classifies  the  decision  variables  in  a  three  ways  as  follows.  In  the  strategic  level,  the  main 
decisions  concern  the  location  of  facilities  like  plants,  warehouses,  and  demand  points.  In 
the  tactical  level,  the  main  decisions  concern  the  fleet  size  and  its  attributes  such  as 
vehicle  capacity.  In  the  operational  level,  the  main  decisions  concern  the  routing  and 
scheduling  of  the  vehicles  in  order  to  service  the  customers  with  corresponding  quantities 
to  be  delivered  (Herer,  1997:69).  Since  we  solve  our  problem  in  the  operational  level,  we 
try  to  decide  which  vehicle  to  use  and  how  many  parts  to  deliver  from  node  i  to  node  j 
including  both  repaired  items  at  repaired  shops  and  inventoried  items  at  depots  for  each 
node.  In  other  words,  the  number  of  vehicles,  the  locations  of  depots  are  given  in  our 
problem.  As  Tersine  states  out  that  the  organizations  may  have  some  constraints  such  as 
management  policies  and  administrative  decisions  as  well  as  space  constraints  related  to 
vehicle,  depot,  and  personnel,  capital  constraints.  Management  policies  such  as  never 
being  stock  out  on  certain  items,  and  administrative  decisions  such  as  reciprocal 
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purchasing  agreements  can  be  another  example  to  the  type  of  constraints  we  have  in  the 
inventory  models  (Tersine,  1994:13). 

The  objective  function  of  the  inventory  models  is  what  to  purchase  or 
manufacture,  when  to  take  action,  and  in  what  quantity  as  Bowersox  et  al.  point  out 
(Bowersox  et  al.  1996:250).  By  assuming  that  we  know  what  to  keep  in  the  inventory; 
namely,  skipping  what  to  purchase  or  manufacture  stage,  we  can  revisit  the  objective  of 
inventory  models  as  follows.  The  objective  function  of  inventory  models  is  to  have  the 
appropriate  amounts  of  materials  in  the  right  place,  at  the  right  time,  and  at  low  cost  as 
Tersine  states  (Tersine,  1994:13).  Notice  that  the  minimization  of  cost  is  not  the  sole 
objective  in  the  inventory  models.  The  objective  is  to  minimize  the  inventory  costs  in  the 
deterministic  models  whereas  the  objective  is  to  minimize  the  expected  cost  of  inventory 
in  the  stochastic  models  (Tersine,  1994:210). 

The  sixth  one  is  the  cost  related  to  transportation  and  inventory.  We  usually 
minimize  the  cost  in  the  objective  function  of  DRAI  problems  provided  we  satisfy  the 
demand.  Inventory  cost  can  occur  in  holding  the  items  in  the  depot  in  the  form  of 
holding  cost.  It  occurs  in  stocking  the  items  less  than  or  greater  than  the  demand  in  the 
form  of  shortage  cost,  and  surplus  cost,  respectively  as  well.  We  can  add  the  fixed 
ordering  cost  to  the  objective  function  in  order  to  be  minimized.  Setup  cost  is  another 
issue  to  consider  in  considering  the  repairing  an  item  at  the  facility. 

Before  going  any  further,  I  would  like  to  define  the  various  types  of  costs  within 
the  inventory  models.  Anderson  et  al.  states  the  definition  of  holding  cost  as  follows.  It 
is  the  cost  associated  with  maintaining  an  inventory  investment;  holding  cost  includes  the 
cost  of  the  capital  investment  in  the  inventory,  insurance,  taxes,  warehouse  overhead,  and 
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so  on.  It  can  be  interpreted  as  the  surplus  cost.  If  the  demand  is  less  than  what  we  have 

then  we  will  come  up  with  the  surplus  cost.  Order  cost  includes  such  items  as  making 

requisitions,  analyzing  vendors,  writing  purchase  orders,  receiving  materials,  inspecting 

materials,  following  up  orders,  and  doing  the  processing  necessary  to  complete  the 

transaction  (Tersine,  1994: 14).  Setup  cost  concludes  the  fixed  charge  associated  with  the 

placement  of  an  order  or  with  the  initial  preparation  of  a  production  or  repair  facility 

(Taha,  1976:389).  We  can  think  of  setup  cost  and  order  cost  as  an  independent  of 

demand  size  or  replenishment  size.  As  soon  as  we  place  an  order,  those  costs  will  occur. 

Shortage  (stockout)  cost  is  the  penalty  cost  running  out  of  stock  when  there  is  a  demand 

for  it.  There  exist  two  cases  in  situation  of  shortage  as  Taha  states  out  as  follows: 

They  generally  include  costs  due  to  loss  in  customers’  goodwill  and 
potential  loss  in  income.  In  the  case  where  unfilled  demand  can  be 
satisfied  at  a  later  date  (backlog  case),  these  costs  usually  vary  directly 
with  both  the  shortage  quantity  and  the  delay  time.  If  the  unfilled  demand 
is  lost  (no  backlog  case),  shortage  costs  become  proportional  to  shortage 
quantity  only  (Taha,  1976:390). 

We  can  even  extend  the  shortages  with  respect  to  origins.  It  can  be  internal  shortage 
where  an  order  of  the  department  within  the  organization  is  not  filled.  By  the  same 
token,  it  can  be  external  shortage  where  customer  order  is  not  filled.  We  should  recall  the 
well-known  problem  in  the  literature  named  newsboy  problem  when  we  consider 
shortage  and  surplus  cost  all  together  across  various  levels  of  demands. 

The  transportation  costs  can  show  a  wider  variety  as  Baita  et  al.  points  out.  They 
can  have  a  fixed  component  per  trip  in  particular,  possibly  depending  on  the  load  level  of 
the  trucks,  plus  variable  components  such  as  distance  the  vehicles  traveled,  and  the 
number  of  customers  visited  (Baita,  1998:587). 
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The  seventh  one,  the  last  one,  is  the  solution  technique  to  the  DRAI  problems. 

We  can  solve  these  problems  as  treating  them  as  two  sub-problems.  Herer  et  al.  points 
out  these  sub-problems  as  follows.  Finding  the  optimal  assignment  of  customers  to 
trucks  corresponding  to  clustering  process,  and  finding  the  best  route  corresponding  to 
route  process  (Herer,  1997:70).  Baita  et  al.  presents  two  decomposition  of  time  horizon 
either  of  the  type  cluster-first,  route-second,  or  vice  versa.  The  second  concern  related  to 
solution  technique  is  how  customers  are  aggregated  into  routes  based  on  their 
geographical  position,  or  on  their  estimated  visit  frequency,  or  on  the  estimated  time  to 
run  off  (Baita,  1998:587).  What  kinds  of  algorithms  should  we  use  for  the  solution?  We 
have  two  choices;  it  can  be  either  exact  or  approximate  algorithm.  If  we  decided  to  use 
any  mathematical  programming  (MP)  for  the  formulation,  we  may  come  up  with  the 
situation  such  that  we  may  choose  our  solution  technique  accordingly  based  on 
characteristic  of  the  problem  such  that  it  can  be  linear,  integer,  non-linear,  or  even  mixed 
integer  (Baita,  1998:587).  It  would  be  good  to  present  the  basics  of  these  problems  and 
the  relationships  between  them.  Winston  explains  the  linear  function  and  linear  problems 
as  follows: 

A  function /(xi,  X2, ...,  xn)  of  xi,  X2, ...,  xn  is  a  linear  function  if  and  only  if 
for  some  set  of  constants  c\,  ci,  .. .,  cw 

/(x  1,  X2,  ...,  Xn)  =  CiXi  +  CiX]  +...+CnXi 

For  any  linear  function /(xi,  X2, ...,  xn)  and  any  number  b,  the  inequalities 
/(x i,  X2, ...,  x„)  <  b  and/(xi,  X2, ...,  x„)  >  b  are  linear  inequalities. 
(Winston,  1994:53). 

For  example,  f(x\,  x2)=  2xi  +  X2  is  a  linear  function  but/(xi,  X2)=xi  X2  is  not  linear 
function,  it  is  non-linear  function.  2xi  +  X2  <  5  is  an  example  for  linear  inequalities. 
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Xi.X2  <  2  is  an  example  for  nonlinear  inequalities.  Linear  problems  consist  of  linear 
objective  function,  linear  equalities  and  linear  inequalities  and  sign  restricted  sings  that 
can  be  whether  nonnegative  any  value  such  as  integer  and  continuous,  or  unrestricted  in 
sing.  The  integer  programming  problems  is  the  same  as  linear  programming  problems 
except  the  restriction  on  the  values  of  variables.  An  Integer  Programming  (IP)  in  which 
all  variables  are  required  to  be  integers  is  called  pure  integer  programming  problems  and 
an  IP  in  which  only  some  of  the  variables  are  required  to  be  integers  is  called  a  mixed 
integer  programming  problems  (MIP)  (Winston,  1994:  464).  A  nonliner  programming 
problem  in  which  the  objective  function  or  some  of  the  constraints  may  not  be  linear 
(Winston,  1994:639). 

We  have  not  talked  about  dynamicity  attributes  of  routing-inventory  problems. 
Dynamicity  of  demand  is  much  more  important  than  that  of  vehicles  or  customers  in 
DRAI  problems  at  the  operational  level. 

When  we  consider  the  demand,  which  is  a  major  concern  of  our  interest,  we 
would  analyze  it  as  Baita  et  al.  state  out: 

In  order  to  apply  any  stochastic  model  to  a  concrete  real  situation,  a  preliminary 
analysis  of  the  demand  is  generally  required,  which  in  turn  must  be  based  on  at  least  four 
elements: 

1 .  A  suitable  large  database  of  historical  records. 

2.  Appropriate  assumptions  about  the  future  behavior  of  the  demand. 

3.  Appropriate  assumptions  about  the  form  of  stochastic  processes  that 
could  suitably  represent  the  demand. 

4.  The  estimation,  using  statistical  tools  of  the  parameters  of  those 
stochastic  processes  (Baita,  1998:596). 


As  Baita  et  al.  state  out,  it  is  not  easy  to  carry  out  these  steps  from  a  practical 
point  of  view  due  to  the  stochastic  assumptions  that  could  adversely  affect  the  possibility 
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of  obtaining  a  model  which  completely  adheres  to  the  actual  elements  of  the  real  situation 

under  study.  Because  of  the  above  reasons,  Baita  et  al.  states  out: 

Unknown  demand,  which  only  requires  to  know  the  bounds  within  which 
it  could  lie,  seemed  a  sensible  way  to  get  more  manageable  models  from  a 
practical  point  of  view,  at  least  in  some  situations,  as  far  as  the  process  of 
assessing  the  demand  is  concerned  (Baita,  1998:596). 

An  Overview  to  Inventory:  Definitions 

I  would  like  to  present  a  definition  of  inventory  and  to  whom  inventory  is 
important.  We  can  use  the  inventory  in  different  meanings  as  Tersine  states  out  as 
follows. 


1.  the  stock  on  hand  of  materials  at  a  given  time  (a  tangible  asset 
which  can  be  seen,  measured,  and  counted); 

2.  an  itemized  list  of  all  physical  assets; 

3.  (as  a  verb)  to  determine  the  quantity  of  items  on  hand; 

4.  (for  financial  and  accounting  records)  the  value  of  the  stock  of 
goods  owned  by  an  organization  at  a  particular  time 
(Tersine,  1994:3). 


We  would  use  the  term  inventory  refers  to  the  first  definition.  Tersine  present  a 
generalized  definition  of  first  one  as  follows.  Inventory  as  material  held  in  an  idle  or 
incomplete  state  awaiting  future  sale,  use,  or  transformation  (Tersine,  1994:3). 

After  we  lay  down  the  several  definitions  of  inventory,  we  are  ready  to  explain 
the  organizations  consider  the  inventory  function  as  a  part  of  their  daily  life.  As  Tersine 
presents  that 


The  problems  of  inventory  do  not  confine  themselves  to  profit-making 
institutions  but  likewise  are  encountered  by  social  and  nonprofit 
institutions.  Inventories  are  common  to  farms,  manufactures,  wholesalers, 
retailers,  hospitals,  churches,  prisons,  zoos,  universities,  and  national, 
state,  and  local  governments.  Indeed,  inventories  are  also  relevant  to  the 
family  unit  in  relation  to  food,  clothing,  medicines,  toiletries,  and  so 
forth(Tesine,  1994:3). 
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Since  it  is  widely  concerned,  One  wonders  that  why  all  of  the  organizations  above 
are  keeping  the  same  type  of  inventory? 

Types  of  Inventory 

We  have  four  types  inventory  like  supplies,  raw-materials,  in-process  goods,  and 
finished  goods  as  Tersine  states  out.  Supplies  are  inventory  items  consumed  in  the 
normal  functioning  of  an  organization.  Supplies  are  not  a  part  of  the  final  product.  We 
can  give  some  examples  to  supplies  such  as  pencil,  paper,  disks,  facility  maintenance 
items.  Raw-materials  are  items  purchased  from  suppliers  to  be  used  as  inputs  into  the 
production  process.  Lumber,  glue,  screws  are  good  examples  of  raw-materials  to 
furniture  manifacturer.  In-process  goods  are  partially  completed  final  products  that  are 
still  in  the  production  process.  Finished  goods  are  the  final  products,  which  are  available 
for  sale,  distribution  or  storage  (Tersine,  1994:4).  As  an  example,  manufacturers  hold  all 
types  inventory  while  wholesalers  and  retailers  hold  finished  products  in  their  inventory. 

Since  inventory  is  important  and  organizations  can  handle  different  types 
inventory  at  the  same  time,  one  may  wonder  the  functions  of  it  in  order  to  take  control  of 
it.  We  are  concerning  the  finished  items  in  our  thesis  model. 

Functions  of  Inventory 

When  we  examine  the  first  definition  of  inventory,  we  see  that  it  is  all  about  the 
number  of  stock  on  hand  at  a  given  time.  We  may  ask  this  simple  question  without 
thinking.  Why  do  we  have  to  keep  some  number  of  stocks  on  hand?  Or  by  the  same 
token,  what  are  the  reasons  for  existence  of  inventory?  As  Tersine  states  out  that 
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inventory  exists  for  the  following  reasons  and  he  explains  the  functions  associated  with 
inventory  (Tersine,  1994:4).  Supply  and  demand  frequently  differ  in  the  rates  at  which 
they  respectively  provide  and  require  stock  due  to  four  functional  factors  of  inventory 
such  as  time,  discontinuity,  uncertainty,  and  economy  factor. 

Time  factor  of  inventory  involves  the  long  process  of  production  and  distribution 
required  such  as  time  to  develop  production  schedule,  ship  raw  materials  from  suppliers 
before  goods  reach  the  final  customer. 

The  discontinuity  factor  of  inventory  makes  us  deal  with  dependent  operations 
such  retailing,  distributions,  warehousing,  manufacturing  and  purchasing  in  an 
independent  and  economical  manner. 

The  uncertainty  factor  of  inventory  includes  errors  in  demand  estimates,  variable 
production  yields,  equipment  breakdowns,  strikes,  acts  of  God,  shipping  delays,  and 
unusual  weather  conditions. 

The  economy  factor  of  inventory  let  the  organization  use  cost  reducing 
alternatives  such  as  bulk  purchases  with  quantity  discounts  instead  of  ordering  separately 
without  regard  to  transportation  and  lot  size  economies. 

If  we  want  to  generalize  what  we  covered  so  far,  we  find  out  that  to  determine  the 
number  of  items  on  hand  at  a  given  time  is  a  challenging  task.  Even  though  we  are 
dealing  with  one  type  of  inventory  within  the  organization  such  as  finished  goods  in 
retailer,  we  have  four  functions  of  it  to  consider.  In  our  thesis  problem,  we  are  focusing 
on  the  time  and  uncertainty  function. 

Functional  classification  is  another  way  to  look  at  the  inventory  itself.  It  shows  us 
a  different  aspect  of  inventory  to  consider.  On  the  way  to  answering  the  reason  of 
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existence  of  inventory,  we  covered  the  four  function  of  it  so  far.  We  will  touch  on  the 
classification  aspect  of  inventory  in  the  following  section. 

Functional  classification 

We  can  place  all  inventories  in  one  or  more  of  the  following  based  on  its  utility  as 
Tersine  states  and  explains  six  functional  classifications. 

1 .  Working  stock 

2.  Safety  stock 

3.  Anticipation  stock 

4.  Pipeline  stock 

5.  Decoupling  stock 

6.  Psychic  stock  (Tersine,  1994:7) 

Working  stock  is  inventory  acquired  and  held  in  advance  of  requirements  so  that 
ordering  can  be  done  on  a  lot  size  rather  than  on  an  as  needed  basis. 

Safety  stock  is  inventory  held  in  reserve  to  protect  against  the  uncertainties  of 
supply  and  demand. 

Anticipation  stock  is  inventory  built  up  to  cope  with  peak  seasonal  demand, 
erratic  requirements  (strikes,  vacation  shutdowns),  or  deficiencies  in  production  capacity. 

Pipeline  stock  is  inventory  on  trucks,  ships,  and  railcars  or  in  a  literal  pipeline 
externally.  On  the  other  hand,  if  it  is  being  processed,  waiting  to  be  processed,  or  being 
moved  internally,  then  it  corresponds  to  pipeline  internally. 

Decoupling  stock  acts  as  lubrication  for  the  supply-production-distribution  system 
that  protects  it  against  excessive  friction.  Remember  that  we  can  treat  each  system 
individually. 
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Psychic  stock  increases  the  chance  an  item  is  seen  and  considered  for  purchase. 
Full  shelves  increase  sales  whereas  understock  shelves  as  well  as  stockouts  can  lead  to 
lost  sales  and  lost  customers. 

Looking  at  inventory  in  terms  of  functional  classification  can  let  us  focus  on 
specifically  one  or  more  inventory  class.  We  should  manage  every  element  of  functional 
classification  of  inventory  in  such  a  different  way  that  we  should  take  consideration  into 
their  specific  attributes.  We  can  say  that  existence  of  inventory  is  based  upon  functional 
factors  and  functional  classifications.  We  are  dealing  with  safety  stock  and  working 
stock  classifications  in  the  inventory  side  of  the  thesis  problem. 

After  we  figure  out  the  existence  of  inventory,  we  may  question  ourselves  in  such 
a  way  that  what  kinds  of  problems  we  will  encounter  while  determining  the  numbers 
items  on  hand  at  a  given  time.  We  will  touch  on  the  inventory  problem  classifications  in 
the  following  section. 

Inventory  problem  classification 

We  classified  the  DRAI  problems  in  seven  categories.  Here,  we  will  focus  on  just 
inventory  problems.  By  putting  side  by  side  these  two  different  categories,  we  will 
mainly  explain  the  areas  related  to  only  inventory  problems. 

Tersine  classify  the  inventory  problems  as  follows.  The  inventory  problems  can 
be  organized  according  to  the  repetitiveness  of  inventory  decision,  the  source  of  supply, 
the  knowledge  of  demand,  the  knowledge  of  the  lead-time,  and  the  type  of  inventory 
system  (Tersine,  1994,  9-11). 
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According  to  DRAI  problem  classification,  we  cover  the  repetitiveness  under  the 
decision  type  section,  the  knowledge  of  demand  under  the  demand  section.  One  point 
that  we  did  not  consider  about  knowledge  of  demand  in  DRAI  is  independent  and 
dependent  demand.  Demand  for  finished  goods  are  characteristically  independent  while 
demand  for  raw  materials,  components  and  subassemblies  are  dependent. 

The  other  classification  elements  are  specific  to  inventory  problems.  We  have 
two  kinds  of  source  of  supply  such  as  outside  and  inside  supply.  One  unit  of  an 
organization  supplies  the  items  to  the  other  unit  within  inside  supply.  On  the  other  hand, 
outside  resource  supplies  the  items  with  a  corresponding  purchase  amount  within  outside 
supply. 

We  can  see  the  knowledge  of  the  lead-time  as  a  constant  or  variable.  If  it  is 
variable,  we  may  determine  its  distribution  empirically  or  specified. 

The  last  classification  is  about  the  type  of  inventory  system.  Some  of  the  most 
common  inventory  systems  are  perpetual,  periodic,  material  requirements  planning, 
distribution  requirements  planning  and  single  order  quantity  inventory  system  as  Tesine 
explains  (Tersine,  1994:11). 

Perputual  inventory  system  orders  stock  whenever  the  inventory  position  reaches 
the  reorder  point. 

The  periodic  inventory  system  orders  stock  at  discrete  points  in  time. 

The  material  requirement  planning  (MRP)  system  orders  stock  only  to  meet  pre¬ 
planned  production  requirements.  Anderson  et  al.  define  it  as  follows.  It  is  a 
computerized  inventory  system  whose  functions  are  to  schedule  production  and  to  control 
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the  level  of  inventory  for  components  with  dependent  demand  (Anderson  et  al, 

1994:504). 

The  distribution  requirements  planning  (DRP)  system  orders  stock  to  meet 
distribution  center  requirements  in  multi-echelon  network. 

The  single  order  quantity  system  orders  stock  to  meet  unique  or  short-lived 
requirements. 

Relation  between  Inventory  Problems  and  Newsboy  Problems 

As  Winston  pointed  out  the  objective  of  basic  inventory  models  is  to  minimize  the 

costs  associated  with  maintaining  inventory  and  meeting  customer  demand  by  answering 

the  following  questions.  When  should  an  order  be  placed  for  a  product?  How  large 

should  each  order  be?  The  answer  for  the  first  question  is  related  to  demand  function  as 

well  as  shortage  and  surplus  cost.  The  answer  for  the  second  question  is  related  to 

vehicle  capacity  that  carries  the  cargo  from  main  depot  to  demand  points,  and  base,  depot 

capacity  that  keeps  the  cargo  associated  with  a  holding  cost.  On  the  other  hand,  the  best 

way  to  comprehend  shortage  and  surplus  cost  is  to  examine  the  well-known  news  vendor 

problem  (NV).  Gallego  et  al.  defines  the  classical  NV  as  follows. 

The  newsboy  problem  is  to  decide  the  stock  quantity  of  an  item  when 
there  is  a  single  purchasing  opportunity  before  the  start  of  the  selling 
period  and  the  demand  for  the  item  is  random.  The  trade-off  is  between 
the  risk  of  overstocking  (forcing  disposal  below  the  unit  purchasing  cost) 
and  the  risk  of  understocking  (losing  the  opportunity  of  making  a  profit) 
(Gallego,  1993:  826). 

In  other  words,  we  will  have  a  shortage  cost  if  the  random  demand  is  greater  than 
order  quantity.  Alternatively,  we  will  have  a  surplus  cost  if  the  random  demand  is  less 
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than  order  quantity.  Khouja  states  out  relation  between  inventory  and  newsboy  problems 
as  follows.  The  single  period  problem  (SPP),  also  known  as  newsboy  or  news-vendor 
problem,  is  to  find  the  order  quantity  which  maximizes  the  expected  profit,  or  minimizing 
the  expected  costs  of  overestimating  and  underestimating  demand,  in  a  single  period 
stochastic  framework  (Khouja,  1999:537).  Remember  the  objective  is  to  minimize  the 
cost  in  a  single  deterministic  framework.  The  relation  between  NV  and  inventory 
problems  is  always  valid  if  we  are  dealing  with  a  single  period.  One  may  wonder  what 
the  level  of  relation  would  be  if  we  have  multiple  periods.  Even  though  we  used 
uniformly  distributed  demand  in  the  newsboy  inventory  cost  function,  we  want  to  lay 
down  the  possible  extension  out  of  NV.  Since  NV  has  wide  applicability  in  economy,  we 
can  analyze  the  NV  literature  with  more  depth  as  Khouja  does.  The  extensions  to  NV  as 
follows: 


1 .  Extensions  to  different  objectives  and  utility  functions 

2.  Extensions  to  different  news-vendor  pricing  policies  and  discounting 
structures. 

3.  Extensions  to  different  supplier  pricing  policies. 

4.  Extensions  to  random  yields. 

5.  Extensions  to  different  states  of  information  about  demand. 

6.  Extensions  to  constrained  multi-product. 

7.  Extensions  to  multi-product  with  substitution. 

8.  Extensions  to  multi-echelon  system. 

9.  Extensions  to  multi-location  models. 

10.  Extensions  to  models  with  more  than  one  period  to  prepare  for  the 
selling  season. 

11.  Other  extensions  (Khouja,  1999:538). 


We  will  mention  these  extensions  roughly  based  upon  Khouja  literature  review. 
The  ideas  we  presented  under  each  extension  are  based  upon  several  authors  studies 
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covered  by  Khouja.  We  aim  to  give  a  general  idea  where  we  can  extend  the  single  period 
inventory  problems  and  relate  some  of  these  extension  to  our  thesis  problem. 

The  first  extension,  extensions  to  different  objectives  and  utility  functions,  says 
that  maximizing  the  probability  of  achieving  a  target  profit,  denoted  by  Pb,  is  more 
consistent  with  the  actions  of  organizations’  managers  than  maximizing  the  expected 
profit.  Deriving  a  closed  form  for  an  optimal  order  quantity,  denoted  by  Q*,  for  various 
kinds  of  demand  distributions  such  as  Normal,  Beta,  Gama  is  a  matter  of  problem.  The 
key  point  here  Q  provide  us  with  the  objective  of  maximizing  profit  and  achieving  PB. 
The  objective  of  achieving  PB  can  be  extended  under  two  cases.  First  case  is  all  about 
shortage  cost  which  can  be  either  zero  or  greater  than  zero.  Second  case  is  a  NV  problem 
of  two-product  instead  of  one.  Maximizing  expected  utility,  and  analyzing  demand  risk, 
which  can  be  either  risk  averse  or  risk  neutral,  are  another  issue  of  objective  function. 

The  second  extension,  extensions  to  different  news  vendor  pricing  policies  and 
discounting  structures,  is  about  determining  Q*  when  suppliers  offer  quantity  discounts. 
We  can  apply  discount  policy  in  three  different  ways.  All  unit  quantity  discounts  are 
applicable  for  every  item  we  purchased.  On  the  other  hand,  incremental  quantity 
discounts  are  applicable  for  items  purchased  after  break  point.  Third  discount  type  is 
about  the  level  of  loading  items  into  truck.  The  more  we  close  to  the  capacity  of  truck, 
the  more  discounts  we  have. 

The  idea  of  the  larger  container  we  used  in  transportation  instead  of  smaller  ones 
in  NV  the  smaller  units  cost  we  have  is  true  for  the  case  where  Q  is  made  up  of  n 
standard  size  containers.  Remember  what  Stenger  said  about  the  tradeoff  between 
inventory  and  transportation.  We  will  rephrase  his  idea  for  the  completeness.  As  we  do 
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things  to  reduce  transportation  costs  (bigger  shipment  sizes  and/or  slower,  less  reliable 
models),  we  increase  the  need  for  inventories  (cycle  stocks  and  safety  stocks)  at  the 
receiving  point  (and  possibly  at  the  shipping  point)  (Stenger,  1999).  The  key  point  in  this 
context  is  simply  we  do  not  wait  for  demand  for  fill  up  to  the  capacity  of  vehicle  as 
Stenger  implies;  instead,  we  have  a  chance  to  load  the  demanded  item  into  larger  or 
smaller  containers.  After  we  clarify  these  two  ideas  in  order  to  show  that  there  is  no 
contradiction  between  two  claims,  we  can  touch  on  the  situation  where  a  proportion  of 
customers  are  willingly  to  wait  for  emergency  supply  provided  we  have  a  shortage  of 
demand.  This  is  an  another  type  for  this  context. 

The  third  one,  extensions  to  different  supplier  pricing  policies,  explains  the 
relationship  between  customer  demand  and  selling  price  of  an  item.  The  option  can  be 
decrease  the  selling  price  in  order  to  increase  the  demand.  The  other  option  can  be 
multiple  discounts  until  the  organization  sells  excess  inventory. 

The  fourth  one,  extensions  to  random  yields,  considers  the  defective  items  in  Q*. 
By  taking  into  consideration  of  these  defective  items,  we  can  calculate  Q*  under 
uniformly  or  exponentially  distributed  demand. 

The  fifth  one,  extensions  to  different  states  of  information  about  demand,  is  all 
about  parameters  of  demand  function.  It  can  be  either  simple  demand  distribution  such 
as  uniform  or  complex  demand  such  that  no  known  distribution  can  fit  the  demand 
function.  Another  promising  option  is  more  promising.  It  is  sufficient  to  know  the  mean 
and  variance  of  demand  to  be  able  to  derive  Q*  to  maximize  expected  profit. 

The  sixth  one,  extensions  to  constrained  multi-product,  is  focusing  on  multi¬ 
product  instead  of  single  product.  When  we  think  that  a  typical  newsstand  will  have  a 
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large  number  of  products  to  sell,  corresponding  to  n  variables  in  the  model,  whereas  a 
small  number  of  constraints,  corresponding  to  m  constraints  in  the  model.  There  is  a 
possibility  that  we  can  convert  the  primal  models  into  dual  and  solves  the  dual  of  it  with 
fewer  variables. 

The  seventh  one,  extensions  to  multi-product  with  substitution,  is  related  to  the 
following  situation  where  there  is  shortage  of  demanded  item.  When  a  customer 
demands  another  item,  if  substitution  his  demand  with  another  item,  or  goes  to  another 
newsstand  to  meet  his  demand. 

The  eight  one,  extensions  to  multi-echelon  system,  is  valid  for  selecting  optimal 
components  stock  levels  in  an  assemble-to-order  system.  As  Tersine  defines  that  when 
an  order  is  received,  the  product  is  assembled  from  group  standard  subassemblies  already 
in  inventory.  Since  the  product  has  numerous  optional  features  and  each  customer  may 
desire  a  unique  configuration,  finished  goods  are  not  usually  available  at  the  department. 
The  assembly  lead-time  is  usually  quite  short  for  most  products  (Tersine,  1994:24). 

The  ninth  one,  extensions  to  multi-location  models,  propose  multiple  locations  in 
the  problem  setting.  Each  location  has  either  same  selling  season  or  selling  seasons  of 
the  different  location  lag  each  other.  We  may  examine  the  former  case  as  follows,  we 
may  have  centralized  system  which  consists  of  central  depot  and  n  retailer  or 
decentralized  system,  which  consists  of  independent  sellers  that  keep  inventory  at  their 
own  expense.  We  can  examine  the  latter  case  as  follows.  Khouja  gives  an  example  to 
explain  this  case  such  that  the  ending  of  summer  season  in  United  States  of  America  may 
let  the  organization  sell  his  remaining  in  Australia  where  summer  is  about  to  begin. 
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The  tenth  one,  extensions  to  models  with  more  than  one  period  to  prepare  for  the 
selling  season,  is  related  to  selling  items  in  more  than  periods  in  a  single  season.  Here, 
we  assumed  that  we  have  more  than  one  period  in  the  selling  season.  The  problem  we 
have  in  this  type  of  extension  can  be  in  two  fashions.  The  first  one  is  determining  the  Q* 
or  production  quantity  towards  the  ends  of  season.  The  second  one  determining  the 
probability  that  customer  will  buy  the  specified  product  in  the  beginning  of  season.  As 
season  rolls  on,  the  probability  of  customers’  demand  to  that  item  gets  identified. 

The  last  extension  is  not  focusing  on  one  area.  Yield  management  can  be  an 
interest  of  this  section.  Khouja  gives  an  example  to  clarify  this  type  of  management. 
Airlines  offer  discount-fare  tickets  if  you  buy  the  tickets  before  flight  departures  such  as 
two,  three  weeks  before.  What  would  be  the  number  of  seats  available  at  discount  rate 
can  be  the  one  of  the  main  questions  in  the  yield  management.  Another  example  is 
advanced  booking  of  orders.  Buying  seasons  tickets  before  athletics  games  season 
started  is  an  example  for  advanced  booking  of  orders. 

As  a  result  of  these  extensions  we  see  that  NV  is  veiy  applicable  in  both  retailing 
and  pure  service  organizations  such  as  air  transportation.  Also,  the, reduction  in  product 
life  cycles  brought  about  by  technological  advances  makes  the  NV  more  relevant 
(Khouja,  1999:550). 

Without  considering  the  distribution  of  demand  function  itself,  it  seems  that  the 
inventory  problem  can  contribute  a  significant  amount  of  cost  to  total  cost  due  to  various 
kinds  of  costs  like  holding  cost,  and  shortage  cost.  This  result  brings  our  attention  to  a 
point  such  that  we  should  take  care  of  these  costs  in  our  formulation  if  we  want  to  model 
real  world  issues. 
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The  Relationship  between  Inventory  and  Repair  Process 


It  is  clear  that  inventory  plays  a  key  role  in  SCM.  If  we  have  high  tech  equipment 
like  an  aircraft  in  our  inventory,  then  repairing  process  becomes  as  much  important  as  the 
inventory.  Because  aircraft  is  consist  of  many  high  tech,  high  cost  items  (Sherbrook, 
1992:xxi).  Before  going  through  inventory  with  the  repairable  items  we  should  take  a 
look  at  some  more  definitions  such  as  reorder  point,  cycle  time,  and  backorder. 

Anderson  et  al.  states  these  definitions  in  a  clear  way.  Inventory  position  is  what  we  have 
on  hand  plus  the  inventory  on  order.  Reorder  point  is  the  inventory  position  at  which  a 
new  order  should  be  placed.  Cycle  time  is  the  length  of  time  between  the  placing  of  two 
consecutive  orders.  Backorder  is  the  receipt  of  an  order  for  a  product  when  there  are  no 
units  in  inventory  (Anderson,  1994:503-504). 

By  using  the  definitions  above,  it  would  be  easy  to  comprehend  the  foundations  of 
system-based  repairable  item  inventory  models  that  are  currently  used  in  United  States 
Air  Force  (USAF).  These  principles  are  (S-l,  S)  ordering  policy,  Poisson  processes, 
Poisson  random  variable,  Palm’s  Theorem,  and  the  nature  of  the  demands  (Kliger, 
1994:9). 

(S-l,  S)  Ordering  Policies.  If  the  inventory  theory  objective  is  to  compute 
optimal  stock  levels  for  each  item,  then  what  is  the  stock  level  or  inventory  position? 
(Sherbrook,  1992:23).  We  should  define  the  inventory  position  as  follows. 

(1)  S  =  OH  +  DI-BO 

where  stock  level  denoted  by  S,  which  is  a  constant  is  equal  to  the  number  of  units  of 
stock  on  hand  (OH),  plus  the  number  of  units  of  stock  due  in  from  repair  and  re-supply 
minus  the  number  of  backorders  (BO)  (Sherbrook,  1992:23).  Since  repairable  items  tend 
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to  be  high  tech,  high  cost,  and  low  demand  at  a  base,  we  can  define  the  reorder  point  at  a 
fixed  stockage  objective  level  like  S-l.  In  other  words,  as  soon  as  inventory  level  falls 
down  by  S,  we  place  an  order  for  an  equal  number  of  units  that  have  been  demanded 
(Klinger,  1994:10). 

Poisson  Processes. 

We  will  explain  some  definitions  to  be  able  to  understand  Poisson  processes.  The 
definition  and  the  properties  of  counting  processes,  the  definition  of  independent 
increments  are  the  key  elements  on  the  way  to  Poisson  processes.  A  stochastic  process 
[N(t),  t  >0}  is  said  to  be  a  counting  process  if  N(t)  represents  the  total  number  of 
“events”  that  have  occurred  up  to  time  t  (Ross,  1997:249).  The  counting  process  must 
satisfy  the  following  conditions: 

i.  N(t)  >  0. 

ii.  N(t)  is  an  integer  valued. 

iii.  If  s  <  t,  then  N(s)<  N(t). 

iv.  For  s  <  t,  N(t)  -  N(s)  equals  to  number  of  events  that  have  occurred 
in  the  interval  ( s ,  t )  (Ross,  1997:249). 

The  only  missing  part  in  definition  of  Poisson  process  is  independent  increments.  A 
counting  process  is  said  to  posses  independent  incerements  if  the  numbers  of  events 
which  occur  in  disjoint  time  intervals  are  independent  (Ross,  1997:250).  The  number  of 
events  have  occurred  by  time  5,  that  is  N( 5),  must  be  independent  of  the  number  of  the 
events  occurring  between  times  5  and  10,  that  is  iV(10)  -  N( 5),  would  be  a  good  example 
for  independent  increments.  By  combining  these  three  elements  we  would  define 
Poisson  processes  as  follows. 

The  counting  process  {N(t),  t  >0}  is  said  to  be  a  Poisson  process  having 

rate  X.,  X  >  0,  if 

i.  N{ 0)  =  0. 
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ii.  The  process  has  independent  increments. 

iii.  The  number  of  events  in  any  interval  of  length  t  is  Poisson 
distributed  with  mean  Xt.  That  is,  for  all  s,  t  >  0  (Ross,  1997:251). 


(2)  P {N(t  +  s)  -  N(s )  =  n)  =  exp j-  A  t  ^ j,  n  =  0, 1 , 2, . . . 


The  Poisson  Random  Variable. 


A  random  variable  X,  taking  on  one  of  the  values  0,  1, 2. . .  is  said  to  be  a  Poisson 
random  variable  with  parameter  A,  if  for  some  A  >  0, 


(3)  p(i)  =  P{x  =i}  =  exp 


-X* 

i\ 


We  can  use  Poisson  random  variable  as  an  approximation  of  binomial  random  variable 
when  the  binomial  parameter  n  is  large  and  p  is  small  (Ross,  1997:30). 

Compound  Poisson  Processes. 

We  will  define  this  concept  and  then  give  a  complementary  example  for  it. 

A  stochastic  process  {X(t),  t  >0}is  said  to  be  a  compound  Poisson  process 
if  it  can  be  represented  as 

NO) 

(4 )X(t)  =  jjYi,  t>  0 

;=i 

where{ N(t),  t  >0}is  a  Poisson  process,  and  { 7, ,  i  >0}is  a  family  of 
independent  and  identically  distributed  random  variables  which  are  also 
independent  of  \N(t),  t  >0} .  The  random  variable  X(t)  is  said  to  be  a 
compound  Poisson  random  variable  (Ross,  1997:282). 


The  example  for  the  compound  Poisson  processes,  which  are  composed  of  many 

Poisson  distributed  random  variables,  can  be  as  follows. 

Suppose  that  buses  arrive  at  a  sporting  event  in  accordance  with  a  Poisson 
process,  and  suppose  that  the  numbers  of  customers  in  each  bus  are 
assumed  to  be  independent  and  identically  distributed.  Then  {X(t),  t  >0}is 
a  compound  Poisson  process  where  X(t)  denotes  the  number  of  customers 
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who  arrived  by  t.  y,  represents  the  number  of  customers  in  the  /th  bus 
(Ross,  1997:282). 


Palm's  Theorem. 

It  enables  us  to  estimate  the  number  of  items  in  repair  at  the  steady  state 

probability  if  we  know  the  probability  distribution  of  the  demand  process  and  the  mean  of 

the  repair  time  distribution.  The  formal  definition  of  Palm's  Theorem  is  as  follows: 

If  demand  for  an  item  is  a  Poisson  process  with  annual  mean  m  and 
if  the  repair  time  for  each  failed  unit  is  independently  and 
identically  distributed  according  to  any  distribution  with  mean  T 
years,  then  the  steady-state  probability  distribution  for  the  number 
of  units  in  repair  has  Poisson  distribution  with  mean  mT 
(Sherbrook,  1992:21). 

Nature  of  demand. 

Since  we  talked  about  the  demand  in  great  detail  under  the  DRAI  problems,  we 
would  just  say  that  we  assumed  the  demand  for  repair  items  is  independent.  Demands 
can  be  deterministic  or  stochastic  in  repairable  item  inventory  models  as  well. 

Measurement  Criteria. 

We  have  three  measurement  criteria  while  evaluating  repairable  item  inventory 
systems  like  fill  rate,  backorders,  and  availability. 

Fill  rate. 

It  is  the  percentage  of  demands  that  can  be  met  at  time  they  are  placed 
(Sherbrook,  1992:24).  The  other  definition  for  it  is  the  probability  that  at  least  one  spare 
item  is  available  on  the  warehouse  shelf  when  a  demand  for  an  item  occurs;  it  is  the 
probability  that  the  number  of  demands  during  the  re-supply  time  are  strictly  less  than  the 
spare  stock  level  (Klinger,  1994:13). 
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Backorders. 


They_are  the  number  of  unfilled  demands  that  exist  at  a  point  in  time  (Sherbrook, 
1992:24).  They  are  the  numbers  of  “holes”  in  an  aircraft,  or  the  numbers  of  missing 
items  on  an  aircraft  (Klinger,  1994:13). 


Aircraft  availability  rate  is  the  percentage  of  aircraft,  which  are  available  or  fully  mission 


capable  (Klinger,  1994:13). 


Integrated  Logistics 

Integrated  logistic  idea  includes  the  idea  of  combining  vehicle  routing  and 
inventory.  The  reason  that  we  mention  supply  chain  management  in  the  literature  review 
is  to  understand  the  importance  of  looking  at  individual  processes  together.  After  we 
present  three  main  areas  like  vehicle  routing,  inventory  and  repair,  we  should  come  up 
with  the  following  figure. 


Figure  1:  The  supply  chain  management 
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When  we  look  at  the  graph  we  see  that  inventory  and  information  flow  are 
passing  through  all  three  components.  To  be  fully  effective  in  today’s  competitive 
environment  where  a  necessary  condition  for  success  is  to  have  integration  between 
internal  functions,  we  should  expand  their  integrated  behavior  to  incorporate  customers 
and  suppliers  (Bowersox,  1996:34).  Inventory  flow  helps  us  to  determine  how  many 
items,  when  and  where  needed  (Bowersox,  1996:34).  Information  flow  is  composed  of 
two  parts:  coordination  flows  and  operational  flows.  The  former  one  includes 
coordination  between  strategic  objectives,  capacity  constraints,  logistical  requirements, 
inventory  deployment,  manufacturing  requirements,  procurement  requirements,  and 
forecasting.  The  latter  one  includes  coordination  between  order  management,  order 
processing,  distribution  operations,  inventory  management,  transportation  and  shipping, 
and  procurement  (Bowersox,  1996:37). 

Physical  distribution  is  concerning  about  movements  of  a  finished  goods.  The 
overall  concern  of  manufacturing  support  is  not  how  production  occurs  but  rather  what, 
when,  and  where  products  will  be  manufactured  (Bowersox,  1996:35).  We  are  dealing 
with  materials  that  are  under  the  control  of  the  manufacturing  enterprise.  The  term 
procurement  in  government  sector,  purchasing  in  acquisition,  or  buying  in  general  is 
concerns  with  purchasing  and  arranging  in-bound  movement  of  materials,  parts,  and/or 
finished  inventory  from  suppliers  to  manufacturing  or  assembly  plants,  warehouses,  or 
retail  stores  (Bowersox,  1996:35). 

We  have  known  that  integrated  logistics  approach  will  yield  a  lower 
objective  cost  value.  We  have  already  known  that  Federgruen  and  Zipkin  proved  that 
combining  the  inventory  and  vehicle  routing  problems  together  would  yield  a  lower 
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objective  cost  value.  What  we  are  trying  to  accomplish  in  Chapter  Three  is  to  show 
combining  the  vehicle  routing  and  inventory  in  the  new  problem  definition,  which  is 
extension  of  Federgruen  and  Zipkin  model,  will  yield  weather  a  lower  objective  cost 
value  or  not. 
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Chapter  3:  Methodology 


Introduction 

The  starting  point  of  our  thesis  process  is  the  problem  (Federgruen  and  Zipkin, 
1983:1020)  which  was  presented  in  the  literature  review.  First,  we  will  list  the 
parameters,  variables,  and  sets  of  notations  used  for  inventory/allocation  and  VRPs. 
Second,  we  will  present  an  instance  of  this  problem  (Chan,  1999  draft:9-9).  Third,  we 
will  show  how  to  solve  the  starting  problem  either  by  using  generalized  Benders 
decomposition  method  or  without  using  it.  After  we  solved  it  without  using  generalized 
Benders  decomposition,  we  will  talk  about  the  extension  from  the  initial  model  including 
new  variables,  new  parameters  and  set  of  notations.  We  will  try  to  apply  the  generalized 
Benders’  decomposition  technique  to  the  thesis  problem  in  the  last  phase. 

A  Starting  Point:  VRP  and  Allocation/Inventory  Problem 

We  will  list  the  common  variables  and  sets  of  notations  to  be  used  for  formulating 
inventory  and  vehicle  routing  problems  together  as  follows  (Chan,  1999  draft:9-7): 
xhij is  a  binary  variable  equal  to  1  if  aircraft  with  tail  number  h  flies  the  arc  from  i  to  j,  and 
equal  to  0  otherwise. 

Z,  is  a  continuous  variable  that  denotes  the  amount  of  allocation  among  available 
products  to  node  i. 

y  i  is  a  binary  variable  equal  to  1  if  aircraft  with  tail  number  h  allocates  to  node  i  and 
equal  to  0  otherwise. 

£is  a  random  variable  associated  with  a  demand  with  any  kinds  of  distribution. 
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H  is  the  set  of  all  aircraft  where  h  =  1, •••,//"  and  H"  =  2. 


I  is  the  set  of  all  nodes  where  I  =0,  1 ,  2,  and  3.  The  main  depot  is  node  0. 

Additionally,  we  need  to  know  the  following  parameters: 
djj  is  the  distance  (in  time  units)  between  nodes  i  and  j. 

/?,  is  an  initial  inventory  level  of  node  i. 

c,  is  a  unit  carrying-cost  or  surplus  cost  function  in  location  i. 

C,  shortage  cost  in  location  i. 

Vh  is  the  vehicle  capacity  constraint  associated  with  the  aircraft  with  the  tail  number  h. 
P  is  the  amount  of  available  resources  at  the  central  depot. 

Chan  lays  down  the  model  formulation  explicitly  as  follows: 

The  objective  function 

(5)  Minimize  4  +Z$i(Zi) 

iel  jel  h^H  i 

The  objective  function  is  to  minimize  transportation  cost  and  inventory  costs.  The 
inventory  cost  function  is: 

oo 

(6)  q, (z, )  =  f  CM- P,  - z,  Wt (£)  +  f  Ci (P,  +  z,  ~ QdFt i  e  I 

Pi+z,  0 

where  Fj(.)  is  the  cumulative  distribution  function  of  demand  occurred  at  node  i. 
We  provide  an  example  for  the  usage  of  the  inventory  cost  function  in  the  following 
sections  in  Chapter  Three. 
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We  will  first  lay  down  the  sets  of  constraints  found  in  the  classical  VRP  and  then, 
incorporate  the  additional  constraints  related  to  inventory  in  order  to  combine  VRP  and 
inventory/allocation  problem.  Here  are  the  classical  VRP  constraints: 


|#iifi  =  0 
1  if  i  =  1, 2,...,|l| 


(8)  SExJ 

iel  heH' 


jMif  j=0 
tl  if  j  =  lA...J 


Constraint  (7)  and  constraint  (8)  ensures  that  we  used  all  available  vehicles  stationed 
at  depot  to  serve  the  demand  points. 

(9)  ZxJ-  Sx‘=0  Vh,  Vp  e  I. 

ieMp  jeMp 

Constraint  (9)  maintains  the  route  continuity  for  all  vehicle  types  where  Mp 
denotes  the  out-of-nodes  (source  nodes)  and  M  denotes  the  into-nodes(sink  nodes). 

Note  that  Mp  and  M  p  are  the  same  sets  of  nodes  for  each  constraint. 

(10)  E  <  1  Vh,  and  Ex-o  <  1  Vh 

Constraint  (10)  consists  of  two  constraints.  They  ensure  that  vehicle  availability 
is  not  exceeded  at  depot  0  for  vehicle  h  where  M0  denotes  the  out-of-node  0  and 
M0  denotes  the  into-node  0. 

(11)  EE*J<  \L\~l  Lc{2,...,|/|},  2<|L|<|/|-1;  h  e  H 

teL  jeL 

Constraint  (1 1)  is  the  sub-tour  breaking  constraint  where  L  is  the  non-empty  set  of 
1-1.  The  reason  we  are  subtracting  one  out  of  the  number  of  nodes  in  the  problem  is  due 
to  an  having  one  central  depot.  All  vehicles  begin  their  tour  from  the  main  depot  and 
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come  back  to  it.  This  ensures  that  all  tours  in  the  solution  be  legitimate;  namely,  no  cycle 
will  occur  in  the  solution. 

The  constraints  related  to  allocation  (inventory)  side  of  the  problem  are  now 
presented. 

(12)  Zz,5>*<V,  heH 

iel 

Constraint  (12)  ensures  that  the  load  assigned  to  each  vehicle  h  is  within  its  capacity. 

(13)  Zz,.  <F 

iel 

Constraint  (13)  guarantees  that  the  total  amount  shipped  is  available  at  the  depot. 

(14)  2*5  ~y?*°  v/>  heH  (Chan,  1999 draft:9-7). 

jeMj 

Constraint  (14)  satisfies  the  linking  constraints  between  vehicle  routing  variable  xhy 
and  allocation  variable  y“j.  If  we  place  an  order  to  node  i,  one  arc  out  of  node  i  should 
equal  to  one,  otherwise  zero. 


A  Numerical  Example  for  the  Starting  Model 

We  mentioned  in  literature  review  that  Federgruen  and  Zipkin  combined  the  VRP 
and  IRP  in  their  study.  For  the  sake  of  completeness,  we  would  like  to  present  the 
problem  description  again. 

The  initial  inventory  (perhaps  supply  remaining  from  the  previous  day)  for 
each  location  is  reported  to  the  depot.  This  information  is  used  to 
determine  for  the  following  day  the  allocation  of  the  available  product 
among  the  locations.  The  assignment  of  locations  to  the  vehicles  and  the 
routes  are  set  at  the  same  time.  After  deliveries  are  made  (say  at  the  end 
of  day)  the  demands  occur,  and  inventory-carrying  and  shortage  costs  are 
incurred  at  each  location  proportional  to  the  end-of-the-day  inventory 
levels.  (Federgruen,  Zipkin,  1984:1020). 
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Johnson  et  al.  and  Moore  et  al.  presented  a  four-node  example  for  the  problem  of 
Federgruen  and  Zipkin  which  is  covered  in  great  detail  in  Chan’s  book  (Chan,  1999 
Draft:9-9).  The  problem  is  a  combination  of  the  shortest  routing  problems  for  multiple 
vehicles  and  the  allocation/inventory  problems. 

The  objective  is  to  minimize  the  inventory  costs  incurring  at  each  demand  point 
and  the  travel  costs.  Inventory  cost  is  the  newsboy  inventory  cost  where  there  are 
shortage  and  surplus  costs.  Travel  cost  is  related  to  the  distance  traveled  between  nodes. 

The  objective  function  above  is  subject  to  the  following  constraints.  There  are 
maximum  and  minimum  inventory  levels  for  each  demand  point.  We  cannot  deliver 
more  than  depot  capacity.  Available  supply  amount  at  the  depot  and  the  vehicle  capacity 
constraint  are  additional  limitations  on  the  problem.  We  should  assume  that  the  vehicle 
capacity  is  not  exceeded  during  the  period.  We  are  considering  delivering  only  one  type 
of  supply  materials  to  demand  points,  but  not  delivering  and  picking-up  at  the  same  time. 
Besides  that  we  can  visit  each  demand  point  with  at  most  one  vehicle  without  creating  a 
cycle.  We  have  two  real  vehicles  and  one  dummy  vehicle  in  the  model.  We  used  the 
dummy  vehicle  for  two  reasons.  The  first  reason  is  to  satisfy  the  TSP  requirement,  which 
is  visiting  each  node  exactly  once.  The  second  reason  is  to  prohibit  delivering  to  nodes 
which  do  not  need  any  replenishments.  In  short,  if  the  demand  point  does  not  need  any 
replenishment,  we  send  a  dummy  vehicle  to  that  node  so  as  to  visit  but  not  deliver.  After 
the  demand  points  declare  the  amounts  at  hand,  the  central  depot  prepares  a  delivery  plan 
to  minimize  the  objective  cost,  that  is,  to  minimize  the  inventory  and  traveling  cost  for 
one  period,  which  can  be  a  day,  a  week,  or  a  month. 
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Figure  2  -  A  Starting  Model  -  Single  depot/multiple  route/tour  network 


Node  0  denotes  the  main  depot,  and  node  1,  2,  and  3  denote  the  demand  points  in 
the  basic  network  as  shown  above. 

We  displayed  the  cost  matrix  between  each  demand  points  and  inventory 
allocation  information  in  Table  3  and  4,  respectively. 


Table  2:  Traveling  cost  matrix 


Node  0 

Node  1 

Node  2 

Node  3 

Node  0 

- 

3 

2 

Node  1 

3 

- 

4 

5 

Node  2 

3.5 

4 

3.8 

Node  3 

2 

5 

3.8 

The  traveling  cost  matrix  is  symmetric.  Going  from  node  i  to  node  j  is  the  as 
going  from  node  j  to  node  i.  It  can  be  interpreted  as  the  time  required  from  getting  from 
demand  point  1  to  2,  or  the  cost  associated  with  that  traveling  time.  It  is  not  a  distance  in 
terms  of  mileage,  however,  the  distance  is  in  terms  of  time.  Even  though  the  distance 


between  two  nodes,  let  us  say  node  i  to  node  j,  are  the  same,  it  takes  different  time  units 
to  get  from  node  i  to  node  j  due  to  weather,  air  traffic,  cargo  load  etc.  Therefore,  we 
preferred  using  time  instead  of  mileage  in  the  distance  matrix. 


Table  3  -  Inventory-allocation  data 


Node 

1 

2 

3 

Max.  Inventory 

500 

100 

500 

Min.  Inventory 

0 

0 

0 

Initial  Inventory 

100 

50 

150 

Demand  (PDF) 

0.002 

0.010 

0.002 

Newsboy  cost  function 

2600 -8zl 

437.5 -7.5z2 

3012.5  -  6.5z3 

On  the  other  hand,  we  calculate  the  inventory  cost  function  according  to  Uniform 
distribution.  All  the  information  we  used  in  order  to  set  up  these  two  tables  above  is 
corresponding  to  parameter  settings  except  demand  in  terms  of  probability  distribution 
function  (pdf)  and  newsboy  cost  function  in  the  inventory/allocation  data  table.  Note  that 
the  inventory  cost  functions  is  a  nonincreasing  functions.  It  means  that  the  more  we 
deliver  provided  that  the  depot  capacity  is  available,  the  less  we  would  have  inventory 
cost.  The  second  important  observation  in  the  inventory  table  is  that  there  is  no  inventory 
cost  function  associated  with  main  depot.  Because  there  is  no  demand  associated  with 
main  depot.  The  role  of  the  main  depot  is  to  satisfy  the  demand.  The  third  important 
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observation  is  that  there  is  no  lateral  re-supply  between  nodes.  It  means  that  we  can  not 
supply  node  1  from  any  demand  nodes  but  only  main  depot. 

A  Look  to  the  Inventory  Cost  Function  at  Demand  Points 

We  would  like  to  explain  how  we  get  the  inventory  cost  functions  (i.e.  2600  -  8z0 
for  demand  point  1  by  giving  an  example. 

(15)  qi(zi)=  J  Ci(t-\li-zi)dFi({)  + 

Pi+z,  0 

Let  us  Ci=$20  per  item,  which  is  shortage  cost,  ci=$5  per  item,  which  is  unit  carrying 
cost,  /?,  =1  unit,  which  is  the  initial  level  that  is  reported  to  main  depot,  and  the  depot 

maximum  capacity  is  equal  to  4  unit.  By  taking  the  integral  of  the  inventory  cost 
function  with  the  given  parameters,  we  would  come  up  with  the  following  result,  that  is 
12.8-3.75zj. 

After  we  understand  how  we  develop  the  inventory  cost  functions,  we  can  discuss 
the  pdfs  of  the  demands.  We  applied  the  pdf  formulation  of  Uniformly  distributed 
demand  (Wackerly,  Mendenhall,  Scheaffer,  1996:150)  which  can  be  equally  likely 
between  zero  and  500  unit  for  demand  point  1,  which  corresponds  to  0.002  in  Table  4. 

The  Solution  Approach  to  the  Starting  Point 

We  have  two  alternatives  to  solve  the  initial  problem.  One  of  them  is  to  use 
generalized  benders  decomposition  method,  and  the  other  one  is  to  formulate  and  solve 
the  problem  without  considering  any  issues.  In  the  latter  way,  we  let  the  solver  pick  the 


AjtZf 

c(P,.  +  7.-#MF;(0,  iel 
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appropriate  switch  on/off  variables,  yhj,  that  correspond  to  a  placing  order  decision 
variables. 

Generalized  Benders  Decomposition. 

Benders  partitioning  procedure  can  be  applied  for  solving  mathematical 
programming  problems  that  includes  either  integer,  linear,  a  nonlinear  in  terms  of 
variables  or  functions. 

(16)  Min  /(*>  y)  st  G(x ’  y)  -  °»  x  e  X,  y  e  Y 

*.y 

where  the  functions/and  g  are  differentiable  function  (Sahinidis,  Grossman,  1991:481). 
We  can  apply  this  technique  for  the  following  situations  (Geoffrion,  1972): 

For  fixed  y,  (16)  separates  into  a  number  of  independent  optimization  problems, 
each  of  them  has  their  own  solution  vector  x.  Multi-period  design  problems  can  be  an 
example  to  this  situation. 

For  fixed  y,  (16)  yields  a  special  structure  that  we  can  solve  it  efficiently.  Mixed 
integer  linear  (MIP)  and  mixed  integer  nonlinear  problems  (MINP)  can  be  for  this 
condition. 

Problem  (16)  is  not  in  a  convex  region  if  we  keep  x  and  y  all  together  in  the 
model.  But  fixing  y  can  lead  us  to  have  a  convex  set. 

If  you  remember  the  starting  point  problem  you  will  notice  that  there  is  a  binary 

u 

variable,  y  „  such  that  we  can  fill  the  demand  coming  from  the  node  if  it  is  equal  to  one 
otherwise  we  cannot  fulfill  the  demand  if  it  is  equal  to  zero.  We  realized  that  if  we  fix 
the  placing  order  constraints,  we  would  have  a  problem  with  a  special  structure  that  we 
can  deal  with  a  known  algorithm. 


56 


Before  going  further,  we  will  explain  what  we  understand  about  fixing  the  placing 
order  constraints.  We  know  that  we  can  visit  the  nodes  at  most  one  vehicle  from  problem 
assumptions.  We  have  three  nodes  and  three  vehicles  in  the  starting  model,  and  each 
vehicle  can  visit  each  node.  It  means  we  can  choose  any  vehicle  to  visit  any  node.  For 
example,  vehicle  0,  which  is  a  dummy  vehicle  has  a  chance  to  visit  node  1,  2,  and  3.  This 
is  true  for  the  other  vehicle  1,  and  2  as  well.  It  means  we  have  three  combinations  for 
each  node  to  visit,  that  corresponds  to  the  combination  of  27.  Some  of  the  combinations 
are  certainly  infeasible  due  to  the  fact  that  we  should  send  two  real  vehicles  out  of  and 
into  the  main  depot. 

As  Geoffrion  stated,  we  would  have  a  MIP  by  fixing  feasible  y’s  such  that  we 
would  have  two  subproblems  to  solve  in  our  starting  model.  These  are  a  traveling 
salesmen  problem  with  k  persons  (k-TSP),  and  an  inventory/allocation  problem. 

The  original  model  written  by  Chan  states: 

Recall  that  xh=  Xjjh  and  y11  =  yj  are  vectors  of  decision  variables  for  routing  and  coverage 
respectively  (Chan,  1999  Draft:9-8). 

MinSESdux5 

;  j  h 

subject  to 

(17)  xheX,  yeY,heH 

(18)  ^h(xh,yh)  >  0 

(19)  2>iy^Vh  h  e  H 

iel 

(20)  £Zi<P 

iel 

where  $  represents  all  linear  inequalities  defining  the  hth  TSP  polytope.  The  set  Y11 
represents  all  the  possible  assignment  of  demands  i  to  vehicle  h  provided 
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(21)  Zy-=  1  for  Vh. 

iel 

Having  presented  the  original  model,  we  are  ready  to  discuss  the  master  problem, 
which  is  equivalent  to  the  original  model  in  y-space.  The  master  can  be  stated  as  follows. 

(22)  Min  {z|  z  >  z'(y) ,  y  e  Y"} 

We  should  define  the  master  Benders  cut  and  the  constraints  associated  with 
ye  Y". 

We  obtained  the  Benders  cut  after  a  couple  iterations  in  the  following  form.  For 
those  who  are  interested  in  the  derivation  stages,  (Federgruen  and  Zipkin  1983:  1020), 
and  (Chan,  1999  Draft,  9-8)  are  good  references. 

(23)  z>QP+Z( phvh  +  Kh  +  *£ )  +  X  P/°y?  +  Z  I(*f  +  (pi  )y.h • 

heH  iel- 0  heHiel- 0 

Q  is  the  dual  variable  corresponding  to  the  main  depot  constraint,  and  ph  is  the  dual 
variable  to  account  for  the  hth  delivery- vehicle  capacity  constraints  in  the 
inventory/allocation  supbproblem.  Kh  is  the  total  travel  cost  on  tour  h,  and  Khj  is  the  unit 
arc  cost  (or  marginal  cost)  to  reaching  node  i. 

(phi  is  the  inventory  (newsboy )  cost  in  the  form  of  Eq.(24). 

(24)  q  (Zi)-(n-ph).z(. 

These  dual  variables  are  playing  a  major  role  while  forming  the  Benders  cut. 

Both  of  them  have  an  impact  on  the  intercept  and  the  slope  of  the  Benders  cut.  If  O  is 
greater  than  zero,  then  the  intercept  point  will  be  higher  but  slope  will  decrease.  It  means 
that  we  do  not  have  a  very  deep  Benders  cut.  Deep  Benders  cut  means  we  divide  the 
optimal  region  in  such  a  way  that  we  become  so  close  to  the  optimal  solution.  We  can 
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think  of  the  same  idea  for  p.  If  p  is  greater  than  zero,  then  the  intercept  point  will  be 
higher  but  the  slope  will  decrease.  In  our  numerical  example,  we  reached  the  optimal 
solution  within  three  iterations.  We  can  perform  a  sensitivity  analysis  as  well.  For  the 
first  iteration  the  Q.  is  equal  to  zero,  p  associated  with  dummy  vehicle  and  vehicle  1  is  8, 
and  6.5,  respectively.  It  means  if  we  increase  the  capacity  of  vehicle  one  by  one  unit,  we 
improve  the  objective  function  by  6.5  unit.  Larger  vehicle  capacity  denotes  higher 
intercept,  but  not  steep  slope.  It  means  we  can  cut  the  optimal  region  mildly  so  that  we 
can  reach  the  optimal  solution.  In  other  words,  it  may  take  a  couple  of  more  iteration  to 
get  the  optimal  solution. 

If  both  of  them  is  zero,  then  we  can  generate  many  cuts  without  reducing  the  size 
of  the  optimal  region.  Since  we  have  only  operating  cost  to  form  up  the  Benders  cut,  we 
cannot  generate  a  useful  cut  to  reduce  the  feasible  region.  It  means  we  cannot  move 
toward  the  optimal  solution.  We  stay  at  the  same  point  in  the  same  region,  but  the  y^’s 
are  changing. 

After  we  define  the  Benders  cut  in  the  master  problem  we  are  ready  to  explain 
two  sets  of  constraints  associated  with  y  e  Y" .  The  first  set  of  constraint  is  the  Eq.(21) 
that  satisfies  the  condition  that  each  node  is  visited  by  exactly  one  vehicle.  The  second 
set  of  constraints  is  equation  (19)  that  satisfies  the  condition  that  we  allocated  vehicles 
within  their  capacities.  Here  we  should  note  that  since  we  are  solving  the  master  problem 
in  terms  of  y’s,  we  must  substitute  the  values  of  z,  ‘s  found  in  the  previous 
allocation/inventory  subproblem  into  the  Eq.(19).  Solving  the  master  problem  in  y  give 
us  a  new  vector  of  y’s  to  pass  to  k-TSP  and  inventory/allocation  subproblem.  We  redo 
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the  procedure  until  we  generate  no  longer  cuts  for  the  master  problem.  The  stopping 
criterion  holds  if  we  reach  the  point  where  the  two  successive  iterations  are  all  the  same. 

An  example  for  the  application  of  this  procedure  will  help  to  clarify  it.  Let  an 
initial  feasible  yh;  be  y°1=y22=y13=l  and  the  rest  of  yYs  equal  to  zero.  When  we  solve  the 
k-TSP  problem  yielding  the  following  result:  x'o3  =x13o=l,  and  x2o2=x22o=l,  and 
x°oi=x°io=l  with  an  objective  function  value  of  17.  The  formulation  and  solution  to  the 
TSP  sub-problem  is  at  Appendix  1  and  2,  respectively.  By  substituting  the  fixed  yYs 
into  the  delivery  allocation  constraints  (Eq.19),  we  solve  the  inventory/allocation 
subproblem  in  terms  of  zfs  that  yields  zi=0,  Z2=50,  z3=150  with  an  objective  function 
value  4700.  The  formulation  and  the  solution  to  the  Inventory/allocation  sub-problem  is 
at  Appendix  3,  and  4,  respectively.  The  necessary  values  to  set  up  the  Benders  cut  are  as 
follows.  The  dual  variables  associated  with  main  depot  supply  constraint  and 
delivery/allocation  constraints  for  each  vehicle  is  Q=0,  p'=8,  p2=6.5,  and  p3=0, 
respectively.  k°=6,  k‘=4,  k2=7  are  the  operating  cost  for  each  vehicle.  For  example, 
vehicle  0  goes  to  node  0  and  comes  back,  therefore  its  operating  cost  is  3+3=6.  We  can 
calculate  the  other  operating  costs  in  the  same  manner.  k°i=3  t  k*3=2  _  and  k22  =3.5  is  the 
marginal  cost  servicing  the  node  1,  2,  and  3,  respectively.  We  give  another  example  to 
explain  this  situation.  If  Vehicle  0  goes  to  node  1,  its  marginal  cost  to  servicing  node  1  is 
the  cost  where  vehicle  0  comes  from.  In  this  solution  set,  k°i=3.  But  if  vehicle  0  visits 
node  1  after  visiting  node  2,  then  the  marginal  cost  would  be  equal  to  the  cost  of  going 
from  node  2  to  node  1,  that  is  k°i=4.  By  using  the  results  of  the  two  sub-problems,  the 
Benders  cut  is 
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(25)  z  -  5203 *y°i  +  675*y°2  -  1675*y°3  -  2600*y\  +  262.5*yJ2  -  1065*y13 

-  2600*y2i  -  66*y22  -  2038*y23  >  1000.5 

and  the  yhi  constraints  are 

(26)  50*y°2  +  150*y°3  <  0 
50*y12  +  150*yJ3  <  150 
50*y22  +  150*y23  <  200 

(27)  y°i  +  yS  +  y2!  =  1 
y°2  +  yli  +  y22  =  l 

y°3  _j_  -|-  y23  —  \ 

where  y11;  is  a  binary  variable,  Eq  (26)  is  the  vehicle  capacity  constraints,  and  Eq.(27)  is 
the  placing  an  order  constraints  for  each  demand  point.  The  calculation  of  Benders  cut  is 
at  Appendix  IF.  The  master  problem  above  yields  the  new  set  of  y^  such  that  z=4731, 
and  y’3=  y2i=  y22=l-  The  formulation  and  the  solution  to  the  relaxed  master  problem  is  at 
Appendix  5  and  6,  respectively.  We  continue  to  solve  the  k-TSP  and  inventory/allocation 
problem  with  these  yh (s.  We  generate  a  new  cut  for  the  relaxed  master  problem  and  this 
process  continues  until  the  two  successive  iterations  yields  the  same  yi  ‘s  and  the  same 
cut  for  the  master  problem.  Table  4  shows  the  result  of  each  iteration. 

Table  4:  The  results  of  Generalized  Benders  Algorithm 


7 

Zi 

operating 

cost 

inventory 

cost 

total  cost 

Y0,=l 

y22=i 

Y,3=l 

Z1=0 

Z2=50 

Z3=150 

17 

4717 

Yi3=1 

y2,=i 

ES5I 

Z3=150 

14.5 

3475 

y12=i 

Yi3=1 

EB 53 

Z3=100 

15.3 

3425 
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We  may  think  of  the  solution  process  in  a  different  manner.  After  we  have  fixed 
then's  in  the  model,  we  will  have  a  MIP.  We  may  get  a  solution  in  terms  of  xhjj  and  zfs 
out  of  MIP  while  fixing  the  y*Y  s.  By  the  same  token  as  we  did  in  the  previous  solution 
method,  we  generate  a  Benders  cut  to  solve  the  relaxed  master  problem  which  is 
equivalent  to  original  one  in  terms  of  yhj' s.  After  we  get  the  new  set  of  yhj' s,  we  solve  the 
new  MIP  with  new  y^,  and  we  form  a  new  Benders  cut  and  solve  the  relaxed  master 
problem  for  the  new  set  of  y*Y  s  until  we  come  to  a  point  that  we  cannot  improve  the 
relaxed  master  problem  objective  function.  At  that  point  we  assumed  that  we  reached  the 
optimal  solution.  I  attached  the  first  iteration  of  the  MIP  with  fixed  y°i=  y22=  y*3=l  at 
Appendix  8.  We  got  the  same  result  as  the  one  with  Federgrun  and  Zipkin  procedure  in 
the  first  iteration  as  we  expected.  Since  the  result  of  the  first  iteration  is  the  same  as  the 
previous  method,  it  means  we  will  get  the  same  Benders  cut  for  the  following  iterations. 
Therefore;  we  just  show  the  first  the  result  of  MIP  at  Appendix  9. 

In  the  first  iteration  we  have  no  control  on  the  solution  due  to  the  fact  that  we 
selected  the  y* jj’s  arbitrarily.  But  in  the  second  iteration,  we  come  up  with  the  y’Yj’s  such 
that  the  resulting  TSP  model  yields  the  shortest  tour  with  an  objective  function  14.5,  and 
associated  total  cost  including  inventory/allocation  objective  value  is  3489.5.  In  the  third 
iteration,  we  assigned  the  vehicles  to  the  nodes  with  a  little  bit  slower  assignment  like 
15.3  but  the  total  cost  goes  down  to  3440.3.  This  result  shows  us  that  combining  vehicle 
routing  and  inventory  yields  better  (lower)  result  than  treating  TSP  and 
inventory/allocation  individually. 
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Why  do  we  Need  an  Alternative  Approach  to  the  Starting  Point  Problem 


We  used  the  Generalized  Benders  Decomposition  technique  to  take  an  advantage 
of  special  structure  of  the  problem  and  to  solve  it  effectively  in  terms  of  computational 
time.  But  we  know  that  we  can  reformulate  the  same  problem  as  if  we  do  not  use  any 
decomposition  technique. 

Firstly,  since  we  want  to  extend  the  initial  problem  in  several  aspects,  we  want  to 
formulate  the  problem  without  paying  attention  to  the  solution  technique.  If  we  omit  the 
dummy  vehicle  from  the  original  model  formulation,  and  we  let  the  model  fix  the 
appropriate  yhi  ‘s  in  the  equation  (20)  without  dummy  vehicles,  we  come  up  with  the 
same  optimal  solution.  We  delete  the  dummy  vehicle  to  decrease  the  number  of  decision 
variables  at  hand. 

The  second  reason  for  deleting  the  dummy  vehicle  out  of  thesis  problem  is  related 
to  nature  of  the  problem.  When  we  attempt  to  solve  the  starting  point  model  without 
fixing  yVs,  we  obtain  a  solution  other  than  optimal  solution.  In  other  words,  we  let  the 
model  pick  the  appropriate  y*1  fs.  It  yields  the  following  result:  xV  =x1io=l,  and 
x202=x223=  x230=1,  Zi=150,  z2=50,  z3=150  and  y\  =  y22=  y2^  =  1  with  a  total  cost  of 
35 15.3.  Since  the  starting  point  model  is  mixed  integer  non-linear  model,  even  Hyper 
Lingo  cannot  find  the  real  optimal  solution.  It  gives  one  of  the  local  results. 

After  we  obtained  the  same  results,  we  decided  to  extend  the  model  in  different 
ways.  We  do  focus  on  modeling. 

The  first  milestone  in  the  thesis  process  is  to  show  whether  it  is  possible  to 
combine  inventory/allocation,  VRP  and  repair  problems  in  a  unique  model.  Regardless 
of  the  solution  technique  like  Benders  decomposition  we  focus  on  the  formulating  the 
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problem.  The  second  milestone  is  the  derivation  of  generalized  Benders  decomposition 
equations  for  thesis  problem. 


A  General  Look  at  the  Possible  Extensions  from  the  Starting  Model  -  New 
Decision  Variables 


After  we  determined  that  we  have  a  possibility  to  solve  the  starting  point  problem 
without  using  Bender’s  decomposition  method,  we  feel  comfortable  to  deal  with  the 
starting  point  problem  to  extend  in  several  directions.  One  of  the  possible  extensions  can 
be  about  focusing  on  the  newsboy  properties.  The  other  one  is  related  to  the  VRP  side  of 
the  toy  problem,  especially  focusing  on  the  delivery  part. 

We  should  define  the  extra  decision  variables  and  parameters  needed  for 
extension  of  the  problem. 


if  L.B<XyUp<U.B.  and/or 
otherwise 


L.B<Xydel<U.B 


Hj  is  a  continuous  variable  that  denotes  the  amount  of  allocation  among  available 
products  of  node  i  to  node  j. 

(28)  Zij  =  xhijup  +  xhijdel 

Zu  is  a  continuous  variable  that  denotes  the  amount  of  allocation  among  available 
products  of  node  i  to  supply  itself. 

(29)  Zu  =  xhjup  +  xhidel 

yhij  is  a  binary  variable  equal  to  one  if  vehicle  h  allocates  the  supply  from  node  i  to  node  j 
and  equal  to  zero  otherwise. 
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y,-  is  a  binary  variable  equal  to  one  if  node  i  is  such  a  sufficient  node  that  it  can  deliver 
outside  but  not  accept  repaired  items,  equal  to  zero  if  node  i  is  not  sufficient  itself  so  that 
it  can  accept  supplies  but  not  deliver  outside. 

A  is  the  amount  of  time  vehicle  type  h  spends  at  demand  point  I  (Chan,  1999  Draft:8-20). 
U(h)  is  the  time  that  a  vehicle  h  spends  “on  the  road”  (Chan,  1999  Draft:8-21). 

Here  are  some  examples  for  each  new  decision  variables: 

x{  is  the  vehicle  routing  decision  variable  as  described  earlier,  which  denote 

vehicle  h  goes  from  node  i  to  node  j  if  it  is  equal  to  one;  otherwise,  zero.  We  add  two 
component  to  the  original  vehicle  routing  decision  variable  such  that  these  are  related  to 
what  quantity  we  repaired  at  node  i,  and  what  quantity  we  hold  in  the  inventory,  xhijup 
and  xhjjdel,  respectively.  The  two  new  decision  variables  are  lying  within  the  upper 
bound  (UB)  and  lower  bound  (LB). 

We  give  two  examples  for  each  to  understand  what  those  are  such  that  x^up 
denotes  that  vehicle  1  pick  ups  the  repaired  items  at  node  2  and  delivers  to  node  3,  x202up 
denotes  that  vehicle  2  pick  ups  the  repaired  item  at  node  0  and  delivers  to  node  2.  x^up 
is  the  decision  variable  based  upon  the  repair  capabilities  of  node  2,  which  is  uniformly 
distributed  between  [0-9].  x202up  is  the  decision  variable  based  upon  the  repair 
capabilities  of  node  0,  which  is  uniformly  distributed  between  [0-15]. 

On  the  other  hand,  x1 23 del  denotes  that  vehicle  1  pick  ups  the  inventoried  items  at 
node  2  and  delivers  to  node  3,  and  x  02del  denotes  that  vehicle  2  pick  ups  the  inventoried 
items  at  node  0  and  delivers  to  node  2.  x^del  is  the  decision  variable  based  upon  the 
remaining  inventory  of  node  2  at  the  end  of  the  day,  and  x202del  is  the  decision  variable 
based  upon  the  remaining  inventory  of  node  0  at  the  end  of  the  day. 
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The  delivery  variable  in  the  toy  problem,  Zj,  has  two  purpose.  One  of  them  is  the 
amount  of  demand  at  node  i,  and  the  other  one  is  the  amount  of  delivery  from  main  depot 
to  node  i.  Actually,  we  have  one  index  of  z\  associated  with  one  main  depot  and  it  is 
enough  for  the  formulation  of  toy  problem.  When  we  assumed  that  every  node  can 
supply  the  other  node  and  itself  as  well,  then  we  add  the  second  index  to  the  delivery 
decision  variable,  zoi  is  the  delivery  amount  from  node  zero  to  node  one.  This  delivery 
amount  is  the  sum  of  repaired  and  inventoried  items  that  shipped  from  node  zero  to  node 
one  by  any  vehicle. 

Since  we  are  dealing  with  repaired  items  there  should  be  some  possibility  to  re¬ 
supply  itself.  We  defined  a  new  decision  variable,  zu,  to  show  that  each  node  can  re¬ 
supply  itself,  zu  is  the  amount  that  node  1  re-supplied  itself,  and  Z22is  the  amount  that 
node  2  re-supplied  itself,  zu  is  the  sum  of  xlup  and  xldel.  xiup  denotes  that  the  repair 
capabilities  are  used  for  satisfying  the  local  failures  whereas  xidel  denotes  that  we  satisfy 
the  failures  by  using  the  local  inventory. 

We  can  examine  the  placing  order  decision  variable  yht  in  the  same  manner  as  we 
did  to  Zj.  We  add  the  third  index  to  placing  order  constraints  yhjj  due  to  the  fact  that  each 
node  in  the  model  is  treated  like  a  depot.  yV  is  equal  to  one  if  vehicle  1  allocates  the 
supply  from  node  zero  to  node  one;  otherwise  zero.  In  other  words,  if  y'o^l,  then  z0i 
can  take  value,  otherwise  it  is  zero. 

The  last  new  decision  variable,  y„  is  a  binary  variable  that  allows  us  to  satisfy  the 
assumption  that  either  we  can  send  an  item  to  the  other  nodes  or  we  can  accept  items. 
Since  we  are  dealing  with  repaired  items,  it  would  make  no  sense  during  the  period,  say  a 
day,  send  three  items  to  outside,  and  accept  one  item.  We  add  this  variable  to  prevent 
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from  happening  this  situation.  If  yi  is  equal  to  one,  then  it  means  that  node  1  will  be 
sufficient  enough  to  re-supply  itself  and  can  supply  the  other  nodes.  If  yi  is  equal  to  zero, 
then  it  means  that  node  will  use  up  all  of  its  available  resources,  which  are  the  amount  of 
available  inventory  at  hand,  and  available  repair  capabilities  and  accept  some  additional 
items  from  other  nodes  as  well. 

Thesis  Problem 

Since  we  are  using  the  same  distance  matrix  as  we  used  in  the  toy  problem,  we  do 
not  need  any  update  in  the  thesis  problem.  In  addition  to  that  we  do  not  extend  the  toy 
problem  in  terms  of  VRP  perspective  in  depth  except  we  add  crew  duty  time  constraints 
including  some  constant  loading  and  unloading  time.  In  thesis  problem,  we  keep  the  other 
assumptions  of  the  VRP  at  the  same  level. 

We  mainly  focused  on  the  inventory/allocation  perspective  of  the  toy  problem. 
Notice  that,  the  main  depot  is  the  only  resource  to  deliver  to  the  other  nodes  in  the  toy 
problem  whereas  we  set  up  the  thesis  problem  in  such  a  way  that  we  can  deliver  from  any 
nodes  to  any  nodes  with  respect  to  the  inventory  cost  function. 
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Table  5:  The  Inventory/Allocation  Table  for  Thesis  Problem 


Node 

0 

1 

2 

3 

Max.  Inventory  (Dl) 

10 

4 

5 

5 

Min.  Inventory 

0 

0 

0 

0 

Initial  Inventory  (p,) 

2 

1 

2 

1 

Demand  (PDF) 

0.1 

0.25 

0.4 

Shortage  cost  (c* ) 

$5 

$20 

$25 

$15 

Surplus  cost  (Q) 

$10 

$5 

$5 

%5 

Newsboy  cost  function 

18  +  0  .5zio  + 

48.75  -  3.75zoi  - 

51  -  4z02  -  4zi2  - 

39  -  2z03  - 

0.5z2o  +  0.5z30 

3.75z2i  -  3.75z3i 

4z32 

2zi3  -  2z23 

Failure 

- 

[0-9] 

[0-10] 

[0-9] 

Repair 

[0-15] 

[0-6] 

[0-9] 

[0-9] 

The  objective  function 

(30)  Minimize  XX  X^,>^  +  X<7,  (z, ) 

ie/  jelheH  i 

The  objective  function  of  thesis  problem,  equation  (29)  is  the  same  as  the 
objective  function  of  starting  problem,  which  is  equation  (4).  It  is  to  minimize  the 
inventory  cost  and  traveling  cost.  It  is  noteworthy  to  explain  the  effect  of  our  new 
assumption;  that  each  node  can  supply  the  other  nodes  on  the  inventory  cost  function. 

JPi+zij 

C,(i-p,-z,l)dF,(0+  iel-j,  Vj 

Pi+Zy  0 
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For  example;  the  inventory  cost  function  at  node  1  is 
(32)  qi(zj)  =  48.75  -  3.75z0i  -  3.75z2i  -  3.75z3i 

We  will  show  how  we  get  this  inventory  newsboy  cost  function  for  node  one  step  by  step. 
Notice  that  the  demand  coefficients  of  node  1  are  the  same  regardless  to  where  it  was 
supplied.  You  can  see  Table  4  for  the  information  of  the  other  nodes.  As  you  recall,  we 
can  laterally  supply  from  any  node  to  any  node  in  our  model.  It  is  clear  that  we  use 
uniformly  distributed  demand  for  each  node. 


D\ 

c 

1  , 

A+Z01 

r 

1  , 

65 

15 

[ 

cl  uG  + 

f 

Cl — ds  — 

zm 

J 

A+Zoi 

Dl 

J 

0 

Dl 

4 

4  01 

Dl 

r 

1  , 

A +221 

c 

1  , 

65 

15 

f 

cl  ds  + 

[ 

Cl — ds  — 

j 

A  +  *21 

Dl 

J 

0 

Dl 

4 

4  21 

Dl 

r 

1  , 

A +^31 

r 

1  , 

65 

15 

| 

cl - ds  + 

| 

Cl — ds  — 

j 

a+z3i 

Dl 

J 

0 

Dl 

4 

4  31 

Node  one  can  demand  the  same  amount  from  every  other  node.  When  we  sum  up  three 
demand  functions,  we  come  up  with  equation  (32).  We  can  derive  the  inventory  cost 
functions  for  the  remaining  nodes  in  the  model  as  shown  below. 
qi(zi)  =  48.75  -  3.75zoi  -  3.75z2i  -  3.75z3i 
q2(z2)  =  51-  4z02  -  4zi2  -  4z32 
q3(z3)  =  39  -  2z03  -  2zb  -  2z23 
q0(z0)  =  18  +  0.5zio  +  0.5z20  +  0.5z30 

The  demand  function  of  node  zero  is  positive  because  its  shortage  cost  is  less  than  its 
surplus  cost.  It  is  not  a  non-decreasing  function. 
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The  constraints  are  as  follows: 

Equation  (7),  equation  (8),  equation  (9)  and  equation  (10)  are  unchanged.  We 
added  complete  set  of  sub-tour  breaking  constraints,  equation  (10),  to  the  thesis  problem 
just  to  be  safe. 

When  we  consider  allocation/inventory  constraints  we  put  a  lot  of  effort  to  modify 
into  the  new  settings. 

(33)  'ZZzyy-j  <V(h)  h  e  H 

iel  jel 

i*j 

(34)  'Zzy  <pt  Vi  €  / 

j*I 

(35)  ^jZy  <  Capacity J  V  j 

iel 

(36)  Exj-yjkO  Vi, h  e  H 

jeM  j 

(37)  VJ 

heH  iel 
i*i 

The  equivalent  of  equation  (12)  is  equation  (33)  in  the  thesis  problem.  It  ensures 
that  load  capacity  assigned  to  each  vehicle  h  is  within  its  capacity.  Note  that  we  can 
assign  load  to  vehicle  h  from  node  i  to  node  j  if  yhy  is  equal  to  one;  otherwise  we  can  not. 

The  equivalent  of  equation  (13)  is  equation  (34)  in  the  thesis  problem.  Since  we 
have  one  main  supplier  in  the  starting  point  problem,  we  had  one  constraint  for  main 
depot.  The  number  of  depots  in  the  thesis  model  will  be  as  many  as  the  number  of  nodes 
in  the  system.  Therefore,  we  would  have  four  supply  constraints  for  each  demand  point. 
The  number  of  units  available  plus  the  repair  capability  sum  up  to  the  supply  quantity  of 
each  node.  Node  one  has  one  unit  at  the  end  of  the  period,  and  it  can  repair  most  6  units. 
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Node  one  can  then  supply  at  most  seven  units  to  the  other  nodes,  including  re-supplying 
itself  as  well. 

There  is  a  demand  constraint  for  each  demand  point  besides  the  supply 
constraints.  Equation  (35)  ensures  that  demand  cannot  be  exceeded  at  each  demand 
point.  The  initial  level  of  node  one  is  100  units  and  the  maximum  inventory  level  of  node 
1  is  500  units  in  the  initial  model.  This  means  node  one  demand  is  at  most  400  units.  In 
the  thesis  model,  the  computation  of  demand  for  each  demand  point  will  be  different  for 
sure  due  to  the  fact  that  each  node  can  re-supply  itself.  If  the  number  of  failures  at  each 
node  minus  available  local  resources,  which  the  sum  of  inventoried  items  and  repaired 
items  is  greater  than  zero,  then  we  can  send  such  a  quantity  that  we  satisfy  the  failures 
and  maximize  the  inventory  level,  otherwise  we  can  only  maximize  the  inventory  level 
associated  with  that  node.  For  example,  failures  can  be  between  zero  and  nine  at  node 
one.  The  initial  level  and  maximum  inventory  level  of  node  one  is  one,  and  four, 
respectively.  The  repair  capability  is  up  to  six  units.  Failures,  which  is  nine  minus  one  at 
hand  plus  six  units  at  repair  equals  to  two.  We  have  two  units  shortages,  but  node  one 
can  demand  up  to  six  units  which  is  four  units  due  to  depot  capacity  plus  two  units 
shortages.  We  can  calculate  the  demands  of  other  nodes  in  the  same  fashion. 

It  is  time  to  talk  about  the  new  additional  constraints  to  satisfy  our  assumptions. 
(38)  XyUp  +  Xydel  - Mx~  <  0 

{x^p  +  Xjdel  if  i  =  j 

x-jup  +  Xydel  if  i  *  j 
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(40) 


Zy  +  My j  >  0  (i) 

Zy  ~ My }<0  (ii) 

Zy  +  My y  <  M  (iii) 

Zy  -  My  y  >  -M  (iv) 

z ,,  -  Local  Resource;  +  My .  >  0 

(41)  11  7 

Zjj  -Myj  >  -M 

(42)  xtup  +  XyUp  <  Rt 

(43)  Xjdel  +  Xydel  <  fii 

(44)  Y.A  Y.A + SI Vh 

iel  jel  iel  jel 

(45)  Zy-M'Zy-j^  0  Vi,  j  e  I 

i*j  heH 

where  big-M  is  the  largest  number  that  satisfies  the  related  constraints. 

Equation  (38)  assumes  that  we  can  deliver  if  the  vehicle  routing  arc  is  open.  We 
put  all  z’s  in  one  place  in  equation  (39)  so  that  we  can  understand  the  relation  between 
them. 

Equation  (40)  and  (41)  are  connected  with  each  other.  While  the  first  one 
satisfies  the  condition  that  we  either  deliver  the  items  to  the  other  nodes  or  get  the  items 
out  of  the  other  nodes,  the  second  one  ensures  that  if  we  get  the  items  out  of  other  nodes 
we  make  sure  that  that  node  has  used  up  all  its  available  resources.  I  would  like  to  give 
an  example  for  these  two  constraints.  Think  about  these  two  sets  of  constraints: 

(46)  zOl  +  z21  +  z31  +  50yl  =>  0  (i) 
zlO  +  zl2  +  zl3  -  50yl  =<  0  (ii) 
zOl  +  z21  +  z31  +  50yl  =<  50  (iii) 
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zlO  +  zl2  +  zl3  -  50yl  =>  -50  (iv) 

(47)  zll  -  7  +  50yl  => 0 
zll  -50yl  =>  -50 

where  big-M  is  50,  i  =  0,  2,  3  and  j  =  1,  and  7  in  equation  (47)  corresponds  to  possible 
local  resources,  the  sum  of  inventory  at  hand  and  repair  capability,  at  node  one.  If  yi  is 
equal  to  zero,  then  equation  (i)  and  (iii)  will  hold  which  means  we  can  get  the  items  out 
of  other  nodes,  equation  (ii)  and  (iv)  will  be  relaxed  which  means  we  cannot  deliver  to 
the  other  nodes.  If  yl  equal  to  one,  then  equation  (ii)  and  (iv)  holds,  which  means  we  can 
deliver  to  the  other  nodes,  and  equation  (i)  and  (iii)  will  be  relaxed,  which  means  we 
cannot  get  the  items  out  of  the  other  nodes.  Equation  (44)  satisfies  that  if  yi  is  equal  to 
zero,  it  forces  to  use  all  available  resource,  which  are  7  in  this  case. 

Equation  (42)  is  the  constraint  related  to  repair  capacity  in  the  node.  This  repair 
capacity  includes  what  it  repairs  for  itself  and  for  the  other  nodes. 

Equation  (43)  is  the  constraint  related  to  available  inventory  at  each  node.  This 
constraint  has  the  same  the  logic  as  the  repair  capacity.  It  includes  what  it  supplies  to 
itself  and  to  the  others  out  of  depot  inventory. 

We  modified  the  VRP  constraints  by  adding  equation  (42),  which  satisfies  the 
condition  that  we  cannot  exceed  available  crew-duty  hours  associated  with  each  crew 
team.  We  assumed  that  loading  and  unloading  takes  a  constant  time  such  as  half  an  hour 
in  the  model. 

We  accomplish  the  first  milestone  that  is  corresponding  to  formulating  the  thesis 
model.  In  the  second  milestone,  we  simplified  our  thesis  problem  so  that  we  can  apply 
generalized  Benders’  decomposition  technique. 
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The  Second  Milestone 


The  second  milestone  includes  the  application  of  generalized  Benders’ 
decomposition  technique  to  the  thesis  problem.  We  will  do  some  simplification  to  the 
thesis  problem  so  that  we  can  apply  the  technique.  We  will  keep  the  lateral  supply 
attribute  of  the  thesis  problem  in  the  second  milestone.  Every  node  can  supply  to  every 
other  node.  We  omit  the  some  constraints  out  of  thesis  problem  such  as  having  “accept  or 
deliver”  option  of  each  node,  and  tracking  the  inventoried  and  repaired  items 
individually. 

The  challenging  part  of  the  generalized  Benders’  decomposition  technique  is  to 
define  the  relaxed  master  problem,  which  is  equivalent  to  the  original  problem.  The 
Benders’  cut  for  the  starting  point  model  is  as  follows: 

(48)  z > QP+  Z(/>X  +  xh  +  *•*)+  Z  y°  +  Z  Z(*f  +  <P-) yf • 

heH  iel-O  heHiel-  0 

The  first  two  summations  are  referring  to  the  intercept  of  the  Benders’  cut,  and  the  last 
two  summations  are  representing  the  slope  of  the  Benders  cut.  Since  we  have  four  depots 
in  the  second  milestone,  we  should  write  down  four  additional  cuts.  The  only  difference 
from  the  starting  model  is  adding  more  indexes  to  the  cut. 

The  first  cut  is  generated  based  upon  main  depot  zero,  which  is  shown  in  equation 

(49) .  The  other  cuts  are  generated  based  upon  node  one,  two  and  three  in  equation  (50), 
(51),  and  (52),  respectively. 


(49) 


z  >  Q0P0  ■*" 


Eg 

heH 


)hVh  +  Kh 


+  K0  )  +  X 


0  0  , 
<Pi  y<  + 


/e/-0 


EEotf +  *>,*)  y,h 


heHiel-Q 
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(5°)  z>n,p,  +  Z(p\  +  k"  +  *•„*)+ E  v°y°  +  Z  Z(*f  +  «>‘)yf  ■ 

heH  /e/-l  fce//ie/-l 

(51)  z  > q2p2  +  Z(phvh  +  /ch  +  ^)+yZ  9i  y,°  +  E  Z(*i*  +  )y.h • 

heH  ieI-2  heHieI-2 

(52)  z  >n3p,  +  2(/<>V,  +  K*  +  *■„*)+ Z  wV  +  Z  Z(«f  +  n*)yf  • 

heH  ie/-3  hsHiel-  3 

We  want  to  explain  the  inventory  cost  portion  of  the  Benders’  cut.  In  equation  (49),  the 
slope  portion  of  the  equation  is  based  upon  the  operating  (Khj),  and  inventory  cost  (<phi). 
We  explain  the  operating  cost  in  the  generalized  Benders’  section  of  Chapter  Three 
explicitly.  We  modified  the  inventory  cost  function  to  the  modified  thesis  problem.  We 
can  define  the  inventory  cost  functions  associated  with  the  first  section  of  the  equation 
(49)  as  follows. 

<Pi  =  48.75 - 3.75  x  (z0,  +  z21  +  z31)  -  (D0  -  p0)z0l  -  (Q0  -  p0)z2l  -  (Ci0  -  p0)z3 , 

(53)  (p2  —  51  —  4  X  (z02  +  Zj2  +  Z32  )  —  (f^o  _  /^o)^02  _  _  /^o)^12  —  (^0  _  )^32 

(Pj,  =  39  —  2  X  (Zq3  +  Zj3  +  Z23  )  —  (f^o  —  Po  )^03  —  (^0  _  Po  )^13  —  (^0  _  Po)Z  23 
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The  second  portion  of  slope  equation  can  be  defined  as  follows. 


(px  =  48.75  —  3.75  x  (z01  +  z2\  +  z31 ) —  (^o  ”  Ao)^oi  “  (^o  ”  Ao)^2i  —  (^o  —  Ao)^ 31 
(p2  =  51  —  4  x  (Zq2  +  Z\2  +  Z32 )  “  (f20  —  Pq  )  Zq2  (f^o  Ao)^12  “  (^0  “  Ao  )^32 

^3  =  39  —  2x  (z03  +  +  Z23)  —  (fi0  “  Ao)^03  “  (fi0  ~~  Ao)^13  “  (^0  ”  Ao)^23 

W  =  48.75  -  3.75  x  (z01  +  z21  +  z31 )  -  (O0  -  A )z01  -  (^o  “  A  )*2i  “  (Qo  “  A )*3i 

(54)  (p2  =  5  1  —  4  X  ( Zq2  +  Z12  +  Z32  )  “  (^0  _  A  )^02  “  (&0  ~  Pi  )^12  ~  (Q0  _  Al  )^32 

(p3  =  39  —  2  X  (Z03  +  Z13  +  Z23  )  ”  (^0  “  Pi  )^03  “  (^0  “  Al  )^13  ”  (^0  “  Al  )^23 

V\  =  48.75  —  3.75  x  (z01  +  z2\  +  z3i )  “  (^o  A2 )^oi  ”  (^0  —  Pi  )^2i  “  (^0  Pi  )^3i 
<Pi  =  5 1  ■ -  4  x  (z02  +  Z12  +  z32)~  (do  ~  Pi )^02  -  (Q0  -  Pi )zn  “  (^0  "  Pi )z32 

(p3  —  39  —  2  X  (z03  +  Z13  +  Z23  )  “  (f^o  ~  Pi)Zq3  ~  (fi0  “  p2  )Z13  “  ""  Pi  )^23 

We  can  generate  Benders’  cut  for  the  node  one  in  the  same  fashion  as  follows. 

The  first  portion  of  slope  in  equation  (50)  is 

(pQ  =  18  +  0.5  X  (Z1Q  +  Z20  +  Z30)  —  (^1  ”  Ao  )^10  ”  (^1  “  Ao)^20  “  (^1  Ao  )^30 

(55)  (p2  =  5 1  —  4  x  (z02  +  Z12  +  Z32 )  “  (^1  —  Ao  )^02  —  (^1  “  Ao  )^i2  ~  (^1  —  Ao  )^32 

^ V 3  “  39  —  2  X  (z03  +  Zj3  +  Z23  )  “  (fij  —  Po  )Z03  ”  (^1  ”  Ao  )^13  “  (^1  “  Ao  )^*23 

and  the  second  portion  of  slope  in  equation  (50)  is 

(p§  —18  +  0.5  x  (Z|Q  +  Z2o  ^30  )  _  (^1  _  Ao  )^10  ”  (^1  —  Ao  )^20  “  (^1  *"  Ao  )^30 
(p2  =  5 1  —  4  X  (z02  +  Z12  ^32  )  —  (^1  “  Ao  )^02  “  (^1  —  Ao  )^12  ”  (^1  —  Ao  )^32 

(p3  =  39  —  2  X  (z03  +  z13  +  Z23)  “  (^1  “  Ao)^03  “  (^1  —  Ao)^13  —  (^1  —  Ao)^23 

(p\  =  18  +  0.5  x  (z10  +  z20  +  z30)  -  (d,  -  px  )z10  “  (Qi  -  A  )^2o  -  (Qi  -  A  )^so 

(56)  (p\  =  5 ^ 1  - 4  x  (z02  +  z12  +  z32 )  -  (Qj  -  a  K2  “  (4  “  A  Ki  "  (Qi  “  A  )^32 
^3  =  39  —  2  x  (z03  +  z13  +  z23 )  ~  (fi]  ”  P!  )z03  —  (f^i  “  Ai  )^i3  “  (^1  “  Ai  )^23 

\ Vo  =  18  +  0.5  x  (z10  +  z2q  +  Z3q)  “  (^1  “  p2 ) ^10  ~  (^1  ~  Pi)z>2q  ~  (^1  “  Pi)Z3q 
(Pl  =  5 1  —  4  X  (z02  +  Zj2  +  z32  )  “*  (f^!  “  p2  )Zq2  ~  (f^i  “  Pi  )zn  “  (f^i  “  p2  )^32 
y V\ 3  =  39 “ 2 x  (z03  +  Zj3  +  z23 ) —  (f^j  —  Pi)Zq3  ~~  (f^i  “  Pi )Zj3  ~ (fij  —  p2 ) z 23 
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If  we  generate  a  Benders’  cut  based  upon  node  zero,  we  show  that  we  deliver  to  node 
one,  two  and  three  explicitly  in  equation  (54).  If  we  generate  a  Benders’  cut  based  upon 
node  one,  we  show  that  we  deliver  to  node  zero,  two  and  three  explicitly  in  equation  (56). 
We  can  generate  the  cuts  based  upon  the  other  nodes  as  well.  Notice  that  the  intercept  of 
each  cut  regardless  of  supply  node  is  the  same.  Since  we  ran  out  of  time,  we  could  not 
provide  the  numerical  results. 

Now  we  are  ready  to  talk  about  the  solution  strategies  we  come  up  with  to  solve 
the  thesis  model. 
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Chapter  4  -  Analysis 


Introduction 

We  showed  the  formulation  of  thesis  problem  in  the  methodology  section.  We 
will  show  the  results  of  a  thesis  problem  without  using  Generalized  Benders 
decomposition  technique  that  corresponds  to  the  result  of  first  milestone  in  the  thesis 
process. 

First  Milestone 

We  would  like  to  give  you  some  information  about  our  original  model  that  we 
show  the  mathematical  formulation  of  it  at  Appendix  10.  We  have  a  total  of  124 
variables  and  117  constraints  in  the  model.  52  of  them  are  binary,  and  the  rest  of  them 
are  continues.  Binary  variables  are  placing  order  variables  (y^),  vehicle  routing 
variables  (xhy)  and  switch  on-off  variables  related  to  either  sending  or  receiving  items 
(yO.  We  have  only  two  non-linear  constraints  in  the  model  associated  with  vehicle 
capacity  constraints  that  correspond  to  the  90th  and  91st  constraints.  We  have  four  nodes, 
and  two  vehicles. 

If  we  do  not  fix  the  yh y’s,  the  branch-and-bound  can  procedure  can  take  252 
solution  trees.  Each  time  the  software  Lingo  can  follow  the  different  path  due  to  the 
different  storage  memory  produced  to  solve  the  given  problem  (Lingo  Technical 
assistance).  The  result  of  the  first  milestone  can  be  summarized  as  follows. 
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Table  6  :  The  result  of  the  first  milestone 


First  solution 

Second  solution  (Alternative) 

Objective  value 

Z  =  133.55 

Z*=  133.55 

Routes 

Y 1  ~I  T7I 

X  02  ~  X  23  -  X  30 

2  2 

X  01  -x  jo 

1  I  I 

X  03  -  X  32  -  X  20 

2  2 

X  01  -  X  io 

Initial  position 

- 

Node  0 

Initial  level: 2 

Max. level:  10 
Repair  Cap:  [0-1 5] 

Zoo  =  7 
(x0up=7) 

Zoo  =11 

( Xodel  =2,  x0up=9) 

Node  1 

Initial  level:  1 
Max.level:4 

Repair  Cap.: [0-6] 
Failure  :[0-9] 

Zn  =  7 

( xidel  =1,  xiup=6) 

Zoi=6 

( x20iup  =  4,  x2oidel=2) 

Z„  =  7 

( xidel  =1,  xiup=6) 

Zoi=6 

( x2oiup  =  6) 

Node  2 

Initial  level:2 
Max.level:5 

Repair  Cap:  [0-9] 
Failure  :[0-10] 

Z22=ll 

( X2del  =2,  X2up=9) 

Zo2=4 

(x102up=4) 

Z22=ll 

( x2del  =2,  x2up=9) 

Z32  =4 

(x132up=3,  x!32del=  1) 

Node  3 

Initial  level:  1 
Max.level:5 

Repair  Cap:  [0-9] 
Failure:  [0-9] 

Z33  =  10 

( x3del  =0.99999,  x3up=9) 

yzoi  =  y‘o2=y1o3=i 
y3  =  yo  =  l 

y2oi  =  y132=y103=i 
y3  =  yo  =  l 

We  can  go  backward  to  analyze  the  situation  such  that  we  send  six  serviceable 
parts  to  node  1.  There  can  be  at  most  nine  failures  at  node  one,  but  node  one  can  meet 
seven  of  them  by  using  our  local  resource.  It  means  node  1  use  all  of  his  available  repair 
and  inventory  to  meet  the  failures.  We  are  a  shortage  of  two  serviceable  parts.  Since  the 
maximum  inventory  capacity  is  four,  we  can  demand  at  most  six  serviceable  parts  to 
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stock  up  to  the  maximum  level.  Node  one  got  four  serviceable  items  out  the  repaired 
shop  and  two  serviceable  items  out  of  the  inventory  out  of  node  zero. 

For  node  two,  we  can  have  at  most  ten  failures  during  that  period.  But  we  can 
satisfy  eleven  of  them  by  using  the  local  resource.  Node  two  uses  all  of  its  repair 
capabilities  and  inventory  to  satisfy  the  failures.  Since  our  maximum  inventory  level  is 
four,  we  can  demand  at  most  four  serviceable  items  from  any  node  in  the  system.  Node 
two  got  four  units  out  of  the  repaired  shop  of  node  zero. 

For  node  three,  we  have  nine  failures  but  it  meets  all  the  failures  by  using  its  own 
local  resource,  besides  that  it  has  a  surplus  of  one.  Since  y3  is  equal  to  one,  it  is  ready  to 
send  it  to  the  other  nodes. 

The  placing  order  variables  are  consistent  with  the  result.  Since  we  are  sending 
out  of  node  zero,  5^01  =  y'o2=  y’o3  are  equals  to  one.  Besides  that,  we  deliver  to  node  one 
by  vehicle  two,  therefore,  placing  an  order  constraint  y2oi  and  associated  vehicle  routing 
variables  such  as  x2oi  -  x2iois  equal  to  one.  y^is  equal  to  one,  because  we  are  using 
vehicle  one  goes  out  of  node  two  that  is  consistent  with  the  constraint  set.  We  do  send 
nothing  from  node  three  to  any  other  node  even  though  y3  is  equal  to  one,  because  we 
satisfied  the  demand  requirements  of  node  one  and  two.  We  cannot  deliver  to  node  zero, 
because  it  supplies  to  node  one.  On  the  other  hand,  the  switch  on/off  variables  yi,  and  y2 
are  equal  to  zero  mean  that  node  one  and  node  two  can  only  demand  some  send 
serviceable  items  from  the  other  nodes.  In  this  result,  y3  did  send  nothing  because  node 
one  and  two  are  satisfied  by  node  zero.  If  node  zero  did  not  have  enough  to  send  them, 
then  node  three  would  be  ready  to  satisfy  their  demands.  Node  three  has  only  one  surplus 
item  to  deliver  to  the  other  nodes. 
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Vehicle  two  goes  to  node  one  and  come  back.  The  total  time  spent  of  vehicle  two 
is  the  sum  of  traveling  time  from  node  zero  to  node  1,  plus  0.5  hour  for  unloading  and 
plus  traveling  time  from  node  one  to  node  zero,  which  is  equal  to  six  and  half  hour.  The 
that  of  vehicle  one  is  sum  of  traveling  time  from  node  zero  to  node  two,  plus  unloading 
time  at  node  two,  plus  traveling  time  from  node  two  to  three,  and  traveling  time  from 
node  three  to  node  zero,  which  is  equal  to  nine  hours  and  eight  minutes.  Note  that 
vehicle  two  delivers  nothing  to  node  three  because  node  three  is  a  candidate  to  deliver  to 
the  other  nodes  due  to  its  low  shortage  cost  compare  to  node  one  and  two. 

We  have  an  alternative  solution  to  this  problem.  We  have  the  same  objective 
function  with  a  value  of  133.55,  but  we  satisfy  the  demands  of  node  one  and  two  from 
node  zero  and  node  three. 

For  node  one,  we  are  sending  six  repaired  items  out  of  node  zero  to  node  one  to 
satisfy  the  demand.  This  case  node  three  delivers  three  repaired  and  one  inventoried  item 
to  node  two  to  satisfy  the  demand  of  node  two.  Node  one  and  two  have  the  same  number 
of  failures  as  before  such  that  they  got  the  same  number  of  serviceable  items.  But  there  is 
a  big  difference.  At  this  time,  node  three  can  be  able  to  send  six  more  units  to  the  node 
two  if  there  is  a  need. 

The  switch  on/off  constraints  for  delivery,  yo  and  y3,  are  consistent  with  the 
results.  Remember  yo  is  equal  to  one  means  that  we  can  deliver  out  of  node  zero.  y3  is 
equal  to  one  means  that  we  can  deliver  out  of  node  three. 

The  placing  order  constraints,  y2ory! 32^03  are  consistent  with  the  result.  Since 
vehicle  two  goes  out  of  node  two,  y  01  arc  is  open.  By  the  same  idea,  vehicle  one  visits 
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node  three  and  two.  The  other  placing  order  constraints  are  consistent  with  the  result,  as 
well.  See  the  complete  list  of  non-zero  decision  variables  at  Appendix  1 1 . 

We  may  interpret  these  two  results  as  follows.  In  the  first  run,  we  satisfy  the 
demand  out  of  node  zero.  It  means  node  zero  it  is  the  supply  source  for  the  model  in 
order  to  get  the  optimal  result.  In  the  second  run,  the  supply  source  is  the  node  three  and 
node  zero.  In  this  problem  instance  we  have  two  options  to  satisfy  the  demands.  We 
may  think  of  it  in  terms  of  changing  the  place  of  supply  sources.  This  happens  usually  in 
the  deployment  in  military.. 

With  the  current  newsboy  inventory  costs,  we  may  think  that  we  may  deploy  all 
the  aircraft  from  node  three  to  node  one  and/or  two  and  treat  the  node  three  as  secondary 
main  depot. 

What  would  happen  if  we  change  the  demand  quantities  of  nodes  one,  two,  or 
three.  We  want  to  test  the  model  whether  it  has  a  debug  or  not.  The  new  set  of  demand 
quantities  is  such  that  node  one  requires  four  serviceable  items  in  stead  of  six,  node  two 
requires  the  same  amount,  that  is  four  serviceable  items,  and  node  three  requires  six 
serviceable  items  in  stead  of  four.  The  demand  constraints  for  nodes  one,  two,  and  three 
are  the  87th,  88th,  and  89th  ones  at  Appendix  10. 

We  displayed  the  results  in  Table  2.  We  would  like  to  talk  about  only  the  first 
run.  Since  we  changed  the  demand  quantity  of  node  one  from  six  to  four,  we  can  supply 
up  to  four  units.  Node  one  takes  only  a  total  of  four  units  out  of  node  zero.  These  are 
3.07  repaired  items  and  0.92  inventoried  item.  Note  that  we  defined  these  variables  as 
continuous.  Node  two  takes  a  total  of  four  serviceable  items  out  of  node  three  such  that 
one  of  them  is  from  inventory  and  the  rest  is  from  repaired  shop. 
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Table  7:  The  Demands  of  Node  one  and  three  Modified 


First  solution 

Second  solution  (Alternative) 

Objective  value 

Z*=  141.05 

Z*=  141.05 

Routes 

—~2  2  2  ” 

A  03  -  X  32  ~  X  20 
v1  v1 

X  01  -  X  io 

rr  ~r 

X  02  -  X  23  -  X  30 

2  2 

X  01  -  X  10 

Initial  position 

- 

Node  0 

Initial  level: 2 

Max. level:  10 
Repair  Cap:[0-15] 

Zoo  =13 

(xoup=11.93,  xodel=1.07 ) 

Zoo  =  9 

( x0del  =2,  x0up=7) 

Node  1 

Initial  level:  1 
Max.level:4 

Repair  Cap.: [0-6] 
Failure :  [0-9] 

Zu=7 

( xidel  =1,  X]Up=6) 

Zo,=4 

( xVup  =  3.07,  xVde^O^) 

Zn  =  7 

( xidel  =1,  xiup=6) 

Zoi=4 

( x2mup  =  6) 

Node  2 

Initial  level: 2 
Max.level:5 

Repair  Cap:  [0-9] 
Failure  :[0-10] 

Z22=ll 

( X2del  =2,  x2up=9) 

Z32  =4 

(x202up=3,  x202del=l ) 

Z22=ll 

( x2del  =2,  x2up=9) 

Zo2=4 

(x132up=3,  x132del=  1) 

Node  3 

Initial  level:  1 
Max.level:5 

Repair  Cap:  [0-9] 
Failure:  [0-9] 

Z33  =  10 

( X3del  =0.99999,  X3Up=9) 

y'oi  =  y232 
y3  =  yo  =  l 

y2oi  =  yJ02 
y3  =  y0  =  1 

Note  that  yV,  y232  are  consistent  with  the  routes  and  delivery  variables.  Vehicle 
one  visits  and  delivers  to  node  one,  vehicle  two  visits  and  delivers  to  node  two.  y3,  yo  are 
consistent  with  the  results  as  well.  We  deliver  out  of  node  three  and  node  zero. 

We  generate  three  solutions  to  illustrate  the  importance  of  the  combining  the  inventory 
and  vehicle  routing  problems. 
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Table  8  :  The  result  of  the  thesis  problem 


options 

Zi 

Operating  cost 

Inventory  cost 

Total  cost 

1 

Zo3  =  4 
Z2,=6 

14.5 

126.25 

140.75 

2 

Zoi  =  6 

Zo2  =  4 

17 

118.25 

135.25 

3 

Zoi  =  6 

Z32  =  4 

15.3 

118.25 

133.55 

In  the  first  iteration,  we  follow  the  shortest  route  to  satisfy  the  demands.  If  we 
want  to  satisfy  the  demands  in  the  quickest  way,  we  will  apply  the  first  option.  Notice 
that  the  total  cost  is  the  highest  among  the  possible  options  for  this  numerical  example. 
The  second  and  third  option  has  the  same  inventory  cost  function,  but  different  operating 
costs.  Notice  that  in  the  second  and  third  option,  the  demand  points  get  what  they  need. 
Node  two  gets  4  units  out  of  node  zero  in  the  first  option,  or  out  of  node  three  in  the 
second  option.  The  third  option  results  in  a  lowest  cost.  Even  though  we  are  satisfying 
the  customer  demands  a  little  bit  slower  than  we  do  in  the  first  option,  we  have  the 
optimal  results.  This  table  shows  us  that  we  will  have  a  lower  cost  if  we  consider 
inventory/allocation  and  vehicle  routing  together.  Even  though  we  extend  Federgruen 
and  Zipkin  in  the  thesis  model,  we  come  up  with  the  same  conclusion  as  they  found  in 
their  study. 


84 


The  impact  of  changing  the  number  of  nodes  and  vehicles  on  the  problem 


The  impact  of  adding  one  more  node  into  model  can  be  summarized  as  follows: 

It  is  for  sure  that  we  expand  the  number  of  constraints  while  adding  the  number  of  nodes 
into  the  model.  We  try  to  mention  the  impact  of  it  on  the  each  type  of  constraint.  We 
add  one  more  constraint  for  out/in  of  main  depot  constraints,  which  are  equation  (7)  and 
(8),  respectively.  We  add  one  for  each  vehicle  type  for  the  route  continuity  constraint, 
which  corresponds  to  equation  (9).  Since  we  are  not  changing  the  number  of  vehicles  in 
the  model,  we  will  have  the  same  number  of  constraints  in  terms  of  tracking  the  vehicles 
on  the  way  to  main  depot,  which  corresponds,  to  equation  (10).  We  add  four  additional 
sub-tour  breaking  constraints  for  each  vehicle  type,  which  is  originally  equation  (11). 
Since  equation  (33)  is  related  to  vehicle  capacity  constraints,  it  will  have  the  same 
number  of  constraints.  We  add  one  more  constraint  for  equation  (34),  and  (35).  We  will 
add  one  placing  order  constraint  for  each  vehicle  type  for  equation  (36).  Since  the 
equation  (37)  is  based  on  the  number  of  nodes,  we  will  add  one  more  constraint.  The 
situation  for  equation  (38)  will  be  different.  We  add  one  for  each  previous  node,  and  five 
for  the  new  coming  node  and  two  for  the  equation  (39).  The  number  of  equations  for 
(40)  will  remain  the  same,  but  we  add  two  more  constraints  for  equation  (41),  (42),  and 
(43).  The  equation  (45)  will  be  the  same. 

When  we  think  of  multi-item  in  our  model,  we  add  one  more  index  to  the  delivery 
variables,  and  add  a  newsboy  cost  function  for  each  item.  The  new  delivery  variable  will 
be  Zijk,  denotes  item  k  goes  from  i  to  j.  The  common  resources  will  be  vehicle  routing 
constraints,  and  vehicle  capacity  constraints.  Every  item  can  be  analyzed  independently. 
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Every  possible  solution  can  pass  up  to  the  common  resources.  As  soon  as  we  satisfy  the 
common  resources,  we  assumed  that  we  reached  the  optimal  solution. 

We  may  think  of  adding  more  nodes  in  a  generic  way  as  well.  We  will  present 
the  equations  again  for  the  sake  of  completeness.  We  will  assume  that  we  add  n  nodes 
more,  which  can  be  denoted  like  II7 1. 


(7) 

(8) 


\H"\  if  i  =  0 
1  if  i  =  1,2,...,  |I| 

\H"\  if  j  =  0 
1  if  j  =  l,2,...,|l| 


Equation  (7)  and  equation  (8)  are  going  to  increase  by  II7 1.  In  other  words  the  new 
number  of  constraints  for  equation  (7)  and  (8)  will  be  III  +  II7I. 

(9)  Vh,  Vp  e  I. 

ieMp  jeMp 


Equation  (9)  is  going  to  increase  by  IH77I  times  II7I.  This  means  we  write  down  the  new  set 
of  constraints  for  each  vehicle  type. 

(10)  Z  xj  <  1  Vh,  and  <  i  Vh 

jeM0  ieM0 


The  number  of  constraints  in  equation  (10)  remains  the  same  due  to  the  fact  that  it 
depends  on  the  vehicles. 

(11)  ZZ*J<|L|-1  Lc{2,...,|/|},  2<|L|<|/|-1;  he  H 

ieL  jeL 

Constraint  (11)  is  going  to  increase  in  the  following  fashion.  Cr(1, 1-1)  denotes  the 
combination  of  “I  chooses  1-1”.  If  we  add  I7  nodes  to  the  problem,  then  the  number  of 
constraints  is  the  summation  of  Cr(I+I7, 1+I7),  Cr(I+I7, 1+I7-l),  Cr(I+I7, 1+I7-2),  up  to 
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Cr(I+I/,  1).  Remember  node  zero  is  the  main  depot.  In  addition  to  that  multiply  each 
combination  by  the  number  of  vehicles. 

(33)  ZZ zijyl<V(h)  h  e  H 

is/  jel 
•*j 

Equation  (33)  remains  the  same  due  to  the  fact  that  it  based  on  the  number  of  vehicle 
types. 

(34)  Z  Zy  <  P>  Vi  e  I 

jel 

Equation  (34)  is  going  to  increase  by  II7 1.  In  other  words  the  new  number  of  constraints 
for  equation  (34)  will  be  III  +  II7!. 

(35)  Z Zy  <  capacity .  Vj 

i  el 

Equation  (35)  is  going  to  increase  by  the  II7 1. 

(36)  Z*i-y!*0  Vi,h  e  H 

jeMj 

Equation  (36)  is  going  to  increase  by  IH7/I  times  II7I.  In  our  case,  if  we  add  two  more 
nodes,  the  additional  number  of  constraints  to  the  original  would  be  four. 

(37)  ZZyJsi  V / 

heHiel 

i*j 

Equation  (37)  is  going  to  increase  by  II7 1. 

(38)  x^up  +  x*jdel  -  Mx~  <  0 

Suppose  that  the  total  number  of  nodes  would  be  equal  to  II7  _newl  after  adding  II7 1  nodes 
to  the  original  problem.  Then,  we  have  the  following  number  of  constraints:  II7  _new- 
1I.II7  newl.lH77!. 
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{XjUp  +  xtdel  if  i  =  j 

xfjUp  +  Xydel  if  i  ^  j 

We  can  define  zy  for  two  different  situations.  If  i=j,  then  we  will  add  II7 1  additional 
nodes.  If  i*j,  then  II7  _new-ll.II7  _newl  constraints  are  added. 


(40) 


Zy  +  My  j  >  0  (i) 

Zy  -  My}  <  0  (ii) 
Zy  +  Myj  <  M  (iii) 
Zy  -  Myj  >  -M  (iv) 


Equation  (40)  is  composed  of  four  sub  constraints.  It  is  going  to  increase  by  4II7 1. 

z  „  -  Local  Resource:  +  My ..  >  0 

(41) 

^  ;  zjj-Myj>-  50 

Equation  (42)  is  composed  of  two  sub  constraints.  It  is  going  to  increase  by  2II7 1. 

(42)  xt  up  +  Xy  up  < 

(43)  xtdel  +  Xydel  <  J3t 

Equation  (42)  and  (43)  are  going  to  increase  by  II7 1. 

(44)  I >*  Z4  +  HUdyxfj  <  V(h)  Vh 

iel  jel  16/  jel 

Equation  (44)  is  based  upon  the  number  of  vehicles. 

(45)  zrM^<  0  Vi, j e  I 


HeH 


The  new  number  of  constraints  for  equation  (45)  would  be  II7  _new-ll.II7  _newl. 
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We  can  answer  the  impact  of  the  adding  more  vehicles  to  the  problem  as  well. 

IV7!  denotes  additional  vehicles  added  into  the  original  model. 

Equation  (9),  (10),  (33),  (36),  and  (44)  are  going  to  increase  by  IV7!.  Equation  (1 1) 
is  going  to  increase  by  the  multiplication  of  IV7I  with  every  combination. 
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Chapter  5 


Summary  of  This  Research 

The  thesis  process  started  with  the  toy  problem  that  we  studied  throughout  the 
process.  Its  size  was  small  enough  to  capture  the  initial  requirements  of  the  thesis 
problem.  It  was  complex  enough  to  be  a  good  candidate  for  thesis  problem. 

Firstly,  we  solved  and  understand  the  toy  problem  and  its  environment.  It  was 
composed  of  inventory  and  vehicle  routing  problems.  After  we  felt  comfortable  on  the 
toy  problem,  we  decided  to  head  for  the  formulating  the  thesis  problem. 

We  wanted  to  add  the  third  component,  which  is  repair  to  the  toy  problem.  It  was 
a  challenging  task  to  combine  all  three  components  in  a  unique  model.  The  first  question 
that  we  wanted  to  answer  if  it  is  doable  or  to  what  extend  it  is  doable. 

We  draw  the  boundary  of  the  problems  as  follows.  There  is  a  repair  capability  for 
each  base  other  than  the  main  depot  as  well.  There  is  a  newsboy  inventory-cost  function 
associated  with  each  base  to  the  every  other  base.  We  set  up  a  model  for  repairable  items 
only.  Each  base  has  two  options:  accept  the  repairable  items,  or  deliver  the  repairable 
items,  but  not  both.  If  the  base  chooses  the  “accept”  option,  we  make  sure  that  that  base 
has  used  all  its  available  resources,  initial  inventoried  items  and  repair  capabilities,  before 
it  receives  re-supplies.  In  the  transportation/delivery  submodel,  we  place  limitation  on 
the  crew  duty  hours  available. 

The  objective  function  is  to  satisfy  the  demands  by  minimizing  traveling  cost  and 
inventory  cost.  Satisfying  the  demands  is  more  dominant  than  minimizing  traveling  cost; 
therefore,  the  satisfying  the  demand  is  the  driving  factor  in  the  objective. 
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Conclusion 


The  purpose  of  this  model  can  be  summarized  in  three  steps.  In  the  first  step,  we 
analyze  and  comprehend  the  inventory/allocation  model  (Federgruen  and  Zipkin, 
1984:1019).  Here  we  take  an  advantage  of  a  numerical  example  in  Chan’s  textbook 
(Chan,  1999-Draft:  9-6).  In  the  second  step,  we  would  apply  the  generalized  Benders 
decomposition  technique  to  the  problem  as  Federgruen  and  Zipkin  did.  Benders  cut  is 
playing  a  major  role  in  the  solution  process.  Since  the  Benders  cut  is  composed  of  dual 
variables,  we  performed  a  sensitivity  analysis  for  the  dual  variables  associated  with  main 
depot  capacity  and  vehicle  capacity  constraints.  The  starting  point  model  shows  us  that 
combining  vehicle  routing  and  inventory  saves  a  lot  of  resource  in  terms  of  time  or 
money.  In  the  third  step,  we  extend  the  starting  point  model,  which  corresponds  to  our 
thesis  model,  in  such  a  way  that  we  let  the  model  issue  the  lateral  supply  between  every 
nodes  in  the  model.  We  track  the  items  individually;  namely,  they  may  come  from  either 
inventory  or  repaired  shop  of  the  node.  When  we  examined  the  result  of  our  thesis 
problem,  we  saw  that  we  would  have  some  savings  in  terms  of  time  and  money  if  we 
combined  the  inventory  and  vehicle  routing  together  as  Federgruen  and  Zipkin  did  in 
their  study.  Since  we  ran  out  of  time  we  could  not  show  the  result  of  the  generalized 
Benders’  decomposition  equations  for  the  second  milestone.  But  it  is  not  working 
properly. 

Suggestions  and  Recommendations 

We  came  up  with  the  mixed  integer  programming  since  we  used  the  uniformly 
distributed  demand  function. 
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Since  I  run  this  model  a  couple  of  times,  running  an  experiment  on  the  parameters 
of  the  thesis  problem  may  yield  interesting  results.  The  computational  time  can  be  a  huge 
hindrance  if  you  expand  the  nodes  in  the  model.  Therefore,  it  would  be  good  to  apply  the 
Generalized  Benders  decomposition  technique  to  thesis  problem  so  that  you  can  run  the 
model  as  much  as  you  want  in  a  short  computational  time. 

The  second  advantage  of  Generalized  Benders  decomposition  is  to  have  an 
opportunity  to  be  able  to  attack  Exponential  distributed  demands  at  inventory  problems 
and  Poisson  distributed  demands  at  repair  problems. 

This  extension  is  valid  for  Turkish  Air  Force  such  that  it  would  be  useful  to 
include  multi-item  in  stead  of  expanding  the  model  in  terms  of  number  of  nodes.  At  most 
20-25  items  will  be  enough  for  Turkish  Air  Force. 

Another  extension  can  be  the  use  of  heterogeneous  vehicles  in  the  model.  Since 
Turkish  Air  Force  has  three  different  kinds  of  vehicles  for  transportation  in  the  inventory, 
it  would  result  in  a  very  useful  model.  These  three  aircraft  types  are  stationed  at  two 
different  bases.  Formulating  the  problem  by  assuming  that  those  are  stationed  at  two 
different  bases  will  add  a  complexity  to  the  problem. 

The  loading  and  unloading  times  can  be  a  variable  associated  with  a  known 
distribution  instead  of  constant  values.  This  allows  us  to  focus  on  the  land  operations 
side  of  vehicle  routing  problem. 
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Appendix  1 :  The  TSP  Sub-problem  for  Starting  Model 


MIN  3.8  X032  +  5  X031  +  2  X030  +  3.8  X023  +  4  X021  +  3.5  X020  +  5  X013 
+  4  X012  +  3  X010  +  2  X003  +  3.5  X002  +  3  X001  +  3.8  X232 
+  5  X231  +  2  X230  +  3.8  X223  +  4  X221  +  3.5  X220  +  5  X213 
+  4  X212  +  3  X210  +  2  X203  +  3.5  X202  +  3  X201  +  3.8  X132 
+  5  X131  +  2  X130  +  3.8  X123  +  4  X121  +  3.5  X120  +  5X113 
+  4  XI 12  +  3  X110  +  2  X103  +  3.5  X102  +  3  X101 
SUBJECT  TO 

2)  X230  +  X220  +  X210  +  X130  +  X120  +  XI 10  =  2 

3)  X031  +X021  +X001  +X231  +X221  +X201  +X131  +X121  +X101  =  1 

4)  X032  +  X012  +  X002  +  X232  +  X212  +  X202  +  X132  +  XI 12  +  X102  =  1 

5)  X023  +  X013  +  X003  +  X223  +  X213  +  X203  +  X123  +  XI 13  +  X103  =  1 

6)  X203  +  X202  +  X201  +  X103  +  X102  +  X101  =  2 

7)  X013  +  X012  +  X010  +  X213  +X212  +  X210  +  XI 13  +  XI 12  +  XI 10  =  1 

8)  X023  +  X021  +  X020  +  X223  +  X221  +  X220  +  X123  +  X121  +  X120  =  1 

9)  X032  +  X031  +X030  +  X232  +  X231  +X230  +  X132  +  X131  +X130=  1 

10)  X131  +X121  -XI 13  -X112 -X110 +  X101  =  0 

11)  X231  +X221  -X213-X212-X210  +  X201  =  0 

12)  X031  +  X021  -  X013  -  X012  -  X010  +  X001  =  0 

13)  X132  -  X123  -  X121  -  X120  +  XI 12  +  X102  =  0 

14)  X232  -  X223  -  X221  -  X220  +  X212  +  X202  =  0 

15)  X032  -  X023  -  X02 1  -  X020  +  XO 12  +  X002  =  0 

16) - X132  -  X131  -X130  +  X123  +X113  +  X103  =  0 

17) -  X232  -  X23 1  -  X230  +  X223  +  X213  +  X203  =  0 

18) -  X032  -  X031  -  X030  +  X023  +  X013  +  X003  =  0 

19) -  Yll  +  X113  +  XI 12  +  XI 10  >=  0 

20) -  Y21  +  X213  +  X212  +  X210  >=  0 

21)  X013 +X012 +  X010>=  1 

22) -  Y12  +  X123  +  X121  +  X120  >=  0 

23)  X223  +  X221  +  X220  >=  1 

24) -  Y02  +  X023  +  X021  +  X020  >=  0 

25)  X132 +  X131  +X130  >=  1 

26) -  Y23  +  X232  +  X23 1  +  X230  >=  0 

27) -  Y03  +  X032  +  X03 1  +  X030  >=  0 

28)  X132  +  X131  +X123  +X121  +X113  +  X112<=  2 

29)  X232  +  X231  +  X223  +  X221  +  X213  +  X212  <=  2 

30)  X032  +  X03 1  +  X023  +  X021  +  X013  +  X012  <=  2 
END 

INTEY11 
INTEY21 
INTE  Y01 
INTEY12 
INTE  Y32 
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INTE  Y02 
INTE  Y13 
INTE  Y23 
INTE  Y03 


Appendix  2:  The  Solution  to  TSP  Sub-problem  for  the  Starting  Model 

Global  optimal  solution  found  at  step:  24 

Objective  value:  17.00000 

Branch  count:  0 


Variable 

Value 

Reduced  Cost 

X103 

1.000000 

0.0000000 

X130 

1.000000 

0.0000000 

X202 

1.000000 

0.0000000 

X220 

1.000000 

0.0000000 

X001 

1.000000 

0.0000000 

X010 

1.000000 

0.0000000 

Y01 

1.000000 

0.0000000 

Y22 

1.000000 

0.0000000 

Y13 

1.000000 

0.0000000 
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Appendix  3:  The  Inventory/Allocation  Sub-problem  for  Starting  Model 

MIN  -6.5Z3-7.5Z2-8Z1 
SUBJECT  TO 


2)  Z3+Z2  +  Z1  <=  400 

3)  Z1  <= 

400 

4)  Z2<= 

50 

5)  Z3  <= 

350 

6)  Z1  <= 

0 

7)  Z3<= 

150 

8)  Z2<= 
END 

200 
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Appendix  4:  The  Solution  to  the  The  Inventory/ Allocation  Sub-problem  for  Starting 
Model 

Global  optimal  solution  found  at  step:  3 

Objective  value:  -1350.000 

Variable  Value  Reduced  Cost 
Z2  50.00000  0.0000000 

Z3  150.0000  0.0000000 

Row  Slack  or  Surplus  Dual  Price 
4  0.0000000  7.500000 

6  0.0000000  8.000000 

7  0.0000000  6.500000 
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Appendix  5:  The  Calculation  of  Master  Cut  for  the  Master  Problem  (Starting) 


ITERATION  1 


Q  :=0  zl  :=0  z2 =  50  z3:=150 


'  8  ' 

6 

'  3  ' 

i 

o 

o 

CO 

zl 

■  0 

6.5 

k  := 

4 

KhO  := 

2 

Khi  :  =  | 

0  0  3.5 

z  := 

z2 

v  := 

150 

0 

7 

3.5 

0  2  0 

z3 

'200 

v  := 


2600-  8'zl-  (n  -t-  pj-zl 

437.5-  7.5-z2-  (n  -hp,)-z2 

3012.5-  6.5' z3-  (to  -f  pj-z3 


2600-  8zl-  (a  +p2)-zl 

437.5-  7.5-z2-  ([1  +  p2)-z2 

3012.5-  6.5- z3-  (fl  -HP2)-z3 
suml  :=-(v  -h  Khi) 


2600-  8  zl- (Q +p  Vzl 

vO  := 

437.5-  7.5-z2-  (ft  +  p  J-z2 

v0  = 

2.6*10? 

-337.5 

3012.5-  6.5  z3-  (Cl  +pj-z3 

837.5 

2600-  8zl-  -|-  p3^  *zl 

437.5-  7.5- z2-  (a  +p3)-z2 

3012.5-  6.5  z3-  (n  +p3)’z3 


2.6*10? 

2.6*10? 

2.6*10? 

-337.5 

-262.5 

62.5 

837.5 

1.063-103 

2.038*10? 

(-v°>i  j+  (sumlj  ,)  =  "5.203*  10? 


THE  MATRIXv  CORRESPONDS  TO  THE 
INVENTORY  COST  OCCURED  AT  EACH  DEMAND 
POINT.  EACH  ROW  CORRESPONDS  TO  THE 
NODE  1,  2,  and  3,  RESPECTIVELY. 

EACH  COLUMN  CORRESPONDS  TO  THE 
VEHICLE  0,  1,  and  2  RESPECTIVELY. 


THESE  ARE  THE  COEF.  OF  DUMMY  VEHICLE 
(-  v0)2  ,  -t-  suml2  j  -  675  PLACING  ORDER  VARIABLES  (yOl  -y02-y03) 

IN  THE  MASTER  CUT. 

(- v0)3  1  +  suml3  =  -1.675*10? 


(suml)  = 


-2.603*10? 

337.5 


-837.5 


-2.6*10?  -2.6*10? 

262.5  -66 

-1.065*10?  -2.038*10? 


THE  SECOND  AND  THIRD  COLUMN 
CORRESPOND  TO  THE  VEHICLE  1 ,  AND  2. 
FOR  EXAMPLE  SECOND  COLUMN,  FIRST 
ROW  IS  THE  COEF.  OF  yl  1 ,  THIRD  ROW  IS 
THE  COEF.  OF  y13  or  THE  SECOND 
COLUMN,  FIRST  ROW  IS  THE  COEF.OF  y21 
IN  THE  MASTER  PROBLEM. 


z 

h  =  1 


Kh+KhVPh'\ 


1000.5 


THIS  IS  THE  RIGHT-HAND  SIDE  OF  MASTER  CUT. 


98 


Appendix  6:  The  Formulation  of  Master  Problem  for  the  Starting  Model 


MIN  Z 
SUBJECT  TO 

2) -  5203  Y01  +  675  Y02  -  1675  Y03  -  2600  Y1 1  +  262.5  Y12  -  1065  Y13 

-  2600  Y21  -  66  Y22  -  2038  Y23  +  Z  >=  1000.5 

3)  50  Y02  +  150  Y03  <=  0 

4)  50  Y12  +  150  Y13  <=  150 

5)  50  Y22  +  150  Y23  <=  200 

6)  Y01  +  Yll +Y21  =  1 

7)  Y02  +  Y12  +  Y22  =  1 

8)  Y03  +  Y13  +  Y23  =  1 


END 

INTE  Y01 
INTE  Y02 
INTE  Y03 
INTE  Yll 
INTE  Y12 
INTE  Y13 
INTEY21 
INTE  Y22 
INTE  Y23 
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Appendix  7:  The  Solution  to  the  Master  Problem  for  the  Starting  Model 


Global  optimal  solution  found  at  step:  4 

Objective  value:  4731.500 

Branch  count:  0 


Variable 

Value 

Reduced  Cost 

Z 

4731.500 

0.0000000 

Y13 

1.000000 

1065.000 

Y21 

1.000000 

2600.000 

Y22 

1.000000 

66.00000 

100 


Appendix  8:  The  MIP  formulation  to  the  Starting  Model 


MIN  -  6.5  Z3  -  7.5  Z2  -  8  Z1  +  3.8  X032  +  5  X03 1  +  2  X030  +  3.8  X023 
+  4  X021  +  3.5  X020  +  5  X013  +  4  X012  +  3  X010  +  2  X003 
+  3.5  X002  +  3  X001  +  3.8  X232  +  5  X231  +  2  X230  +  3.8  X223 
+  4  X221  +  3.5  X220  +  5  X213  +  4  X212  +  3  X210  +  2  X203 
+  3.5  X202  +  3  X201  +  3.8  X132  +  5  X131  +  2  X130  +  3.8  X123 
+  4  X121  +  3.5  X120  +  5  XI 13  +  4  XI 12  +  3  X110  +  2  X103 
+  3.5X102  +  3X101 

SUBJECT  TO 

2)  X230  +  X220  +  X210  +  X130  +  X120  +  XI 10  =  2 

3)  X031  +X021  +X001  +X231  +X221  +X201  +X131  +X121 +X101  =  1 

4)  X032  +  X012  +  X002  +  X232  +  X212  +  X202  +  X132  +  XI 12  +  X102  =  1 

5)  X023  +  X013  +  X003  +  X223  +  X213  +  X203  +  X123  +  XI 13  +  X103  =  1 

6)  X203  +  X202  +  X201  +  X103  +  X102  +  X101  =  2 

7)  X013  +X012  +  X010  +  X213  +X212  +  X210  +  X113  +  X112  +  X110  =  1 

8)  X023  +  X021  +  X020  +  X223  +  X221  +  X220  +  X123  +  X121  +  X120  =  1 

9)  X032  +  X031  +X030  +  X232  +  X231  +X230 +  X132 +  X131  +X130=  1 

10)  X131  +X121  -  XI 13  -  XI 12  -  XI 10  +  X101  =  0 

1 1)  X231  +  X221  -  X213  -  X212  -  X210  +  X201  =  0 

12)  X031  +  X021  -  X013  -  X012  -  X010  +  X001  =  0 

13)  X132-X123-X121  -  X120  +  X112  +  X102  =  0 

14)  X232  -  X223  -  X221  -  X220  +  X212  +  X202  =  0 

15)  X032  -  X023  -  X021  -  X020  +  X012  +  X002  =  0 

16) - X132  -  X131  -X130  +  X123  +X113  +X103  =  0 

17) -  X232  -  X231  -  X230  +  X223  +  X213  +  X203  =  0 

18) -  X032  -  X031  -  X030  +  X023  +  X013  +  X003  =  0 

19)  XI 13  +  X112  +  X110  >=  0 

20)  X213  +  X212  +  X210  >=  0 

21)  X013  +  X012  +  X010  >=  1 

22)  X123  +  X121  +  X120  >=  0 

23)  X223  +  X221  +  X220  >=  1 

24)  X023  +  X021  +  X020  >=  0 

25)  X132  +  X131  +X130>=  1 

26)  X232  +  X23 1  +  X230  >=  0 

27)  X032  +  X03 1  +  X030  >=  0 

28)  X132  +  X131  +X123  +  X121  +X113  +  X112<=2 

29)  X232  +  X231  +  X223  +  X221  +  X213  +  X212  <=  2 

30)  X032  +  X031  +  X023  +  X021  +  X013  +  X012  <=  2 
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31)  Z3+Z2  +  Z1  <=400 

32)  Z1  <=  400 

33)  Z2  <=  50 

34)  Z3  <=  350 

35)  Z1  <=  0 

36)  Z3  <=  150 

37)  Z2  <=  200 
END 


102 


Appendix  9:  The  Solution  to  MIP  for  the  Starting  Model 


Global  optimal  solution  found  at  step:  15 

Objective  value:  4717.000 


Variable 

Value 

Reduced  Cost 

X103 

1.000000 

0.0000000 

X130 

1.000000 

0.0000000 

X202 

1.000000 

0.0000000 

X220 

1.000000 

0.0000000 

X001 

1.000000 

0.0000000 

X010 

1.000000 

0.0000000 

X 

1.000000 

0.0000000 

Z2 

50.00000 

0.0000000 

Z3 

150.0000 

0.0000000 

Y01 

1.000000 

0.0000000 

Y22 

1.000000 

0.0000000 

Y13 

1.000000 

0.0000000 

And  the  dual  variables  associated  with  vehicle  capacity  constraints  are  as  follows: 

35  0.0000000  8.000000 

36  0.0000000  6.500000 
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Appendix  10  :  The  formulation  of  First  milestone 


MIN  3  X101  +  3.5  X102  +  2  X103  +  3  XI 10  +  4  XI 12  +  5  XI 13  +  3.5  X120 
+  4  X121  +  3.8  X123  +  2  X130  +  5  X131  +  3.8  X132  +  3  X201 
+  3.5  X202  +  2  X203  +  3  X210  +  4  X212  +  5  X213  +  3.5  X220 
+  4  X221  +  3.8  X223  +  2  X230  +  5  X231  +  3.8  X232  -  4  Z32 
-  3.75  Z31  +  .5  Z30  -  2  Z23  -  3.75  Z21  +  .5  Z20  -  2  Z13  -  4  Z12 
+  .5  Z10  -  2  Z03  -  4  Z02  -  3.75  Z01 
SUBJECT  TO 

2)  X101  +  X102  +  X103  +  X201  +  X202  +  X203  =  2 

3)  XI 10  + XI 12  + XI 13  +  X210  +  X212  +  X213  =  1 

4)  X120  +  X121  +  X123  +  X220  +  X221  +  X223  =  1 

5)  X130  +  X131  +X132 +  X230  +  X231  +X232=  1 

6)  XI 10  +  X120  +  X130  +  X210  +  X220  +  X230  =  2 

7)  X101  +X121  +X131  +X201  +X221  +X231  =  1 

8)  X102  +  XI 12  +  X132  +  X202  +  X212  +  X232  =  1 

9)  X103  +  X113  +  X123  +  X203  +  X213  +  X223  =  1 

10)  X101  -  XI 10  -  XI 12  -  XI 13  +  X121  +X131  =  0 

11)  X201  -X210-X212-X213+X221  +X231  =  0 

12)  X102  +  XI 12  -  X120  -  X121  -  X123  +  X132  =  0 

13)  X202  +  X212  -  X220  -  X221  -  X223  +  X232  =  0 

14)  X103  +  X113  +  X123  -  X130  -  X131  -  X132  =  0 

15)  X203  +  X213  +  X223  -  X230  -  X231  -  X232  =  0 

16)  XI 10  +  XI 12  +  XI 13  -  Y101  -  Y121  -  Y131  >=0 

17)  X210  +  X212  +  X213- Y201  -  Y221  -  Y231  >=  0 

18)  X120  +  X121  +  X123  -  Y102  -  Y1 12  -  Y132  >=  0 

19)  X220  +  X221  +  X223  -  Y202  -  Y212  -  Y232  >=  0 

20)  X130  +  X131  +  X132  -  Y103  -  Y1 13  -  Y123  >=  0 

21)  X230  +  X231  +  X232  -  Y203  -  Y213  -  Y223  >=  0 

22)  X101  +  X102  +  X103  -  Y1 10  -  Y120  -  Y130  >=  0 

23)  X201  +  X202  +  X203  -  Y210  -  Y220  -  Y230  >=  0 

24)  X101  +X102  +  X103<=  1 

25)  X201  +  X202  +  X203  <=  1 

26)  XI 10  +  X120  +  X130  <=  1 

27)  X210  +  X220  +  X230  <=  1 

28)  X112  +  X121  <=  1 

29)  X212  +  X221  <=  1 

30)  X113  +X131  <=  1 

31)  X213  +  X231  <=  1 

32)  X123  +  X132  <=  1 

33)  X223  +  X232  <=  1 

34)  XI 12  +  XI 13  +X121  +X123  +  X131  +  X132  <=  2 

35)  X212  +  X213  +  X221  +  X223  +  X231  +  X232  <=  2 

36)  3.5  X101  +  4  X102  +  2.5  X103  +  3.5  XI 10  +  4.5  XI 12  +  5.5  XI 13 
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+  4  X120  +  4.5  X121  +  4.3  X123  +  2.5  X130  +  5.5  X131  +  4.3  X132 
<=  12 

37)  3.5  X201  +  4  X202  +  2.5  X203  +  3.5  X210  +  4.5  X212  +  5.5  X213 
+  4  X220  +  4.5  X221  +  4.3  X223  +  2.5  X230  +  5.5  X231  +  4.3  X232 
<=  12 

38) -  20  X101  +  X101DEL  +  X101UP  <=  0 

39) -  20  XI 02  +  X102DEL  +  X102UP  <=  0 

40) -  20  XI 03  +  X103DEL  +  XI 03 UP  <=  0 

41) -  20  XI 10  +  XI 10DEL  +  XI 10UP  <=  0 

42) -  20  XI 12  +  XI 12DEL  +  XI 12UP  <=  0 

43) -  20  XI 13  +  XI 13DEL  +  XI 13UP  <=  0 

44) -  20  X120  +  X120DEL  +  X120UP  <=  0 

45) -  20  X121  +  X121DEL  +  X121UP  <=  0 

46) -  20  XI 23  +  X123DEL  +  X123UP  <=  0 

47) -  20  XI 30  +  X130DEL  +  X130UP  <=  0 

48) - 20X131  +X131DEL  +  X131UP<=  0 

49) -  20  XI 32  +  X132DEL  +  X132UP  <=  0 

50) -  20  X201  +  X201DEL  +  X201UP  <=  0 

5 1 ) -  20  X202  +  X202DEL  +  X202UP  <=  0 

52) -  20  X203  +  X203DEL  +  X203UP  <=  0 

53) -  20  X210  +  X210DEL  +  X210UP  <=  0 

54) -  20  X212  +  X212DEL  +  X212UP  <=  0 

55) -  20  X213  +  X213DEL  +  X213UP  <=  0 

56) -  20  X220  +  X220DEL  +  X220UP  <=  0 

57) -  20  X221  +  X221DEL  +  X221UP  <=  0 

58) -  20  X223  +  X223DEL  +  X223UP  <=  0 

59) -  20  X230  +  X230DEL  +  X230UP  <=  0 

60) -  20  X23 1  +  X23 1DEL  +  X23 1UP  <=  0 

6 1 ) -  20  X232  +  X232DEL  +  X232UP  <=  0 

66) -  XOUP  -  XODEL  +  ZOO  =  0 

67) -  X1UP  -  X1DEL  +  Z1 1  =  0 

68) -  X2UP  -  X2DEL  +  Z22  =  0 

69) -  X3UP  -  X3DEL  +  Z33  =  0 

70) -  X201DEL  -  X201UP  -  X101DEL  -  X101UP  +  Z01  =  0 

71) -  X202DEL  -  X202UP  -  X102DEL  -  X102UP  +  Z02  =  0 

72) -  X203DEL  -  X203UP  -  X103DEL  -  X103UP  +  Z03  =  0 

73) -  X210DEL  -  X210UP  -  XI 10DEL  -  XI 10UP  +  Z10  =  0 

74) -  X212DEL  -  X212UP  -  XI 12DEL  -  XI 12UP  +  Z12  =  0 

75) -  X213DEL  -  X213UP  -  XI 13DEL  -  XI 13UP  +  Z13  =  0 

76) -  X220DEL  -  X220UP  -  X120DEL  -  X120UP  +  Z20  =  0 

77) -  X221DEL  -  X221UP  -  X121DEL  -  X121UP  +  Z21  =  0 

78) -  X223DEL  -  X223UP  -  X123DEL  -  X123UP  +  Z23  =  0 

79) -  X230DEL  -  X230UP  -  X130DEL  -  X130UP  +  Z30  =  0 

80) -  X231DEL  -  X231UP  -  X131DEL  -  X131UP  +  Z31  =  0 

81) -  X232DEL  -  X232UP  -  X132DEL  -  X132UP  +  Z32  =  0 
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82) 

ZOO 

+  Z03 

+ 

Z02 

+ 

Z01  <= 

17 

83) 

Zll 

+  Z13 

+ 

Z12 

+ 

Z10<= 

7 

84) 

Z22 

+  Z23 

+ 

Z21 

+ 

Z20  <= 

11 

85) 

Z33 

+  Z32 

+ 

Z31 

+ 

Z30  <= 

10 

86) 

Z30 

+  Z20 

+ 

Z10 

< 

= 

10 

87) 

Z31 

+  Z21 

+ 

Z01 

< 

= 

6 

88) 

Z32 

+  Z12 

+ 

Z02 

< 

= 

4 

89) 

Z23 

+  Z13 

+ 

Z03 

< 

- 

4 

90  zl0*Y110  +  z20*Y120  +  z30*Y130  +  z01*Y101  +  z21*Y121  +  z31*Y131  + 
z31*Y102  +  zl2*Y112  +  z32*Y132  +  z03*Y103  +  zl3*Y113  +  z23*Y123  <=  7 

91)  zl0*Y210  +  z20*Y220  +  z30*Y230  +  z01*Y201  +  z21*Y221  +  z31*Y231  + 
z02*Y202  +  zl2*Y212  +  z32*Y232  +  z03*Y203  +  zl3*Y213  +  z23*Y223  <=  7 

92) -  50  Y1 10  -  50  Y210  +  Z10  <=  0 

93) -  50  Y1 12  -  50  Y212  +  Z12  <=  0 

94) -  50  Y113  -  50  Y213  +  Z13  <=  0 

95) -  50  Y1 20  -  50  Y220  +  Z20  <=  0 

96) -  50  Y121  -  50  Y221  +  Z21  <=  0 

97) - 50  Y123  -  50  Y223  +  Z23  <=  0 

98) -  50  Y 130  -  50  Y230  +  Z30  <=  0 

99) -50Y131 -50Y231 +Z31  <=  0 

100) -  50  Y132  -  50  Y232  +  Z32  <=  0 

101) -  50  Y101  -  50  Y201  +  Z01  <=  0 

102) -  50  Y102  -  50  Y202  +  Z02  <=  0 

103) -  50  Y103  -  50  Y203  +  Z03  <=  0 

104)  50  Y1 +Z31 +Z21 +Z01  >=  0 

105) -  50  Y1  +  Z13  +  Z12  +  Z10  <=  0 

106)  50  Y 1  +  Z3 1  +  Z21  +  Z01  <=  50 

107) -  50  Y1  +  Z13  +  Z12  +  Z10  >=  -  50 

108)  50  Y2  +  Z32  +  Z12  +  Z02  >=  0 

109) -  50  Y2  +  Z23  +  Z21  +  Z20  <=  0 

110)  50  Y2  +  Z32  +  Z12  +  Z02  <=  50 

111) -  50  Y2  +  Z23  +  Z21  +Z20>=-50 

112)  50  Y3  +  Z23  +  Z13  +  Z03  >=  0 

1 13) -  50  Y3  +  Z32  +  Z31  +  Z30  <=  0 

1 14)  50  Y3  +  Z23  +  Z13  +  Z03  <=  50 

1 15) -  50  Y3  +  Z32  +  Z31  +  Z30  >=  -  50 

116)  50  YO  +  Z30  +  Z20  +  Z10  >=  0 

117) - 50  Y0  +  Z03+Z02  +  Z01  <=  0 

118)  50  YO  +  Z30  +  Z20  +  Z10  <=  50 

1 19) -  50  YO  +  Z03  +  Z02  +  Z01  >=  -  50 

120)  50Y1  +Z11  >=  7 

121) -  50Y1  +Z11  >=-50 

122)  50  Y2  +  Z22  >=  11 

123) -  50  Y2  +  Z22  >=  -  50 

124)  50  Y3  +  Z33  >=  10 
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125) -  50  Y3  +  Z33  >=  -  50 

126)  50  Y0  +  ZOO  >=  15 

127) -  50  YO  +  ZOO  >=  -  50 

128)  XODEL  +  X203DEL  +  X202DEL  +  X201DEL  +  X103DEL  +  X102DEL 
+  X101DEL<=  2 

129)  XI DEL  +  X213DEL  +  X212DEL  +  X210DEL  +  XI 13DEL  +  XI 12DEL 
+  X110DEL<=  1 

130)  X2DEL  +  X223DEL  +  X221DEL  +  X220DEL  +  X123DEL  +  X121DEL 
+  X120DEL  <=  2 

131)  X3DEL  +  X232DEL  +  X231DEL  +  X230DEL  +  X132DEL  +  X131DEL 
+  X130DEL<=  1 

132)  XOUP  +  X203UP  +  X202UP  +  X201UP  +  X103UP  +  X102UP  +  X101UP 
<=  15 

133)  X1UP  +  X213UP  +  X212UP  +  X210UP  +  XI 13UP  +  XI 12UP  +  XI 10UP 
<=  6 

134)  X2UP  +  X223UP  +  X221UP  +  X220UP  +  X123UP  +  X121UP  +  X120UP 
<=  9 

135)  X3UP  +  X232UP  +  X231UP  +  X230UP  +  X132UP  +  X131UP  +  X130UP 
<=  9 

136)  Y101  +  Y121  +  Y131  +  Y201  +  Y221  +  Y231  <=  1 

137)  Y102  +  Y112  +  Y132  +  Y202  +  Y212  +  Y232  <=  1 

138)  Y103  +  Y1 13  +  Y123  +  Y203  +  Y213  +  Y223  <=  1 

END 

INTEX101 
INTE  XI 02 
INTE  X103 
INTE  XI 10 
INTE  XI 12 
INTE  XI 13 
INTE  X120 
INTEX121 
INTE  X123 
INTE  XI 30 
INTEX131 
INTE  XI 32 
INTE  X201 
INTE  X202 
INTE  X203 
INTEX210 
INTE  X212 
INTE  X2 13 
INTE  X220 
INTE  X221 
INTEX223 
INTE  X230 
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INTEX231 
INTEX232 
INTEY101 
INTE  Y121 
INTEY131 
INTE  Y201 
INTE  Y221 
INTEY231 
INTE  Y102 
INTE  Y1 12 
INTE  Y132 
INTE  Y202 
INTEY212 
INTE  Y232 
INTE  Y103 
INTE  Y1 13 
INTE  Y123 
INTE  Y203 
INTE  Y2 13 
INTE  Y223 
INTE  Y1 10 
INTE  Y120 
INTE  Y130 
INTEY210 
INTE  Y220 
INTEY230 
INTEY1 
INTEY2 
INTE  Y3 
INTE  YO 


Appendix  11:  The  result  (1)  to  the  First  Milestone 


Local  optimal  solution  found  at  step:  442 

Objective  value:  133.5500 

Branch  count:  30 


Variable 

Value  Reduced  Cost 

X102 

1.000000 

0.0000000 

X123 

1.000000 

0.0000000 

X130 

1.000000 

0.0000000 

X201 

1.000000 

0.0000000 

X210 

1.000000 

0.0000000 

Z01 

6.000000 

0.0000000 

Z02 

4.000000 

0.0000000 

XO 

1.000000 

0.0000000 

XI 

1.000000 

0.0000000 

X2 

1.000000 

0.0000000 

X3 

1.000000 

0.0000000 

Y201 

1.000000 

0.0000000 

Y102 

1.000000 

0.0000000 

Y103 

1.000000 

0.0000000 

X102UP 

4.000000 

0.0000000 

X123UP 

0.1665335E-14 

o.ooooot 

X201UP 

4.000000 

0.0000000 

X201DEL 

2.000000 

0.0000000 

ZOO 

7.000000 

0.0000000 

XOUP 

7.000000 

0.0000000 

Zll 

7.000000 

0.0000000 

XI  DEL 

1.000000 

0.0000000 

X1UP 

6.000000 

0.0000000 

Z22 

11.00000 

0.0000000 

X2DEL 

2.000000 

0.0000000 

X2UP 

9.000000 

0.0000000 

Z33 

10.00000 

0.0000000 

X3DEL 

0.9999999 

0.0000000 

X3UP 

9.000000 

0.0000000 

Y3 

1.000000 

0.0000000 

YO 

1.000000 

0.0000000 
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Appendix  12:  The  result  (2)  to  the  First  Milestone 


Local  optimal  solution  found  at  step:  402 

Objective  value:  133.5500 

Branch  count:  23 


Variable 

Value  Reduced  Cost 

X103 

1.000000 

0.0000000 

X120 

1.000000 

0.0000000 

X132 

1.000000 

0.0000000 

X201 

1.000000 

0.0000000 

X210 

1.000000 

0.0000000 

Z01 

6.000000 

0.0000000 

Z23 

0.5551115E-15 

0.0000( 

Z32 

4.000000 

0.0000000 

xo 

1.000000 

0.0000000 

XI 

1.000000 

0.0000000 

X2 

1.000000 

0.0000000 

X3 

1.000000 

0.0000000 

Y201 

1.000000 

0.0000000 

Y132 

1.000000 

0.0000000 

Y103 

1.000000 

0.0000000 

X132UP 

3.000000 

0.0000000 

X132DEL 

1.000000 

0.0000000 

X201UP 

6.000000 

0.0000000 

ZOO 

11.00000 

0.0000000 

XODEL 

2.000000 

0.0000000 

XOUP 

9.000000 

0.0000000 

Zll 

7.000000 

0.0000000 

X1DEL 

1.000000 

0.0000000 

X1UP 

6.000000 

0.0000000 

Z22 

11.00000 

0.0000000 

X2DEL 

2.000000 

0.0000000 

X2UP 

9.000000 

0.0000000 

Y3 

1.000000 

0.0000000 

YO 

1.000000 

0.0000000 

110 


Appendix  13:  The  Formulation  of  Second  Milestone 


MIN  3X101  +3.5X102  +  2X103  +  3X110  +  4X112  +  5X113  +  3.5X120 
+  4  X121  +  3.8  X123  +  2  X130  +  5  X131  +  3.8  X132  +  3  X201 
+  3.5  X202  +  2  X203  +  3  X210  +  4  X212  +  5  X213  +  3.5  X220 
+  4  X221  +  3.8  X223  +  2  X230  +  5  X231  +  3.8  X232  +  3  X001 
+  3.5  X002  +  2  X003  +  3  X010  +  4  X012  +  5  X013  +  3.5  X020 
+  4  X021  +  3.8  X023  +  2  X030  +  5  X031  +  3.8  X032  -  4  Z32 
-  3.75  Z31  +  .5  Z30  -  2  Z23  -  3.75  Z21  +  .5  Z20  -  2  Z13  -  4  Z12 
+  .5  Z10  -  2  Z03  -  4  Z02  -  3.75  Z01 
SUBJECT  TO 

2)  X101  +  X102  +  X103  +  X201  +  X202  +  X203  =  2 

3)  XI 10  +  XI 12  +  XI 13  +  X210  +  X212  +  X213  +  X010  +  X012  +  X013  =  1 

4)  X120  +  X121  +  X123  +  X220  +  X221  +  X223  +  X020  +  X021  +  X023  =  1 

5)  X130  +  X131  +  X132  +  X230  +  X23 1  +  X232  +  X030  +  X03 1  +  X032  =  1 

6)  XI 10  +  X120  +  X130  +  X210  +  X220  +  X230  =  2 

7)  X101  +X121  +X131  +X201  +  X221  +  X231  +  X001  +  X021  +  X031  =  1 

8)  X102  +  XI 12  +  X132  +  X202  +  X212  +  X232  +  X002  +  X012  +  X032  =  1 

9)  X103  +  XI 13  +  X123  +  X203  +  X213  +  X223  +  X003  +  X013  +  X023  =  1 

10)  X101  -X110-X112-X113+X121  +X131  =  0 

11)  X201  -  X210  -  X212  -  X213  +X221  +X231  =  0 

12)  X001  -  X010  -  X012  -  X013  +  X021  +  X031  =  0 

13)  X102  +  X112-X120- X121  -  X123  +  X132=  0 

14)  X202  +  X212  -  X220  -  X221  -  X223  +  X232  =  0 

15)  X002  +  X012  -  X020  -  X021  -  X023  +  X032  =  0 

16)  X103  +  X113  +  X123  -  X130-X131  -  X132  =  0 

17)  X203  +  X213  +  X223  -  X230  -  X231  -  X232  =  0 

18)  X003  +  X013  +  X023  -  X030  -  X031  -  X032  =  0 

19)  XI 10  +  XI 12  +  XI 13  -  Y101  -  Y121  -  Y131  >=  0 

20)  X210  +  X212  +  X213  -  Y201  -  Y221  -  Y231  >=  0 

21)  X010  +  X012  +  X013  -  Y001  -  Y021  -  Y031  >=  0 

22)  X120  +  X121  +X123- Y102- Y112- Y132>=  0 

23)  X220  +  X221  +  X223  -  Y202  -  Y212  -  Y232  >=  0 

24)  X020  +  X021  +  X023  -  Y002  -  Y012  -  Y032  >=  0 

25)  X130  +  X131  +X132-  Y103  -  Y113  -  Y123  >=  0 

26)  X230  +  X23 1  +  X232  -  Y203  -  Y213  -  Y223  >=  0 

27)  X030  +  X03 1  +  X032  -  Y003  -  Y013  -  Y023  >=  0 

28)  X101  +  X102  +  X103  -  Y1 10  -  Y120  -  Y130  >=  0 

29)  X201  +  X202  +  X203  -  Y210  -  Y220  -  Y230  >=  0 

30)  X001  +  X002  +  X003  -  Y010  -  Y020  -  Y030  >=  0 

31)  X101  +X102 +  X103  <=  1 

32)  X201  +  X202  +  X203  <=  1 

33)  X001  +  X002  +  X003  <=  1 

34)  XI 10  +  X120  +  X130  <=  1 
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35)  X210  +  X220  +  X230  <=  1 

36)  X010  +  X020  +  X030  <=  1 

37)  X112  +  X121  <=  1 

38)  X212  +  X221  <=  1 

39)  X012  +  X021  <=  1 

40)  X113+X131  <=  1 

41)  X213  +  X231  <=  1 

42)  X013+X031  <=  1 

43)  X123  +X132  <=  1 

44)  X223  +  X232  <=  1 

45)  X023  +  X032  <=  1 

46)  XI 12  +  XI 13  +X121  +X123  +X131  +X132  <=  2 

47)  X212  +  X213  +  X221  +  X223  +  X231  +  X232  <=  2 

48)  X012  +  X013  +  X021  +  X023  +  X031  +  X032  <=  2 

49)  Z03  +  Z02  +  Z01  <=  17 

50)  Z13  +  Z12  +  Z10  <=  7 

51)  Z23  +  Z21 +Z20<=  11 

52)  Z32  +  Z31 +Z30<=  10 

53)  Z30  +  Z20  +  Z10  <=  10 

54)  Z31  +Z21  +Z01  <=  6 

55)  Z32  +  Z12  +  Z02  <=  4 

56)  Z23  +  Z13  +Z03  <=  4 

57)  Z01  Y001  +  Z21  Y021  +  Z31  Y031  +  Z02  Y002  +  Z12  Y012  +  Z32  Y032  +  Z03 
Y003  +  Z13  Y013  +  Z23  Y023  +  Z10  Y010  +  Z20  Y020  +  Z30  Y030  <  0; 

58)  Z01  Y101  +Z21  Y121  +Z31  Y131  +Z02  Y102  +  Z12  Y112  +  Z32  Y132  +  Z03 
Y103  +  Z13  Y1 13  +  Z23  Y123  +  Z10  Y1 10  +  Z20  Y120  +  Z30  Y130  <=  7 

59)  Z01  Y201  +  Z21  Y221  +  Z31  Y231  +  Z02  Y202  +  Z12  Y212  +  Z32  Y232  +  Z03 
Y203  +  Z13  Y213  +  Z23  Y223  +  Z10  Y210  +  Z20  Y220  +  Z30  Y230  <=  7 

64) -  50  Y1 10  -  50  Y210  +  Z10  <=  0 

65) -  50  Y1 12  -  50  Y212  +  Z12  <=  0 

66) -  50  Y1 13  -  50  Y213  +  Z13  <=  0 

67) -  50  Y120  -  50  Y220  +  Z20  <=  0 

68) -  50  Y121  -  50  Y221  +  Z21  <=  0 

69) -  50  Y123  -  50  Y223  +  Z23  <=  0 

70) -  50  Y1 30  -  50  Y230  +  Z30  <=  0 

71) -  50  Y131  -  50  Y231  +  Z31  <=  0 

72) -  50  Y 132  -  50  Y232  +  Z32  <=  0 

73) -50Y101 -50Y201 +Z01  <=  0 

74) -  50  Y 102  -  50  Y202  +  Z02  <=  0 

75) -  50  Y103  -  50  Y203  +  Z03  <=  0 

76)  Y101  +  Y121  +  Y131  +  Y201  +  Y221  +  Y231  +  Y001  +  Y021  +  Y031  =  1 

77)  Y102  +  Y1 12  +  Y132  +  Y202  +  Y212  +  Y232  +  Y002  +  Y012  +  Y032  =  1 

78)  Y103  +  Y1 13  +  Y123  +  Y203  +  Y213  +  Y223  +  Y003  +  Y013  +  Y023  =  1 

79)  Y1 10  +  Y120  +  Y130  +  Y210  +  Y220  +  Y230  +  Y010  +  Y020  +  Y030  =  1 

END 
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INTEX101 
INTE  XI 02 
INTEX103 
INTE  XI 10 
INTE  XI 12 
INTE  XI 13 
INTE  XI 20 
INTEX121 
INTE  XI 23 
INTE  XI 30 
INTEX131 
INTE  XI 32 
INTE  X201 
INTE  X202 
INTEX203 
INTEX210 
INTEX212 
INTE  X213 
INTE  X220 
INTEX221 
INTE  X223 
INTE  X230 
INTEX231 
INTE  X232 
INTE  X001 
INTEX002 
INTE  X003 
INTEX010 
INTE  X012 
INTEX013 
INTE  X020 
INTEX021 
INTE  X023 
INTEX030 
INTEX031 
INTEX032 
INTE  Y101 
INTEY121 
1NTEY131 
INTE  Y201 
INTEY221 
INTEY231 
INTE  Y001 
INTEY021 
INTEY031 


INTE  Y102 
INTE  Y112 
INTE  Y132 
INTE  Y202 
INTE  Y212 
INTE  Y232 
INTE  Y002 
INTE  Y012 
INTE  Y032 
INTE  Y103 
INTE  Y1 13 
INTE  Y123 
INTE  Y203 
INTEY213 
INTE  Y223 
INTE  Y003 
INTE  Y013 
INTE  Y023 
INTE  Y1 10 
INTE  Y120 
INTE  Y130 
INTE  Y210 
INTE  Y220 
INTE  Y230 
INTE  Y010 
INTE  Y020 
INTE  Y030 


Appendix  14:  The  Solution  to  the  Second  milestone  (without  generalized  Benders’ 
Decomposition) 


Local  optimal  solution  found  at  step:  3625 

Objective  value:  127.5500 

Branch  count:  304 


Variable 

X101 

X110 

X202 

X223 

X230 

Z01 

Z12 

Z13 

XO 

XI 

X2 

X3 

Y101 

Y212 

Y213 

Y230 


Value 

Reduced  Cost 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

6.000000 

0.0000000 

4.000000 

0.0000000 

3.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 

1.000000 

0.0000000 
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Appendix  15:  The  formulation  of  TSP1  for  Second  milestone 


MIN  3.8  X032  +  5  X031  +  2  X030  +  3.8  X023  +  4  X021  +  3.5  X020  +  5  X013 
+  4  X012  +  3  X010  +  2  X003  +  3.5  X002  +  3  X001  +  3.8  X232 
+  5  X231  +  2  X230  +  3.8  X223  +  4  X221  +  3.5  X220  +  5  X213 
+  4  X212  +  3  X210  +  2  X203  +  3.5  X202  +  3  X201  +  3.8  X132 
+  5  X131  +  2  X130  +  3.8  X123  +  4  X121  +  3.5  X120  +  5  XI 13 
+  4  XI 12  +  3  XI 10  +  2  X103  +  3.5  X102  +  3  X101 
SUBJECT  TO 

2)  X203  +  X202  +  X201  +  X103  +  X102  +  X101  =  2 

3)  X013  +  X012  +  X010  +  X213  +  X212  +  X210  +  XI 13  +  XI 12  +  XI 10  =  1 

4)  X023  +  X021  +  X020  +  X223  +  X221  +  X220  +  X123  +  X121  +  X120  =  1 

5)  X032  +  X031  +  X030  +  X232  +  X23 1  +X230  +  X132  +  X131  +X130=  1 

6)  X230  +  X220  +  X210  +  X130  +  X120  +  XI 10  =  2 

7)  X031  +X021  +X001  +X231  +X221  +X201  +X131  +X121  +X101  =  1 

8)  X032  +  X012  +  X002  +  X232  +  X212  +  X202  +  X132  +  XI 12  +  X102  =  1 

9)  X023  +  X013  +  X003  +  X223  +  X213  +X203  +X123  +  X113  +  X103  =  1 

10)  X131  +X121  -XI 13  -X112-X110  +  X101  =  0 

11)  X231  +  X221  -  X213  -  X212  -  X210  +  X201  =  0 

12)  X031  +  X021  -  X013  -  X012  -  X010  +  X001  =  0 

13)  X132  -  X123  -  X121  -  X120  +  XI 12  +  X102  =  0 

14)  X232  -  X223  -  X221  -  X220  +  X212  +  X202  =  0 

15)  X032  -  X023  -  X021  -  X020  +  X012  +  X002  =  0 

16) -  X132  -  X131  -  X130  +  X123  +  X113  +  X103  =  0 

17) -  X232  -  X23 1  -  X230  +  X223  +  X2 1 3  +  X203  =  0 

1 8) -  X032  -  X03 1  -  X030  +  X023  +  XO 1 3  +  X003  =  0 

19)  XI 13  +  XI 12  +  XI 10  >=  0 

20)  X213  +  X212  +  X210  >=  0 

21)  X013  +  X012  +  X010  >=  1 

22)  X123+X121  +X120>=  0 

23)  X223  +  X221  +  X220  >=  1 

24)  X023  +  X021  +  X020  >=  0 

25)  X132 +  X131  +X130  >=  1 

26)  X232  +  X23 1  +  X230  >=  0 

27)  X032  +  X03 1  +  X030  >=  0 

28)  X103  +  X102  +  X101  >=  1 

29)  X203  +  X202  +  X201  >=  0 

30)  X003  +  X002  +  X001  >=  0 

31)  X103 +X102 +  X101  <=  1 

32)  X203  +  X202  +  X201  <=  1 

33)  X003  +  X002  +  X001  <=  1 

34)  X130  +  X120  +  XI 10  <=  1 

35)  X230  +  X220  +  X2 10  <=  1 

36)  X003  +  X002  +  X001  <=  1 
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37)  X121  +X112  <=  1 

38)  X221  +X212<=  1 

39)  X021  +X012  <=  1 

40)  X131  +X113  <=  1 

41)  X231  +X213  <=  1 

42)  X031  +X013  <=  1 

43)  X132  +  X123<=  1 

44)  X232  +  X223  <=  1 

45)  X032  +  X023  <=  1 

46)  X132  +  X131  +X123  +  X121  +X113  +  X112<=  2 

47)  X232  +  X231  +  X223  +  X221  +  X213  +  X212  <=  2 

48)  X032  +  X031  +  X023  +  X021  +  X013  +  X012  <=  2 
END 
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Appendix  16:  The  solution  to  the  TSP  1  for  Second  milestone 


Global  optimal  solution  found  at  step: 
Objective  value:  17.00000 


Variable 

X103 

X130 

X202 

X220 

X001 

X010 

Y001 

Y202 

Y113 

Y130 


Value  Reduced  Cost 
1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 

1.000000  0.0000000 
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Appendix  17:  The  formulation  of  Inventory/Allocation  for  Second  Milestone 

MODEL : 

MIN  = 

-  3.75*z01  -  4*z02  -  2*z03  +  0.5*zl0  -  4*zl2  -  2*zl3 
+  0.5*z20  -  3.75*z21  -  2*z23  +  0.5*z30  -  3.75*z31  -  4*z32 
+  48.75*x0  +  1 8*x  1  +  5 1 ! *x2  +  39*x3  ; 

zOl  +  z02  +  z03  <=  17; 
zlO  +  zl2  +  zl3  <=  7; 
z20  +  z21  +  z23  <=  11; 
z30  +  z31  +  z32  <=  10; 

zlO  +  z20  +  z30  <=  10; 
zOl  +z21  +z31  <=  6; 
z02  +  zl2  +  z32  <=  4; 
z03  +  zl3  +  z23  <=  4; 

zl0*y010  +  z20*y020  +  z30*y030  + 
z01*y001  +z21*y021  +  z31*y031  + 
z02*y002  +  zl2*y012  +  z32*y032  + 
z03*y003  +  zl3*y013  +  z23*y023  < 0; 

zlO*yl 10  +  z20*yl20  +  z30*yl30  + 
z01*yl01  +z21*yl21  +z31*yl31  + 
z02*yl02  +  zl2*yl 12  +  z32*yl32  + 
z03*yl03  +  zl3*yl 13  +  z23*yl23  <  7; 

zl0*y210  +  z20*y220  +  z30*y230  + 
z01*y201  +z21*y221  +z31*y231  + 
z02*y202  +  zl2*y212  +  z32*y232  + 
z03*y203 +  zl3*y213  +  z23*y223  <  7; 

zlO  -  50*yl  10  -  50*y210  <  0; 
zl2  -  50*yll2  -  50*y212  <  0; 
zl3  -  50*yl  13  -  50*y213  <  0; 

z20  -  50*yl20  -  50*y220  <  0; 
z21  -  50*yl21  -  50*y221  <  0; 
z23  -  50*yl23  -  50*y223  <  0; 

z30  -  50*yl30  -  50*y230  <0; 
z31  -  50*yl31  -  50*y231  <0; 
z32  -  50*yl32  -  50*y232  <0; 
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zOl  -  50*yl01  -  50*y201  <0; 
z02  -  50*yl02  -  50*y202  <  0; 
z03  -  50*yl03  -  50*y203  <  0; 


xO  =  1; 
xl  =  1; 
x2=  1; 
x3  =  1; 

ylOl  =  0 
yl21  =  0 
yl31  =  0 
y201  =  0 
y221  =  0 
y23 1=0 
yOOl  =  1 
y021  =  0 
y03 1=0 

yl02  =  0 
yll2  =  0 
yl32  =  0 
y202  =  1 
y212  =  0 
y232  =  0 
y002  =  0 
y012  =  0 
y032  =  0 

yl03  =  0 
yl  13=1 
yl23  =  0 
y203  =  0 
y213  =  0 
y223  =  0 
y003  =  0 
y013  =0 
y023  =  0 

yll0  =  0 
yl20  =  0 
yl30  =  1 
y210  =  0 
y220  =  0 
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y230  =  0; 
yOlO  =  0; 
y020 = 0; 
y030  =  0; 
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Appendix  18:  The  solution  to  the  Inventory/Allocation  for  Second  Milestone 

Global  optimal  solution  found  at  step:  6 

Objective  value:  132.7500 


Variable 

Value 

Reduced  Cost 

Z02 

4.000000 

0.0000000 

Z13 

4.000000 

0.0000000 

XO 

1.000000 

0.0000000 

XI 

1.000000 

0.0000000 

X2 

1.000000 

0.0000000 

X3 

1.000000 

0.0000000 

Y001 

1.000000 

0.0000000 

Y130 

1.000000 

0.0000000 

Y113 

1.000000 

0.0000000 

Y202 

1.000000 

0.0000000 

Row 

Slack  or  Surplus 

Dual  Price 

8  0.0000000  4.000000 

9  0.0000000  2.000000 

10  0.0000000  3.750000 

13  0.0000000  0.0000000 
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Appendix  19:  The  formulation  of  Master  Problem  for  Second  Milestone 

MIN  Z 
SUBJECT  TO 

2) - 100.5  Y001  -  100.5  Y021  -  100.5  Y031  - 101  Y002  -  101  Y012 

-  101  Y032  -  92  Y003  -  92  Y013  -  92  Y023  -  48.75  Y101 

-  48.75  Y121  -  48.75  Y131  -  35  Y102  -  35  Y1 12  -  35  Y132  -  33  Y103 

-  33  Y113  -  33  Y123  -  48.75  Y201  -  48.75  Y221  -  48.75  Y231 

-  38.75  Y202  -  38.75  Y212  -  38.75  Y232  -  31  Y203  -  31  Y213 
-31  Y223  +  Z  >=  25.5 

3) -  101  Y002  -  101  Y012  - 101  Y032  -  92  Y003  -  92  Y013  -  92  Y023 

-  35  Y102  -  35  Y1 12  -  35  Y132  -  33  Y103  -  33  Y1 13  -  33  Y123 

-  38.5  Y202  -  38.5  Y212  -  38.5  Y232  -  31  Y203  -  31  Y213  -  31  Y223 

-  39  Y010  -  39  Y020  -  39  Y030  -  18  Y110  - 18  Y120  - 18  Y130 
- 18  Y210- 18  Y220-  18  Y230  +  Z>=  25.5 

4) -  97.5  Y001  -  97.5  Y021  -  97.5  Y031  -  92  Y003  -  92  Y013  -  92  Y023 

-  48.75  Y101  -  48.75  Y121  -  48.75  Y131  -  33  Y103  -  33  Y1 13 

-  33  Y123  -  52.25  Y201  -  52.25  Y221  -  52.25  Y231  -  31  Y203 

-  31  Y213  -  31  Y223  -  39  Y010  -  39  Y020  -  39  Y030  -  18  Y110 

- 18  Y120  -  18  Y130  -  18  Y210  - 18  Y220  -  18  Y230  +  Z  >=  25.5 

5) -  97.5  Y001  -  97.5  Y021  -  97.5  Y031  -  102  Y002  -  102  Y012  - 102  Y032 

-  48.75  Y101  -  48.75  Y121  -  48.75  Y131  -  53  Y102  -  53  Y1 12 

-  53  Y132  -  52.25  Y201  -  52.25  Y221  -  52.25  Y231  -  51  Y202 

-  51  Y212  -  51  Y232  -  39  Y010  -  39  Y020  -  39  Y030  -  18  Y110 

-  18  Y120  -  18  Y130  -  18  Y210  -  18  Y220  -  18  Y230  +  Z  >=  25.5 

6)  4Y002  +  4  Y013<=  0 

7)  4Y102  +  4  Y113<=  7 

8)  4Y202  +  4  Y213<=  7 

9)  Y001  +  Y021  +  Y031  +  Y101  +  Y121  +  Y131  +  Y201  +  Y221  +  Y231  =  1 

10)  Y002  +  Y012  +  Y032  +  Y102  +  Y1 12  +  Y132  +  Y202  +  Y212  +  Y232  =  1 

1 1)  Y003  +  Y013  +  Y023  +  Y103  +  Y1 13  +  Y123  +  Y203  +  Y213  +  Y223  =  1 

12)  Y010  +  Y020  +  Y110  +  Y120  +  2  Y130  + Y210  + Y220  +  Y230  =  1 

13)  Y001  +  Y021  +  Y031  +  Y002  +  Y012  +  Y032  +  Y003  +  Y013  +  Y023 
+  Y010  +  Y020  +  Y130  <=  2 

14)  Y101  +Y121  +  Y131  +Y102  +  Y112  +  Y132  +  Y103  +  Y113  +  Y123 
+  Y120  +  Y130  +  Y210  <=  2 

15)  Y201  +  Y221  +  Y231  +  Y202  +  Y212  +  Y232  +  Y203  +  Y213  +  Y223 
+  Y210  +  Y220  +  Y230  <=  2 

END 

INTE  Y001 
INTEY021 
INTEY031 
INTE  Y002 
INTE  Y012 
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INTE  Y032 
INTE  Y003 
INTE  YOB 
INTE  Y023 
INTEY101 
INTEY121 
INTE  Y1 31 
INTE  Y102 
INTE  Y1 12 
INTE  Y132 
INTE  Y 103 
INTE  Y1 13 
INTE  Y123 
INTE  Y201 
INTE  Y221 
INTEY231 
INTE  Y202 
INTE  Y212 
INTE  Y232 
INTE  Y203 
INTE  Y2 13 
INTE  Y223 
INTE  YO 10 
INTE  Y020 
INTE  Y030 
INTE  Y1 10 
INTE  Y120 
INTE  Y130 
INTE  Y2 10 
INTE  Y220 
INTE  Y230 


Appendix  20:  The  solution  to  the  Master  Problem  for  Second  Milestone 


Global  optimal  solution  found  at  step: 
Objective  value:  144.0000 

Branch  count:  3 


Variable 

Z 

Y121 

Y232 

Y223 

Y110 


Value  Reduced  Cost 


144.0000 

1.000000 

1.000000 

1.000000 

1.000000 


0.0000000 

48.75000 

38.75000 
31.00000 
0.0000000 
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